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ABSTRACT 


A  comprehensive  study  of  the  problem  of  helicopter  noise  radiation  is  presented. 

A.ter  a  review  of  the  basic  features  of  the  noise,  the  limited  experimental  data  are 
reviewed  in  some  detail,  and  empirical  laws  are  proposed.  An  exact  theoretical 
expression  for  the  noise  radiation  is  derived.  This  expression  has  been  used  as  the 
bnsis  for  the  development  of  a  comprehensive  computer  program  to  calculate  heli¬ 
copter  noise  at  any  field  point,  including  all  effects  of  fluctuating  airloads  and 
rigid  and  flexible  blade  motions.  Details  of  this  program  are  presented  in  a  com¬ 
panion  report.  Under  very  reasonable  approximations,  an  analytic  expression  has 
been  found  for  the  sound  field  far  from  the  helicopter.  Computations  based  on  this 
expression  have  been  made.  The  results  show  that  it  is  the  very  high  harmonics 
of  the  loading  which  contribute  to  the  important  harmonics  of  the  sound  field.  For 
instance,  calculation  of  the  tenth  harmonic  cf  a  four-blade  rotor  requires  a 
knowledge  of  loading  harmonics  up  to  the  sixtieth.  Details  of  such  loadings  are 
not  available  from  theory  or  experiment.  Therefore,  rotor  aerodynamic  loadings 
have  been  reviewed  in  detail,  and  empirical  harmonic  decay  laws  have  been  de¬ 
rived.  Loading  phases  appear  to  be  best  described  as  random,  and  this  introduces 
simplification  in  the  theory,  together  with  the  necessity  for  definition  of  a  correla¬ 
tion  length .  Results  of  a  parameter  study  show  trends  basically  in  agreement  with 
experiment,  with  sound  at  the  higher  harmonics  basically  proportional  to  thrust  times 
disc  loading  times  tip  velocity  squared  .  For  the  lower  harmonics,  the  dependence 
on  tip  velocity  is  to  the  2B  power  where  B  is  the  number  of  blades.  The  effect  of 
forward  speed  is  to  increase  the  sound  radiated  forward  and  to  decrease  that 
radiated  aft,  causing  a  difference  between  fore  and  aft  radiation  of  as  much  as 
20  dB  for  the  second  and  third  harmonics  at  a  forward  Mach  number  of  0.25.  An 
effective  rotational  Mach  number  concept  is  introduced  which  enables  the  effects 
of  forward  speed  to  be  calculated,  to  good  accuracy,  directly  from  results  for  the 
hover  case.  The  overall  sound  directionality  pattern  is  found  to  have  a  minimum 
slightly  above  the  plane  of  the  disc  and  a  broad  maximum  about  20  degrees  below. 
The  effects  of  both  the  near  field  and  the  blade  motion  are  found  to  be  small .  Appen¬ 
dices  present  analyses  of  the  blade  motion  effects  and  of  the  noise  radiation  by  ran¬ 
dom  blade  loadings.  Appendix  III  presents  design  charts  for  the  noise  radiation  based 
or  the  theory  that  reduces  the  calculation  of  noise  to  a  simple  algebraic  procedure. 
The  theory  generally  shows  .fair  agreement  with  experiment  for  overall  levels  and 
good  agreement  for  trends.  The  charts  presented  should  therefore  be  of  direct  use 
fer  design  trade-off  studies. 
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The  work  reported  herein  was  performed  by  Wyle  Laboratories  Research  Staff,  under 
Contract  DAAJ02-67-C-0023,  for  U.S.  Army  Aviation  Materiel  Laboratories, 

Forr  Cusiis,  Virginia.  The  work  was  carried  out  under  the  technical  cognizance  of 
Mr.  William  E.  Nettles  of  the  USAAVLABS  staff.  A  companion  report,  USAAVLABS 
TR  58-61 ,  gives  details  of  the  computer  program  specially  developed  for  helicopter 
noise  calculations  under  this  contract. 

Wyle  Laboratories  personnel  directly  associated  with  the  project  included  Dr.  M.V, 
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1.0 


INTRODUCTION 


The  present  work  is  a  further  step  in  the  U.S.  Army  Aviation  Materiel  Laboratories* 
sponsored  research  into  the  causes  of,  and  means  for  controlling,  helicopter  noise. 

This  overall  research  program  has  now  been  in  progress  for  several  years,  since  it 
first  oecame  apparent  that  noise  radiation  was  on  important  factor  in  helicopter  opera¬ 
tions .  both  military  and  civil .  Several  studies  of  helicopter  noise  have  now  been 
performed  both  from  the  experimental  ^~'7  and  the  theoretical  8-12  points  of  view. 
However,  it  is  unfortunately  still  true  that  the  basic  noise  problem  is  far  from  being 
solved,  or  even  understood. 

This  eport  concentrates  on  the  problems  caused  by  the  far-field  noise  radiation  of 
the  helicopter.  Aural  detection  of  a  helicopter  is  a  significant  factor  in  the  success 
of  tactical  missions.^  The  unique  characteristics  of  helicopter  noise  not  only  cause 
the  noise  to  propagate  for  exceptionally  great  distances  but  also  probably  enable  an 
experienced  listener  to  distinguish  the  type  of  aircraft,  its  range,  and  its  direction  of 
flight.  In  the  commercial  sphere,  high  noise  levels  threaten  the  success  of  the  very 
operations  which  a  helicopter  is  uniquely  able  to  perform;  that  is,  operations  out  of 
terminals  in  built-up  and  often  heavily  populated  areas  which  are  inaccessible  to 
conventional  aircraft.  There  are  already  many  instances  where  the  disadvantages  of 
noise  have  been  found  to  outweigh  other  advantages  of  helicopter  travel .  However, 
it  should  be  pointed  out  that  internal  noise  problems  ere  gIso  severe.  In  addition  to 
the  <  bvious  problem  of  interference  with  communications,  noise  levels  inside  the 
crew  compartment  are  sufficiently  high  to  constitute  a  very  definite  health  hazard 
to  helicopter  ai 'crews.  Figure  1  is  a  1/3-octave  band  analysis  of  the  noise  in  the 
cockpit  of  a  CH-47B  in  cruising  flight.^  Superimposed  on  this  spectrum  are  hazar¬ 
dous  noise  exposure  curves^  and  the  amount  of  sound  attenuation  which  can  be 
expected  from  well-designed  and  well-fitting  helmet/earmuff  combinations.  It  is 
obvious  that  the  noise  level  is  dangerously  high.  Nevertheless,  internal  noise  prob¬ 
lems  will  not  be  specifically  studied  in  this  report,  although  many  of  the  methods 
developed  and  the  results  obtained  are  relevant  to  both  the  internal  and  the  external 
noise  field . 

Helicopter  noise  can  be  divided  into  two  basic  groups:  that  arising  aerodynamical ly 
and  that  arising  mechanically.  However,  with  the  exception  of  piston  engine  exhaust 
noise,  the  "mechanical"  sound  sources,  including  those  due  to  the  gearbox,  transmis¬ 
sion  and  the  various  vibrating  components,  are  important  to  the  internal  and  near- 
exte  nal  fields  only.  In  the  far-external  field,  the  aerodynamical  ly  generated  sound 
is  dominant.  The  latter  is  associated  with  .the  rotor  airloads  and  includes  various 
types  of  noise  which  are  commonly  known  as  rotational  noise,  vortex  noise,  and 
"blade  slap" .  Which  of  these  is  most  significant  depends  on  a  number  of  factors,  the 
most  important  of  which  are  the  location  of  the  observer  and  the  flight  and  configura¬ 
tion  of  the  helicopter.  At  moderate  distances  from  the  helicopter,  the  various  sources, 
listed  in  their  order  of  importance  to  the  subjectively  judged  magnitude  cf  the  sound, 
or  loudness,  are:0 


•  Blade  slap  (when  it  occurs) 

•  Piston  engine  exhaust  noise 

•  Tail  rotor  rotational  noise 

•  Main  rotor  vortex  noise 

•  Main  rotor  rotational  noise 

•  Gearbox  noise 

•  Turbine  engine  noise 

•  Other  sources 

At  very  sho  t  distances  inside  the  helicopter,  gearbox  noise  becomes  more  important, 
as  noted  pre  viously.  At  extreme  distances,  approaching  the  range  at  which  the  sound 
is. barely  audible,  blade  slap  and  main  rotor  rotational  noise  may  be  the  only  compo¬ 
nents  which  are  heard,  because  atmospheric  and  other  sound  absorption  effects  re¬ 
move  the  high-frequency  energy  from  the  spectrum.  Main  rotor  rotational  noise  con¬ 
tains  the  major  part  of  its  energy  at  low  and  suboudible  frequencies. 

1  ]  DEFINITION  of  rotor  noise  sources 

The  basic  subject  of  the  present  study  is  therefore  the  noise  radiated  by  the  rotor. 

This  consists  of  three  of  the  sources  mentioned  above:  blade  slap,  rotational  noise, 
and  vortex  noise.  A  brief  review  of  the  characteristics  of  each  is  given  below.  How¬ 
ever,  the  underlying  source  of  the  noise  radiated  by  a  helicopter  rotor  is  the  fluctu¬ 
ating  forces  upon  it,  and  the  basic  cause  of  these  forces  may  be  found  by  examination 
of  the  rotor  aerodynamics.  Thus,  in  order  to  understand  the  noise,  it  is  important  first 
of  all  to  understand  the  rotor  aerodynamics. 

The  main  feature  of  rotor  aerodynamics  is  the  lack  of  symmetry.  In  forward  flight  the 
advancing  Hade  encounters  substantially  higher  air  velocities  than  does  the  retreating 
blade,  giving  rise  to  cyclic  variations  in  the  resulting  airloads.  The  equilibrium  of 
the  rotor  is  maintained  by  a  combination  of  cyclic  pitch  control  and  the  freedom  of 
the  rotor  bhdes  to  flap. 

However,  the  most  important  feature  of  the  helicopter  aerodynamics  from  the  stand¬ 
point  of  noi ;e  and  vibration  is  the  rotor  wake.  Each  blade  acts  in  the  same  way  as  a 
wing  in  flight,  and  the  lift  on  it  generates  a  vortex  wake  behind  it  which  has  a 
strong  tendency  to  roll  up  into  a  concentrated  vortex  core.  Each  blade  must  there¬ 
fore  pass  over  the  concentrated  vortex  wake  left  by  its  predecessor.  Depending  on 
advance  rat  o,  net  lift  force,  and  so  on,  this  vortex  may  pass  either  extremely  close 
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to,  or  far  away  from,  fhe  following  blade.  If  the  vortex  passes  close  to  the  blade, 
then  a  substantial  local  increase  in  lift  will  occur  temporarily.  These  comparatively 
rapid  increments  in  lift,  caused  by  vortex  interaction,  are  very  efficient  noise  radia¬ 
tors.  It  appears  that  a  very  large  part  of  the  observed  noise  from  a  helicopter  can  be 
attributed  to  these  vortex  effects. 

Blade  Slap 

A  particularly  severe  type  of  noise  is  well  known  in  helicopter  operation  under  the 
name  "blade  slap"  .  If  is  found  that  under  various  conditions  (for  instance,’  during 
low  power  descent),  the  helicopter  produces  a  particularly  loud  slapping  or  banging 
noise,  which  occurs  at  the  blade  passage  frequency.  This  occurs  at  precisely  those 
conditions  where  fhe  vortex  wake  can  be  expected  to  pass  very  close  to  the  rotor , 
and  can  be  particularly  severe  on  a  tandem-rotor  aircraft  where  the  wake  from  the 
first  rotor  can  pass  through  the  second.  Blade  slap  can  also  occur  on  high  speed 
rotors  and  in  this  case  is  associated  with  transonic  flow  over  the  rotor  blades.  Thus 
there  are  two  possible  sources  of  blade  slap,  due  to  either  vortex  interaction  or  tran¬ 
sonic  flow .  It  may  be  noted  that  the  description  of  the  phenomenon  is  fairly 
straightforward  acoustically,  and  both  the  effects  mentioned  above  can  be  readily 
predicted  from  the  theory.  Cases  corresponding  to  blade  slap  are  discussed  in 
Section  6.1  of  this  report.  Reference  should  also  be  made  to  the  studies  by  Leverton 
and  Taylor.^ 

!« 

When  it  occurs,  blade  slap  is  by  far  the  loudest  source  of  noise  observed  on  the 
helicopter.  However,  it  seems  inappropriate  to  consider  blade  slap  as  a  separate 
phenomenon.  The  helicopter  rotor  is  always  undergoing  some  form  of  vortex  inter¬ 
action,  and  blade  slap  is  simply  a  particularly  severe  form.  Perhaps  it  is  more 
realistic  to  suppose  that,  at  least  from  the  acoustic  point  of  view,  the  helicopter  is 
always  flying  under  some  degree  of  blade  slap. 

Rotational  Noise 

Rotational  noise  is  usually  regarded  as  that  sound  which  is  directly  attributable  to  the 
steady  and  fluctuating  lift  and  drag  forces  acting  on  the  rotor.  This  is  the  major 
source  of  noise  from  propellers  15,16,17  an(j  helicopter  tail  rotors,  where  the  funda¬ 
mental  frequency  is  of  the  order  of  100  Hz.*  For  the  helicopter  main  rotor,  the  fun¬ 
damental  frequency  is  significantly  lower,  around  10  Hz,  due  to  the  lower  rotational 
speed.  Rotational  noise  has  a  spectrum  consisting  of  a  number  of  very  narrow  peaks 
occurring  at  integral  multiples  of  ..e  blade  passage  frequency.  However,  its  impor¬ 
tance  is  often  underestimated  because  of  its  very  low  frequency  and  the  fact  that 
most  of  its  energy  is  contained  within  "subaudible"  frequencies.  The  quotation  marks 
here  ere  used  because  the  subajdible  frequency  range  is  in  fact  difficult  to  define, 


Herrz  (Hz)  is  the  new  international  unit  for  frequency.  1  Hertz  =  1  cycle  per  second 
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since  whet  cannot  be  heard  can  often  be  felt  through  other  physiological  sensations. 

The  subjective  character  of  rotational  noise  is  highly  dependent  upon  its  harmonic 
content.  The  relctive  amplitudes  of  the  higher  harmonics  vary  from  helicopter  to 
nelicopte  and  strongly  depend  upon  the  flight  condition.  If  the  harmonic  ampli¬ 
tudes  fall  off  rapidly  (with  harmonic  number),  the  sound  can  best  be  described  as  a 
"thumpin  |"  which  occurs  at  the  blade  passage  frequency.  As  the  harmonic  decay 
becomes  ess.  that  is,  '<<  the  higher  harmonics  increase  in  amplitude,  the  thumps 
sharpen  into  bangs  and  eventuaiiy  into  "blade  slap"  .  An  example  of  this  effect  is 
the  sound  of  a  helicopter  flying  overhead  at  a  foirlv  high  speed  and  low  altitude.  As 
the  mach  ne  approaches,  the  characteristic  sound  is  a  sharp  popping  which  decays  to 
a  thumping  sound  as  the  aircraft  flies  overhead.  Subsequently,  the  rotational  content 
becomes  barely  perceptible  as  the  aircraft  recedes,  with  the  tail  rotor  noise  domi¬ 
nating  th  ;  observed  sound. 

A  sound  component  which  is  very  similar  in  nature  to  rotational  noise  is  termed 
"thicknes  noise"  J'7  This  arises  due  to  the  displacement  of  air  particles  as  a  blade 
passes  through  their  volume.  Like  rotational  noise,  thickness  noise  has  a  frequency 
spectrum  consisting  of  a  number  of  discrete  harmonics  of  the  fundamental  blade  pas¬ 
sage  frequency,  which  is  the  rotational  frequency  multiplied  by  the  number  of  blades. 
However,  this 's  rara!y  an  important  source  of  helicopter  noise  except  at  very  high 
tip  speed,.  Thickness  noise  has  not  been  considered  in  this  report. 

Vortex  Noise 

Vortex  noise  is  the  name  given  to  the  distinctive  "swishing”  sound  which  character¬ 
izes  helicopter  noise  at  short  distances.  It  is  sometimes  thought  to  be  caused  by  the 
turbulence  associated  with  the  blade  boundary  layers,  but  it  is  probably  principally 
due  to  the  interaction  of  the  blades  with  wake  turbulence.  Vortex  noise  is  random 
in  nature  and  contains  sound  energy  which  is  spread  over  a  substantial  portion  of  the 
audible  frequency  range.  Rotor  vortex  noise  differs  from  other,  more  familiar,  forms 
of  random  noise,  such  as  wind  noise  or  jet  noise,  in  that  its  amplitude  and  frequen¬ 
cies  are  modulated.  This  is  due  to  the  varying  relative  distance  and  ‘'elocity  between 
the  source  and  the  observer  because  of  blade  rotation,  which  causes  harmonic  varia¬ 
tions  in  both  amplitude  and  frequency. 

Experimental  data  presented  in  this  report  indicate  tha*,  contrary  to  previous  opinions, 
rotations.  I  noise  may  continue  to  dominate  the  sound  spectrum  up  to  frequencies  of 
•400  Hz  end  higher.  The  pre/ious  belief  that  vortex  noise  becomes  more  significant 
at  frequencies  greater  than  abou  IX  Hz  may  well  be  the  result  of  two  factors:  (1)  one- 
third-octave  analysis  of  helicopter  noise  does  not  distinguish  individual  harmonics 
,at  the  higher  frequencies,  and  (2)  previous  theoretical  results  predicted  that  rota¬ 
tional  noise  decays  mare  rapidly  than  it  does  in  practice.  It  is  likely  that  improved 
narrow  band  analysis  of  helicopter  noise  will  show  that  rotational  noise  is  signifi¬ 
cant  through  the  major  part  of  the  audible  frequency  range  and,  thus,  that  previously 
identified  "vortex"  noise  is  really  part  of  the  rotational  noise  output  of  the  helicopter. 
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ft  is  clear  from  the  preceding  discussion  that  all  significant  forms  of  rotor  noise  are 
due  to  the  fluctuations  and  motions  of  the  pressures  acting  on  the  blades.  For  this 
reason,  it  is  somewhat  difficult,  and  indeed  of  doubtful  justification,  to  separate 
rotational  noise,  vortex  noise,  and  blade  slap.  It  is  shown  in  this  report  that  ob¬ 
server  rotational  noise- harmonics,  even  at  relatively  low  frequencies,  are  most 
likeiy  explained  by  random  phase  variations  in  the  higher  blade  loading  harmonics. 

Main  otor  rotational  noise  harmonics  as  high  as  the  thirty-fifth  are  discernable  in 
experimental  data  for  two— blade  rotors.  To  obtain  theoretical  predictions  of  these 
harmonic  amplitudes  with  any  accuracy  requires  knowledge  of  at  least  100  harmonics 
of  the  rotor  airloads.  Such  fluctuations  correspond  to  pressure  variations  occurring 
within  rhe dimensions  ofa  blade  chord;  a  condition  approaching  the  domain  of  "vortex" 
noise  generation .  Similarly,  blade  slap  is  said  to  occur  when  the  observed  sound 
becomes  highly  impulsive  in  nature,  which  in  turn  is  the  result  of  the  blade  loading 
ha  ,-mo  lies'  reaching  some  critical  level.  It  could  be  argued  that  "blade  slap"  is 
alway.:  present  in  helicopter  noise;  it  is  merely  a  question  of  degree. 

The  fundamental  object  of  the  report  is  to  predict  all  possible  forms  of  the  noise  radiated 
by  the  helicopter  main  rotor.  This  is  achieved  by  applying  the  basic  acoustical  equa¬ 
tions  which  give  the  sound  radiation  from  a  known  fluctuating  force  distribution .  These 
basic  equations  were  recently  derived  In  a  convenient  form  by  Lawson,  '  Thus,  the 
problem  of  predicting  the  noise  radiation  reduces  to  the  problem  of  predicting  the 
rotor  dynamic  loads.  Unfortunately,  the  accurate  calculation  of  the  noise  field  re¬ 
quires  an  extremely  detailed  knowledge  of  the  fluctuating  loads.  For  instance,  it 
will  be  shown  that  calculation  of  the  tenth  sound  harmonic  for  a  four-blade  rotor  (a 
case  v.  hich  is  certainly  of  practical  interest)  requires  a  knowledge  of  the  loadings  up 
to  about  the  sixtieth.  This  knowledge  is  not  at  present  available,  either  from  experi¬ 
ment  or  theory.  Indeed,  the  only  loading  harmonics  which  can  at  present  be  specified 
will  be  shown  to  have  a  negligible  sound  output.  Thus,  in  order  to  calculate  the 
sound  field,  these  high  harmonic  loadings  must  be  predicted,  and  an  important  part  of 
the  present  work  has  been  to  derive  semiempirical  predictions  for  these. 

Because  of  the  complexity  of  the  acoustic  problem,  a  complete  solution  is  possible 
only  via  a  digital  computer.  A  key  secondary  objective  is  therefore  to  reduce  com¬ 
puter  time  to  a  minimum.  The  previous  theoretical  studies  have  been  accom¬ 

plished  using  this  general  approach,  and  reasonable,  although  far  from  complete, 
agreement  with  experiment  was  found.  Computing  time  using  these  approaches  was 
of  the  order  of  minutes  per  field  point.  The  present  work  uses  a  different  basic 
expression  for  the  computer  studies,  which  removes  one  of  the  principal  problems 
that  cause  excessive  computer  time.  Computer  time  for  the  present,  more  complete 
solution  has  been  reduced  to  about  10  seconds  per  field  point,  and  this  solution 
includes  all  flapping  and  lagging  effects  as  desired.  In  addition,  under  very  reason¬ 
able  approximations,  an  analytic  solution  has  been  obtained.  The  solution  is  a 
rathe*  complex  collection  of  Bessel  functions  and  still  requires  computer  calcula¬ 
tion,  but  computer  time  for  this  case  can  be  reduced  to  the  order  of  a  tenth  of  a 
secord  per  field  point.  Furthermore,  this  basic  analytic  solution  has  revealed  several 
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important  features  of  the  noise  generation  process  which  are  of  considerable  help  in 
understanding  the  problem. 

Thus,  two  very  different  approaches,  the  exact  numerical  method  and  the  approximate 
analytic  solution,  have  been  used  in  the  work.  The  two  approaches  give  identical 
answers  if  appropriate,  identical  input  conditions  are  used.  This  self-checking  feature 
of  the  present  calculations  is  important  in  establishing  confidence  in  the  results.  Fair 
•  agreement  with  experiment  is  also  found,  although  experimental  data  are  far  from 
being  either  complete  or  trustworthy.  Previous  theoretical  studies  have  shown  that 
.  rotational  noise  can  be  predicted  with  reasonable  accuracy  at  the  lower  frequencies. 

.  The  present  work  Illustrates  some  of  the  deficiencies  of  earlier  knowledge  and  shows 
how  theory  and  experiment  complement  each  other  to  yield  methods  for  the  accurate 
‘  estimation  of  rotational  noise  throughout  its  important  frequency  range. 
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2.0 


REVIEW  OF  THE  PROBLEM 


2.  1  REVIEW  OF  EXPERIMENTAL  DATA 

One  of  the  major  difficulties  in  establishing  a  proper  understanding  of  helicopter  noise 
is  the  iack  of  good  experimental  data.  Several  studies  have  been  performed,'  ~ ^ ,  but 
all  5*.  ffer  from  deficiencies.  Few  investigators  have  performed  a  sufficiently  detailed 
study  of  the  problem,  and  even  the  detailed  studies  have  often  given  data  which  are  of 
very  limited  application.  Furthermore,  comparison  of  results  from  different  investiga¬ 
tions  on  the  same  helicopter  often  reveals  discrepancies  of  10  dB  or  more.  Never¬ 
theless,  in  order  to  provide  a  basis  for  the  report,  it  is  felt  that  a  review  of  the  data  is 
necessary.  Attempts  have  been  made  in  each  case  to  determine  the  leading  trends, 
and  empirical  formulas  for  noise  characteristics  are  given.  It  should  be  noted  that  all 
the  formulas  given  apply  under  restricted  circumstances,  and  considerable  care  should 
be  exercised  in  any  attempt  to  apply  the  formulas  to  any  real  prediction  case. 

From  many  points  of  view,  the  first  reported  study,  by  Hubbard  and  Maglieri,^  still 
gives  the  most  useful  information  on  the  rotor  noise  problem.  They  measured  noise 
from  a  helicopter  rotor  on  a  test  tower.  The  tower  was  over  a  diameter  high,  so  that 
ground  interference  should  not  be  too  severe.  Noise  was  measured  at  only  two  loca¬ 
tions:  at  the  hub  and  on  the  ground  at  1  diameter  from  the  hub,  but  the  rotor  was 
run  over  a  very  wide  range  of  disc  loadings  ar.d  tip  speeds.  Although  these  ranges 
are  outside  the  useful  practical  limits,  the  results  do  show  important  trends,  which 
would  have  been  difficult  to  pinpoint  with  any  confidence  in  a  more  restricted  exper¬ 
iment.  A  plot  of  the  overall  levels,  as  measured  on  the  ground  I  diamete<  away 
from  therotor,  is  given  in  Figure  2.  The  data  should  be  fairly  realistic,  although  there 
seem,  to  be  a  possibility  of  downwash-induced noise  at  the  microphone  for  the  high 
disc  ioading  cases.  There  are  several  features  that  can  be  observed  on  this  graph. 

First  it  can  be  observed  that  noise  levels  rise  rapidly  with  tip  speed .  However,  for 
any  oecific  tip  velocily,  the  noise  level  reaches  a  minimum  at  an  intermediate  value 
of  disc  loading.  The.  same  effect  is  shown  in  the  one-third  octave  band  spectra  in 
Figu  e  3,  also  taken  from  Hubbard  and  Maglieri's  paper.  At  very  high  disc  loadings 
the  ! ofor  stalls,  and  this  causes  major  increments  in  noise.  Points  corresponding  to 
stall  ;d  operation  are  marked  on  Figure  2.  At  very  low  disc  loadings  the  noise  also 
rises.  This  is  due  to  the  vortex  wake  interaction  effects  discussed  In  the  Introduction. 
Inde  -d ,  Hubbard  and  Maglieri  appear  to  have  been  the  first  to  report  the  blade  slap 
phenomenon  .  Note  the  increase  in  sound  output  at  low  disc  loadings  for  the  900  ft/sec 
cast.  This  corresponds  to  the.  blade  slap  condition. 

The  data  shown  in  Figure  2  may  be  replotted  in  a  form  suitable  for  making  predictions 
of  overall  noise  level.  It  will  be  assumed  that  helicopters  will  not,  in  general,  run 
in  c  editions  close  to  stall  or  blade  slap.  When  such  extreme  points  are  removed  from 
Figure  2,  the  data  may  be  replotted  as  shown  in  Figure  4 .  This  figure  probably  applies  for 
collective  pitches  of  the  order  of  8  degrees .  Variation  of  collective  pitch  is,  of  course. 
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reflected  in  the  thrust  term  .  Jtwill  be  observed  that  fairly  good  collapse  is  obtained  .  The 
doubi'  logarithmic  plot  of  Figure  4  is  convenient  for  establishing  power  law  trends. 

It  wil  be  observed  that  the  sound  level  rises  at  a  much  faster  rate  with  velocity  at 
high  i"p  speeds  than  at  low.  This  is  a  turther  important  feature  of  the  noise.  Most 
helicopters  operate  with  tip  speeds  in  the  range  of  from  500  to  800  fps.  Over  this 
range  lhe  aata  collapse  to  a  straight  line  fairly  well.  From  this,  an  empirical  predic¬ 
tion  lew  for  overall  levels  at  a  500-ft  sideline  may  be  obtained  as 
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where 

i 

T  =  rotor  thrust  (lift) 

Correcrion  of  Equation  (1)  to  other  ranges  r  may  be  readily  accomplished  by 


dBr  =  dB5oo  +  20  lo9,0  (500/r)  (2) 

It  shou;d  be  pointed  out  that,  in  the  derivation  of  Equation  (I)  from  the  data  in 
Figure  4,  t wo  assumptions  have  been  made.  First,  it  was  assumed  that,  for 
constant  disc  loading,  the  sound  output  was  proportional  to  rotor  diameter  squared. 
This  se;ms  physically  reasonable  and  was  found  to  apply  theoretically  in  the  related 
case  of  compressor  noise.  Second,  it  was  assumed  that  Hubbard  and  Maglieri's 
results  could  be  corrected  to  500  ft,  assuming  spherical  spreading,as  in  Equation  (2). 
This  assumption  is  piubobly  accurate  to  within  a  dB  or  so,  as  will  be  shown  later 
in  the  report.  Furthermore,  Hubbard  and  Maglieri  do  not  give  one-third  octave  plots 
below  1 00  Hz,  end  there  is  therefore  some  doubt  as  to  the  lower  frequency  limit  in 
their  overall  response  data.  It  might  also  be  pointed  out  that  there  are  some  theoret¬ 
ical  re -'sons  that  suggest  a  thrust-squared  law  in  Equation  1),  so  that  the  thrust  term 
would  Men  be  20  log,Q  T..  However,  the  thrust  term  shown  here  has  a  good  empir¬ 
ical  fit  ;o  the  limited  data  available,  as  is  shown  by  the  collapse  in  Figure  4,  It 
should  be  noted  that  configuration  changes  may  well  be  important.  Equation  (2) 
applies  essentially  to  a  two-blade  case.  Data  from  Cox  and  Lynn  2  suggest  that  a 
three -b  ade  refer  may  be  about  3  dB  quieter  than  a  two-blade  rotor  with  the  same 
thrust  and  tip  speed. 

Since  thrust  is  proportional  to  tip  velocity  squared,  Equation  (1)  implies  a  velocity  to 
the  eighth  power  law  for  the  noise.  Acousticians  should  note  that  this  variation  is  not 
too  surp  ising.  The  sixth  power  law  for  dipole  noise  applies  only  to  isolated  dipole 
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f  luck  citing  forces.  Distributions  of  force  (for  instance  over  the  rotor  disc)  imply 
higher  order  multipoles  and  associated  higher  order  velocity  powers.  Indeed,  a  simple 
propeller  obeys  a  velocity  to  the  tenth  power  law  over  most  of  its  operating  range  of 
interest . 


Stuckey  and  Goddard'  recently  presented  new  acoustic  data  on  a  rotor  tower  test. 
Their  50-ft-diameter  rotor  was  only  20  feet  from  the  ground,  and  vortex 
reingustioti  with  consequent  increases  in  noise  level  is  to  be  expected.  Theirdataare 
unlikely  to  be  representative  of  hover  cut  of  ground  effect,  but  the  data  may  beoppli- 
cable  to  lift-off  maneuvers .  They  give  fairly  complete  data,  which  justify  analysis.  No 
simple  effect  of  noise  harmonics  is  observed.  The  harmonic  level  varies  almost 
randomly  from  harmonic  to  harmonic  and  from  test  to  test.  This  effect  is  probably  due, 
at  least  in  part,  to  the  vortex  reingestion  effects  mentioned  above.  However, some 
trends  can  be  seen.  Figure  5  gives  harmonic  levels  measured  for  a  high. and  low 
collective  pitch  at  high  and  low  rotational  speeds.  Figure  5  shows  a  fairly  constant 
10-df  increase r  i  going  from  the  low  to  the  high  thrust  condition.  Note  that  spectral 
shapes  in  each  case  are  changed  little.  However,  there  is  a  very  definite  change  in 
spectral  shape  between  the  high  and  the  low  tipspeed  case  .  The  lower  harmonics  are 
considerably  increased.  At  low  tipspeeds,  the  sixth  to  thirteenth  harmonics  are  almost 
5  dE  higher  than  the  first;  at  high  tip  speeds,  the  effect  is  reversed,  with  the  first 
harmonic  about  7  dB  higher. 


As  was  mentioned  above,  considerable  scatter  is  present  in  Stuckey  and  Goddard's  re¬ 
sults  However,  in  order  to  establish  some  general  trends,  the  graphs  shown  in  Figure  6 
have  been  prepared  .  The  highest  harmonic  given  in  the  data  of  Stuckey  and  Goddard , 
the  fifteenth ,  is  taken  for  reference  .  It  was  found  that  the  sound  level  obeyed  a  V2  T2 
lav/  when  V  is  tip  velocity  and  T  is  rotor  thrust.  Figure  6a  shows  the  sound  level 
against  thrust,  corrected  for  velocity.  It  is  clear  that  a  T2  law  gives  a  good  fit  to 
the  data.  Figure  6b  gives  the  sound  level  against  velocity,  corrected  for  thrust. 

Here  ,  considerable  scatter  may  be  observed,  and  the  V2  lew  cannot  be  regarded  as 
definitely  established.  It  is  of  particular  interest  to  note  that  several  points  at  about 
11  t  ;  13  degrees  collective  pitch  are  about  5  to  10  JB  low.  This  probably 
corresponds  to  a  cleaner  flow  condition,  as  was  observed  in  Hubbard  and  Maglieri's 
data,  and  offers  some  hope  for  a  noise  control  method.  If  the  V2T2  law  is  accepted, 
then,  since  T  ~  V2,  an  overall  V6  law  can  be  inferred.  The  V6  law  is,  of  course, 
the  classical  dipole  power  dependency  expected  for  a  force. ^,9 


Figuie  7  shows  the  equivalent  graphs  for  the  first  harmonic.  Far  more  scatter  Is 
apparent  here.  The  velocity  law  is  more  like  V8  or  greater,  as  opposed  to  the  V2 
law  tor  the  fifreenth  harmonic.  Part  of  this  is  undoubtedly  due  to  the  increase  in 
microphone  response  as  the  fundamental  becomes  of  higher  frequency  with  increase  in 
tip  speed.  This  accounts  for  a  V2  increase,  so  that  the  observed  V8  law  implies  an 
actual  V6  law  .  It  may  be  noted  that  use  of  the  assymptotic  form  of  the  Besse!  functions 
suggests  a  V2®  law  at  the  fundamental  frequency  where  B  is  the  number  of  blades. 
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Stuckey  anc:  Goddard's  data  were  taken  on  a  three-blade  rotor,  and  a  V6  law  is  then 
suggested  theoretical ly,  in  broad  agreement  with  the  experimental  results.  It  might 
also  be  noted  that  the  same  arguments,  including  blade  number  and  microphone  re¬ 
sponse,  applied  to  Hubbard  and  Maglieri's  two-blade  rotor  would  suggest  a  V6  T2  law 
for  the  funcamenta!  frequency.  A  V6  law  is  suggested  in  Equation  (1),  but  probably 
little  weight  should  be  placed  on  this.  The  data  plotted  against  thrust  in  Figure  7a 
again  suggest  a  T2  law.  Considerable  scatter  occurs  here,  basically  due  to  the 
scatter  in  the  velocity  cose.  For  any  fixed  tip  velocity,  the  data  correspond  quite 
closely  to  the  T2  law.  No  detailed  prediction  equations  will  be  given  based  on 
these  data,  but  the  general  trends  should  be  noted,  particularly  the  fact  that  the  sig¬ 
nificance  of  the  lower  harmonics  increases  with  tip  speed. 

2 

For  the  complete  helicopter,  there  is  again  little  data.  Sternfeld  et  al,  gave  an 
extremely  hrge  amount  of  acoustic  data  on  many  different  vehicles.  However,  their 
program  was  primarily  aimed  at  internal  noise  levels,  and  most  of  their  reported  data 
are  not  in  c  convenient  form  for  analysis  in  the  present  problem.  The  most  useful  study 
for  the  present  work  was  that  of  Cox  and  Lynn. 3  They  performed  a  detailed  investi¬ 
gation  of  the  noise  from  a  UH-1A  helicopter.  Typical  overall  levels  in  their  tests  are 
also  shew?  in  Figure  4.  The  extremely  good  agreement  between  Hubbard  and  Maglieri's 
results  for  cn  isolated  rotor  and  Cox  and  Lynn's  results  for  the  complete  helicopter 
should  be  noted.  This  reinforces  the  contention  that  the  main  rotor  is  the  principal 
source  of  sound  and  is  support  for  the  use  of  Equation  (I)  for  prediction  purposes.  Un¬ 
fortunately  no  other  data  have  been  found  in  a  convenient  form  for  inclusion  in 
Figure  4,  so  that  the  use  of  Equation  (1)  cannot  be  regarded  as  completely  justified. 
However,  I  may  be  of  some  use  in  noise  predictions  where  parameters  do  not  vary  too 
much  from  those  of  the  source  data . 

Cox  and  Ly  in  also  pe  rformed  some  narrow  band  spectral  analyses  of  their  acoustic 
signal .  However,  an  error  was  apparently  mode  in  the  identification  of  the  funda¬ 
mental  frequency  as  1 1  Hz  in  their  original  graph.  Use  of  a  10-Hz  fundamental  as 
shown  lead  to  much  better  agreement  with  observed  frequency  plots.  Figure  8  is 
based  on  their  results  with  this  correction,  and  it  suggests  that  the  noise  radiation 
from  the  complete  helicopter  can  be  divided  into  four  parts: 

•  Main  rotor  rotational  noise 

•  Tail  rotor  rotational  noise 

•  Main  rotor  vortex  noise 

•  Gearbox  noise 

Levels  of  these  sources,  broadly  as  suggested  by  Cox  and  Lynn,  are  shown  in  Figure  8. 
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The  test  vehicle  was  turbine  powered,  and  this  probably  accounts  for  most  of  the  radia¬ 
tion  above  2000  Hz.  Results  for  a  reciprocating  engine  helicopter  would  show  addi¬ 
tional  noise  peaks  at  the  engine  firing  frequency  (about  100  Hz)  and  its  harmonics  at 
abou;  the  same  level  as  the  main  and  tail  rotor  harmonics  shown.  It  can  be  seen  that 
the  harmonics  associated  with  the  main  rotor  rotational  noise  predominate  at  lower 
frequencies  (<  100  Hz).  At  the  rnidfrequency  range  (100  to  400  Hz),  the  tail  rotor 
harmonics  appear  more  significantly,  and  the  main  rotor  "vortex”  noise  is  important 
at  around  500  Hz.  It  should  be  noted  that  this  particular  order  of  significance  is  not 
necessarily  true  for  all  helicopters.  For  instance,  the  UH-1A  operates  with  a  par¬ 
ticularly  high  blade  tip  speed.  As  was  shewn  by  Stuckey  and  Goddard's  data  (Figure 
5),  reduction  of  tip  speed  will  markedly  reduce  the  levels  of  the  first  few  main  and 
tail  rotor  harmonics,  but  only  reduce  higher  harmonics  and  vortex  noise  to  a  lesser 
degree . 

This  'vortex"  noise  justifies  further  study.  It  is  generally  considered  to  be  defined  by 
the  background  noise  showing  in  the  spectrum  below  the  peaks  due  to  the  tail  rotor 
harmonics,  as  is  shown  in  Figure  8.  Several  attempts  have  been  made  to  predict  this 
vortex  raise  by  empirical  techniques  based  on  the  assumption  that  it  is  entirely  ran¬ 
dom  or  broad  band  in  nature  (for  example  4,5, 12) #  jq  fQct#  ;t  seems  very  likely  that 
the  vortex"  noise  component  as  generally  defined  really  represents  some  of  the 
high  or  harmonics  of  the  main  rotor  rotational  noise.  Figure  9  gives  a  recent  uftra- 
narre  w  (2  Hz)  spectral  analysis  of  noise  from  a  UH-1B  helicopter  under  lift-off  con¬ 
ditions  recorded  at  Wyle  Laboratories.^  Considerable  care  was  taken  to  ensure  that 
the  spectral  analysis  was  as  accurate  as  possible.  To  produce  this  one  figure,  approxi¬ 
mately  2  hours  of  analysis  time  was  required.  Such  detail  has  not  previously  been 
given  in  the  reported  data,  and  Figure  9  shows  several  new  and  interesting  effects. 

The  most  significant  effect  is  that  harmonics  of  the  main  rotor  blade  passage  frequency 
are  :dentifiable  up  to  at  least  400  Hz  or  about  the  36th  harmonic  of  the  noise.  In 
othe  words,  the  background  noise  level  beneath  the  tail  rotor  harmonics  is  due  to  the 
higher  harmonics  of  the  main  rotor,  at  least  from  100  to  400  Hz.  Figure  9  shows  how 
the  (rue  broad-band  background  level  is  about  3  to  5  dB  below  what  is  usually  des- 
crib‘  d  as  "vortex"  noise.  The  small  subharmonic  halfway  between  the  principal  blade 
passage  harmonics  is  also  of  interest.  The  UH-1  has  two  blades,  and  this  subharmonic 
is  clearly  due  to  incomplete  cancellation  of  the  noise  radiated  by  each  blade.  Iden- 
tica  radiation  from  each  blade  is  always  assumed  theoretically,  and  Figure  9  shows 
how  rhis  is  not  necessarily  true. 

Several  authors  have  attempted  to  predict  rotor  noise  levels  above  the  tenth  harmonic. 

As  v  as  shown  above,  sound  at  these  frequencies  contains  contributions  from  both  the 
rote  harmonics  and  the  vortex  noise.  Both  sources  have  been  combined  by  pre¬ 
vious  investigators  as  "vortex"  noise,  and  equations  for  its  sound  contribution 
given .  Davidson  and  Hargest^  quote  results  from  Goddard  and  Stuckey  to  suggest 


11 


an  appro:-  .mate  formula  for  the  vortex  noise,  which  can  be  written  as 

dBi00  --  20  log, 0  Vj  +  20  log]r  T  -  10  log  S  -  25.5  (3) 

t  * 

•where  S  -  total  blade  plan  area 
12 

Schlegel  et  al,  also  give  a  formula  for  vortex  noise  which  can  be  rewritten  as 

dBsoo  “•  20  lo9| o  VT  +  20  1o9|0T  "  10  1  °9 1 o  S  "  43  (4) 


The  fact  that  the  functional  dependence  of  each  empirical  form  is  the  same  is  en¬ 
couraging,  but  it  will  be  observed  that  Davidson  and  Hargest's  equation  (3)  predicts 
17.5  dB  more  noise  than  Schlegel's  equation  (4).  Dovidson  and  Hargest's  result  is 
intended  to  apply  immediately  underneath  the  rotor,  and  they  suggest  that  much  less 
vortex  ncise  is  radiated  sideways.  Their  correction  for  this  effect  would  reduce  their 
prediction  by  about  10  dB  for  the  position  used  by  Schlegel  et  al,  in  their  tests.  The 
argument,  and  data,  put  forward  by  Davidson  and  Hargest  cannot  be  said  to  offer 
any  real  support  for  this  correction,  but  the  correction  does  appear  to  be  included  in 
their  predictions,  so  that  it  must  be  removed  for  comparison  with  Schlegel .  Several 
other  effects  suggested  by  Davidson  and  Hargest  would  tend  to  increase  the  predicted 
noise  level  still  further.  The  divergence  between  the  two  equations  is  not  readily 
explainable.  Possibly,  the  increased  number  of  blades  (five  and  six)  used  by  Schlegel 
et  al,  contributed  to  the  reduced  noise.  Empirically,  perhaps  5  dB  should  be  added 
to  Equation  (4)  for  prediction  purposes.  Stuckey  and  Goddard^  recently  presented 
further  data  on  a  rotor  tower  test,  some  of  which  were  discussed  above.  For  the  "vor¬ 
tex"  noise  effects,  they  find  essentially  the  same  velocity  dependence  as  other  in¬ 
vestigate,  s  but  suggest  a  1  .66  power  lav\'  for  the  thrust  variation.  However,  these 
data  are  thought  to  be  rather  ??  nited  in  application  because  of  the  ground  interference 
effects  noted  previously. 

In  order  o  interpret  nelicopter  noise  data,  it  is  of  extreme  importance  to  consider  the 
comparative  significance  of  the  various  sound  sources  from  the  detectability  point  of 
view.  The  significance  of  helicopter  noise  lies  in  irs  effect  on  a  human  observer.  It 
is  well  known  that  the  ear  does  not  respond  to  all  sounds  equally.  Indeed,  sounds 
which  are  very  lew  (  <  ICO  Hz)  or  very  high  (  >10,000  Hz)  in  frequency  are  not 
heard  wei  I .  Figure  10  gives  a  curve  of  equal  loudness  against  frequency.  In  other 
words,  stands  which  seem  to  have  the  same  loudness  will  actually  have  comparative 
measured  intensities  which  lie  on  the  curve  shown  in  Figure  10. 

The  meatjred  sound  pressure  levels  of  Figure  9  have  been  replotted  in  Figure  Ha. 
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Using  fiie  weighting  curve  of  Figure  10  on  these  data  gives  Figure  11b.  This  figure 
gives  the  apparent  significance  of  the  sound  from  the  observer's  point  of  view.  Note 
that  while  Figure  11a  suggests  that  the  main  rotor  is  the  dominant  source  of  noise, 

Figure  1  lb  suggests  that  from  the  point  of  view  of  a  human  observer,  the  tail  rotor 
can  be  important.  Furthermore,  the  peak  main  rotor  sound  level  occurs  at  the  second 
harmonic  in  Figure  11a.  The  first  harmonic  was  below  the  range  of  adequate  micro¬ 
phone  response  for  these  data.  Figure  11a  therefore  suggests  that  either  the,  first  or 
the  second  harmonic  is  the  largest.  However,  reference  to  Figure  lib  shows  that  the 
principal  rotor  harmonic,  as  far  cs  the  human  observer  is  concerned,  is  the  third. 

Appl  icafion  of  the  same  correction  to  the  Stuckey  and  Goddard  results  suggest  that 
the  ten  h  to  fifteenth  harmonics  were  the  most  important.  In  the  remainder  of  the 
report,  results  will  usually  be  given  in  terms  of  their  absolute  intensities;  but  it  will 
be  important  to  bear  in  mind  the  distorting  effect  of  the  human  ear  in  evaluating 
these  results  from  the  point  of  view  of  detectability. 

A  further  important  effect  is  the  attenuation  of  sound  during  its  propagation  over  long 
distances.  The  subject  is  too  complicated  to  be  studied  in  detail  here,  since  it 
depend:  on  many  parameters,  such  as  the  height  of  tt.e  helicopter  above  the  ground, 
atmospheric  humidity,  type  of  ground  cover,  and  so  on.  Discussions  may  be  found  in 
References  7,  13,  19,  and  20,  However,  the  key  point  is  that  all  the  high-ifrequency 
noise  will  be  rapidly  attenuated,  typically  at  more  than  30  dB/1000  ft  for  frequencies 
above  1000  Hz.  Thus,  from  the  detectability  viewpoint,  there  is  little  significance 
in  any  noise  above  1000  Hz.  For  instance,  the  gearbox  is  often  the  dominant  source 
inside  the  helicopter^  but  will  never  be  significant  for  detectability. 

Finally,  in  any  discussion  of  the  significance  of  various  sources  of  noise, it  is  important 
to  beat  in  mind  the  limitations  imposed  by  the  method  of  presenting  the  data .  It  is 
usual  to  present  helicopter  noise  data  in  the  form  of  narrow  band  spectra,  as  was  shown 
in  Figures  5  and  6.  This  is  reasonable  because  of  the  many  discrete  frequencies 
present,  which  can  be  identified  only  in  such  a  plot.  Consider  the  spectrum  sketched 
below,  which  shows  a  discrete  frequency  spike  at  100  Hz  with  a  lower  broad-band 
level  extending  between  95  and  105  Hz.  It  seems  obvious  that  the  discrete  fre¬ 
quency  peak  is  far  more  important.  However,  although  the  discrete  frequency  spike 

is  10  dB  higher  than  the  broad-band 
level,  the  broad-band  level  extends 
over  a  wider  range,  10  times  wider 
in  fact.  Thus,  the  discrete  frequency 
spike  and  the  broad-band  level  as 
shown  actually  have  equal  contribu- 
I  tions  to  the  overall  power,  and  they  could 

95  100  105  be  regarded  as  having  equal  importance. 

Hz  Unfortunately,  there  is  no  automatic 

way  of  measuring  the  comparative 


dB 

CJ  i  j  —  a~bwaaa aI 
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c.onttibytipos  of  various  sources  on  on  overall  power  basis.  Only  the  ©vercii  contri¬ 
bution -from  t  the -cornbinfid  sources  car  be  measured .  Comparative  contributions  ccn 
b«  calculated  only  by  c  tedious  odd  it  ion  of  the  various  terms  observable  in  narrow- 
barvd.pfeto,  .such  ss  in  figure  5,  These  caicuiati.cn*  m*vs  beer,  .mode  for  *b«':  spectrum 
and  ore  shown  in  Figure  lie. 

— N^e  mevery differ cjitvpi>**i*5=~  a?  Fijjjre  lie  compared  to  Figure  tic.  Figure 
1  jo.  suggests  the?  teti!  rotor  hormcnict  are  predominant  over  me  jTaquency  luaye  of 
-ftomblOO  tp AOQ.ttt;  byLFigore  11c  shews  that  both  main  rotor  rotational  noise  ana 
'moin'rotOf;  ”nrortex"-fvolve  -are -equally  irnpcriSrtf  -This  again  sktsss  the  significance 
of  Ithe  /naip  rpto?  ct  o  npis.e  producer .  If  the  subjective  correction  carve  of  Figure 
10  l  s  .nc*y  applied  to  Figure  lie,  the  result  is  Figure  f  Id  .  Figure  1  Id  gives  yet 
.  another  order  of  magnitude  for  the  noise  field.  It  is  thought  that,  of  the  many  pos¬ 
sible  ways  of  considering  the  data,  the  one-third  octave  plot  (subjectively  corrected), 
os  shown -in  figure  lid,  is  probably  the  most  meaningful  from  the  detectability  view¬ 
point,  "However,  for  the  bos*c  reporting  of  acoustic  data,  the  narrow  -band  analysis 
(uncorreated  fpr  subjective  effects)  of  Figures  8  end  9  is  the  most  helpful . 

In  conclusion  of  this  brief  review  of  availdble  data,  it  is  worthwhile  to  gi  .-e  two 
-further  erjtpiri cal  equations  which  -:e!ote  to  engine  .noise.  A  result  from  Davidson 
ood  Horgest's  pQper^-jnoy  be  quoted  to  give-  the  noise  radiation  of  a  piston  engine  as 

/  '  ,V  -  dB^  *  37.5  -  17.7  logJC  HP  (5) 

—>*Her«  tJF  isthe  instoHed  horsepower.  For  on  axiai  flow  jet  engine  compressor,  a 
result  from  Reference  18  may  be  used. 

dB  =  5Q  togtc  V-  t  20  ;og(C  D  -  S6  (6) 

where  Yy  -  the  com  pressor  tip  velocity  and  D  ~  the  duct  diameter. 

2 ,2  REVIEW  OF  THEORY 

Initiol  attempts  to  predict  helicopter  rotor  noise  utilized  the  large  body  of  analysis 
built  up  to  describe  the  noise  rodiotion  from  propellers,  for  example  References  15  to 
17.  If  was  inevitably  found  that,  although  propeller  noise  theory  was  sometimes  fairly 
accurate  for  the  first  harmonic  of  the  helicopter  rotational  rtoise,  it  was  cyossly  in 
error--  by  aver  100  dB  in  many  cases  -  for  the  higher  harmonics.  A  typical  example 
is  shown  in  Figure  12a,  which  will  be  discussed  in  more  detail  below.  In  fact, 
propeller  noise  theory  was  found  to  be  unacceptable  for  predicting  the  higher  harmonic 
sound  rodiotion  evert  for  propellers. 

An  davious  possible  reason  for  the  discrepancy  between  this  theory  and  helicopter 
experiments  was  the  existence  of  large  fluctuating  forces  on  the  blades  in  rhe  hell- 


cd pter  cose,  era  th<$  was  suggested  by  many  w<x«fJ .  As  p  (esul*  of  two  defatted 
thecreticcl  studies,  by  ScMege! }?  Loewy  ana  Sutton,  l  me  pqt«mt.i®l  a*gm-. ,  ^ 
ficance  of  the  fluctuating  forces  is  now  cieor.  If  It  t he  opinion  of  the  authors of  ’ 
thfc  report  that  oil  toe  significant  higher  harmonic  sound  effects  .{except  possibly  at 
transonic  or  supersonic  speeds)  coo  be  ufrrlbuted  to  .‘mwe  u*t»i*ddy 

The  two  studies  mentioned  above  justify  further  review  .  A  comparison  or  thc1f|^ptp.tio 
-  f  ,  ..  _T  „  . _  T ,  —  i  icv’i *C  T.^  Ft  12  and  13 .  Both— 

LUt  v.t  PV  C^K«5r  »u«v»h;w»  •  w«V‘»*  •-  o  1  \r*  *_  „  ^ _  .  . 

investigations  use  sabstanHoMy  the  some  approach.  The  equations  so t  sound jgeopr^^ 
lion  from  a  point  force  ore  written,  one  expressions  for  the  ,adicfion  from  the  com¬ 
plete  rotor  are  fours  by  InregrcrHon  over  in-  rjnuf  drsC  .  A  d. ffor  sc-co  :r.  fqnr.  - 

between  the  basic  equations  used  results  from  the  use  of  the  Garnet  and  Watkins  ° 
moving-axis  form  by  Loewy  and  Sutton^  ond  the  mere  usual  fixed— Qxis  form  by 
Schfegei  er  o\  J.n  each  approach,  the  necessary  inte^als  are  evaluated  on  o  com¬ 
puter .  Both  cpprooches  retain  the  acoustic  near -field  terms  in  the  po.nf  source  radia¬ 
tion.  (By  wey  of  explanation,  a  fluctuating  point  force  produces  an  acoustic  pressure 
field  which  contains  two  components,  one  of  which  tails  off  as  r2  and  one  as  r,  where  . 
r  is  distance .  Clearly,  sufficiently  far  away  from  the  source,  only  the  last  (acoustic 
far -field)  term  is  significant.  For  calculations  near  the  source  (say  a  wavelength  or 
sot,  the  first  (acoustic  near -field)  term  must  be  retained.  Schlegei  er  at  s  approach  ^ 
does  assume  a  second  ;  geometric"  far-field  approximation,  whose  terms  of  order  (R/r) 
wnere  R  is  rotor  roaius,  can  be  neglected,  thus  simplifying  the  integration.  All 
fer -field  approximations  will  be  valid  sufficiently  for  from  ;.ie  helicopter-  Schleget 
et  ai  use  a  rectangular  distribution  approximation  to  the  chordwbe  loading  pattern, 
while  loewy  and  Sutton  use  an  analytic  approximation  which  leads  to  complex  alge¬ 
braic  expressions.  This  difference  is  thought  to  be  insignificant .  The  effect  of 
chordwise  loading  is  discussed  In  detail  in  Section  5.2. 

Sehlegei  et  al^  give  aetciled  comparisons  of  their  theory  and  experiments  only  for 
the  first  four  harmonics  of  the  noise.  Fair  agreement  is  found  for  the  first  two,  but  it 
is  clear  that  underestimation  of  the  fourth,  and  presumably  higher,  harmonic  occurs. 
However,  it  should  be  noted  that  they  achieve  a  substantial  improvement  over  the 
basic  propeller  noise  calculations  due  to  Gutin.^  This  can  be  seen  m  Figure  12a. 
Their  report  shows  very  clearly  that  the  higher  harmonics  of  the  loading  have  impor¬ 
tant  contributions  to  the  higher  harmonics  of  the  noise.  Loewy  and  Sutton  '  come 
to  the  some  general  conclusions,  but  their  report  apperrs  to  contain  some  errors.  For 
instonce,  they  show  experimental  data  with  peak  frequencies  which  do  not  occur  ot 
the  blade  passage  frequency  and  its  harmonics.  It  is  difficult  to  conceive  of  an; 
mechanism  by  which  this  could  occur.  Certainly  the  explanation  put  forward  by 
Loewy  and  Sutton  of  random  variations  due  to  bewnwash  would  not  explain  if.  Such 
variations  can  cause  only  broadening  of  the  observed  peak  and  not  a  shift  of  its  center 
frequency.  It  appears  that  at  some  point  in  the  onalysis  between  recording  and  final 
pict,  a  frequency  reference  was  mislaid.  The  most  likely  source  of  error  is  that  the 
frequencies  shewn  os  harmonic  intervals  in  their  original  graph  ore  octually  10-Hz 
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intervals.  Ismg  their  observed  fundamental  frequency  of  14.7  Hz,  then,  gives  har¬ 
monics  in  god  agreement  with  the  experimental  data.  This  corrected  version  of 
their  results  is  shown  in  Figure  13.  The  accuracy  of  prediction  of  the  amplitude  of 
the  harmonics  is  improved  slightly  over  their  original  graph,  but  the  levels  are  still 
i  oo  I  ow  . 

Othf  problems  also  arise  in  the  interpretation  of  Loewy  and  Sutton's  results.  The 
very  substantial  near-field  effects  found  by  Loewy  and  Sutton  were  not  observed  by 
Schlegel  et  j|,  nor  are  they  found  in  the  present  work.  It  seems  likely  that  an  error 
occurred  in  Loewy  and  Sutton's  computations,  which  led  to  underestimation  of  the 
far-field  noise.  This  explains,  to  some  extent,  why  predicted  levels  in  Figure  13 
are  lower  than  measured  data.  A  further  problem  arises  in  the  definition  of  higher 
harmonic  airloads  by  Loewy  and  Sutton.  They  use  data  from  Scheiman^  which  gives 
up  to  24  values  per  revolution.  The  basic  Nyquist  rule  for  frequency  analysis  states 
that  only  he  monies  up  to  one-half  the  number  of  available  points  can  be  calcu¬ 
lated.  22  Thus,  in  this  case,  only  harmonics  up  to  the  twelfth  can  be  calculated. 
Harmonics  u o  to  the  tenth  were  given  by  Scheiman .  However,  Loewy  and  Sutton 
plot  harmonics  out  to  the  twentieth.  This  procedure  must  be  in  error  due  to  folding 
frequency  (or  "aliasing")  effects,  and  explains  the  excessively  high  values  of  load¬ 
ing  harmoni;  s  found  by  Loewy  and  Sutton,  which  undoubtedly  went  a  considerable 
way  toward  offsetting  the  far-field  radiation  underestimate  noted  above.  Further 
applicable  theoretical  work  has  been  accomplished  in  References  18,  23,  and  24, 
bid  numerical  estimates  of  direct  application  to  the  helicopter  problem  are  not 
given.  Sec  ion  3.5  gives  a  brief  comparison  of  these  latter  studies  and  the  present 
work . 
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3.0 


THSOKETlCAl.  ACOUSTICS 


tut  G5HHHAL  bOUATiQ^S  50#  SOtfNU  GE NfRATlQN 


.'  ,.cr  :i  ' 


.  —  .xj*  the  ^  e  n  e ra  ^sc.’Ci’Op.^  t/ti  sound  wj il  bft< dfir i.vcd. ffO^it. f. * rs  t. .  ■  •■■■- 

principles.  Throopkouf  this  section,  Tensor  natetfoa.tei fh-ibe„iufB©&t  <  Qfi: ?9i^jR5?l?S*£i?iS 
utili  be  used.  Tais..ftSK*ti«n  is  perrlqpatwy 

mathematics,  For  example.,  the  -ccnnruUy  equation  may  be  wx-ittea  •; 


*£. 

9 1 


a.ipv.i 


wr.ere 


P  = 
t 

Q  = 


v. 

i 


the  density  *  .  ...  ,  •  '  C.-’-V5..-  •;.'. 

time  •■  -  -. *— ■; -  - — ~±r~ 

the  rate  of  introduction  of  masjvper  unit  volume,  .which  moy  vary 
in  any  desired  way  over  space 


a  three-dimensional  Cartesian  coordinate  r eptesenti ng  x.t,.x^,  *j? 
as  i  rakes  on  the  values  1,  2j  3.  £ ...  3  .  . ,„. 

a  velocity  component  representing  v}  ,v2  ,Vy  as  i  tcAes  on  the 
values  1,2,3.  .  '  -  ‘ 


The  Einstein  summation  convention  used  here  requires  that  when,  the  .tome  ind.e*.«*_ 
repeated  in  any  term,  then  the  summation  of  that  term  over  all  values  of  the  codex  is 
necessary.  Thus,the  second  term  in  the  above  equation  can  be  expanded  as  *  3.  . 


3  (pv.)  3(pv^} 

3x.  3x, 


3(pv2)  3(pv,} 

3x„  T  3x„ 


The  simplicity  achieved  by  using  this  notation  is  apparent. 

Similarly,  the  equation  for  the  conservation  of  momentum  can  be  written 


3(pv.) 


3(pv.  v. )  3  p,. 

_ L_i_  +  - UL_  = 

3  x .  3  x.  i 


J 


J 


(8> 


where 


(i  *  1,2,3)  the  components  of  external  force  per  unit  volume 
acting  over  the  fluid. 
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the  nine  component  stress  tensor  which  includes  the  viscousand 
internal  pressure  forces  on  the  fluid. 


In  Equation  (8)  the  suffix  i  can  again  take  on  the  values  1,  2,  3,  but  here  each 
different  vclue  implies  a  separate  equation.  Equation  (8)  is  a  form  of  the  compressible 
Navier  -Stokes  equations.  Both  the  second  and  the  third  terms  on  the  left-hand  side 
of  Equation  (8)  have  a  repeated  index  j  and  therefore  require  summing.  Writing  out 
Equation  (8)  in  the  long,  conventional  notation  would  require  a  total  of  three  equations, 
each  v,'  th  eight  terms.  Again,  the  advantages  of  the  present  Tensor  notation  are  clear. 


Differentiating  Equation  (7) 
and  subtracting  gives 

with  respect  to 

t  and  Equation  (8)  with  respect  to  x. 

- 

oQ  °Fi 

3i2(pv.v.  +  p..) 

,  1  J  >J 

3t2 

3 1  3x. 

! 

3x.  3x. 

1  J 

Again,  differentiation  of  Equation  (8)  with  respect  to  xj  introduces  a  double  suffix 

i,  which  requires  summing.  Doing  the  same  operation  without  Tensor  notation  would 
have  required  differentiation  of  three  original  equations  by  three  different  variables 
followed  by  summing,  but  in  the  present  notation  the  result  is  obvious.  In  order  to 
derive  rhe  equation  for  sound  generation, the  term 

°o  02  P/dx) 


25 

is  subtracted  from  each  side,  first  done  by  L*ghth  i  1 1 ,  giving  finally  the  general 

equation  for  sound  generation  as 


wlv  re 


02 


SQ 


i 

3t2 

'  °°  ax’  =  8t 

J 

T..  ■= 

p  V.  V. 

+  P-.  -  an  P  6:. 

'J 

1  J 

IJ  !J 

6..  = 

1,  i 

II 

II 

O 

> 

\ 

3F.  d2  T, . 

aT  + 


I 


8x.  9x. 
1  J 


(9) 


the  speed  of  sound  in  the  undisturbed  fluid. 


The  lef:  -hand  side  of  Equation  (9)  is  the  wave  equation  (92/3x2  s  y2),  and  the 

right-bend  side  can  be  regarded  as  c  collection  of  acoustic  source  terms.  The  wave 
equation  is  given  in  terms  of  a  (fluctuating)  density  p  but,  if  desired, can  be  easily 
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converted  to  pressure  p  by  putting  p  -  a2  p,  which  will  apply  in  most  practical 
problems. 


It  m  y  be  observed  that  p  (or  equivalently  p)  occurs  on  both  sides  of  Equation  (9), 
so  ti  at  in  principle  Equation  (9)  cannot  be  solved  directly,  in  practice,  the  terms  on 
the  ight-hand  side  may  be  regarded  as  known,  and  expressions  for  the  sound  field  can 
be  obtained  using  the  well-known  solutions  to  the  Wave  Equation. 


Each 
sour 
hel  u 
radi 
mats 
rote: 
in  fr- 
numl 
T;j , 
pen  a 
term 
helic 
poss’ 
thes 
How 


term  on  the  right  hand  of  Equation  (9)  represents  a  different  possible  acoustic 
e  mechanism.  The  first,  9Q/9t,  gives  the  effect  of  mass  introduction.  In  a 
opier,  an  example  could  be  the  tip  jet  rotor.  Mass  sources  are  the  most  efficient 
tors  of  sound  at  low  speeds.  The  second  term,  9F./9x.  (note  the  necessary  sum- 
3n),  gives  the  effect  of  fluctuating  forces  acting  on  the  air.  Since  the  helicopter 
everts  substantial  fluctuating  forces  on  the  air,  this  is  the  term  of  prime  interest 
a  present  study.  The  third  term,  92  T;^/9xj  9x-,  incorporates  a  rather  large 
;er  of  effects,  the  most  important  of  which,  in  general,  is  that  of  turbulence, 
ay  be  regarded  as  an  acoustic  stress  tensor.  Note  that  since  i  and  j  may  inde- 
ently  take  on  three  values,  Tr  actually  has  nine  components.  Fortunately,  this 
is  of  little  interest  in  the  present  work,  and  it  is  unlikely  to  be  significant  in  the 
:opter  noise  problem  unless  the  rotor  is  driven  by  high-speed  turbulent  jets.  It  is 
ble  to  derive  several  quite  interesting  general  features  associated  with  each  of 
>  noise  radiation  mechanisms,  such  as  velocity  dependence  and  directionality . 
aver,  this  will  not  be  done  here.  A  discussion  is  given  in  Reference  9. 


3.2  SOLUTION  OF  THE  GENERAL  EQUATION 


The  solution  to  the  wave  equation  is  well  known.  If  the  right-hand  side  of  Equation 
(9)  is  written  as  g(y),  the  solution  to  (9)  is 


wht :  e 


P 


y 


P 


/[’] 


dy 


a  fluctuating  density 

the  distance  from  source  to  observer 

the  coordinate  of  the  source  position 


00) 


The  symbol  ~  unaer  y  implies  that  y  is  a  vector  quantity.  This  symbol  is  used  because 
it  requires  a  printer  to  use  heavy  (Clarendon)  type,  as  is  usual  for  vectors.  The 
brackets  around  the  g/r  term  are  of  extreme  importance,  since  they  imply  evaluation 
of  * ceir  contents  at  "retarded"  time  t  =  t  -  r/aQ  .  Because  sound  travels  at  a  finite 
Sp,;,  j  through  the  air,  sound  heard  at  the  same  observer  time  from  different  parts  of  an 


19 


exten.ed  source  must  actually  have  been  emitted  at  different  source  times.  Thus, 
unles.  proper  account  is  taken  of  these  retarded  time  effects,  the  acoustic  calculations 
are  invalid. 


'Although  Equation  (10)  does  give  the  solution  in  principle,  the  actual  source  terms  on 
the  ri  ht-hand  side  of  Equation  (9)  contain  differentials,  and  Equation  (10)  thus  re¬ 
quire;  some  further  manipulation  before  it  can  be  used  in  any  simple  manner.  Several 
methods  of  reducing  Equation  (10)  to  a  more  useful  form  are  available  (for  example, 
Refer-  aces  8  and  18),  but  for  the  presi  nt  problems  it  is  mosr  convenient  to  proceed 
along  ihe  following  lines,  which  follow  in  part  a  method  suggested  by  Lighthill 


We  need  to  solve  the  equation 


W(p)  =  - 


where  W  represents  the  wave  equation,  and  where  the  source  terms  on  the  right-hand 
side  c-  Equation  (9)  have  been  specialized  to  include  just  the  force  terms,  which  are 
the  only  terms  of  interest  in  the  present  helicopter  problem.  Consider  the  equation 


W(p.)  =  -  F. 

i  i 


which  ,  by  Equation  (10),  has  the  solution 


Ti/W  * 


Differentiating  (12)  with  respect  to  x.  gives 


3p.\  -  3F. 


Comparison  of  Equations  (11)  and  (14)  gives  the  solution  to  (11)  as 


..  Hi L 

^  3x. 


_ 3 

i _ 2  3x 


;/b 
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The  above  derivation  shows  how  the  solutions  to  the  wave  equqtibrivfar_,'ib<  fprce  Venmy 
may  costly  be  generated  from  those  for  a -simple  source  g  ,jy!nptyb)r  different  iotipnvi 
Thu  method  will  be  used  in  detail  in  Sec'ion  3.4  to  derive -expressions  tot  ./the  noise 
radiation  from  the  fluctuating  farces  on  the  helicopter  rotor,  mw evei ,_iJi  it. of •'«* U_ J 
Equation  15  will  be  used  directly  to  give  o  convenient  computer  toifliuls  .fpfiSP.Tpu^h; 
tat  ion  <4  .*i£_»'iO;;o  from  art  a*t>?tr<M»*y  moving  .-Uit.uw.7.,*  f~"c  T_nn\~- 


_ 3.3  _  THE  MOVING  FLUCTUATING  FORCE 


If  we  now  assume  that  the  forces  are  in  motion,  it  is  convenient  to  specify^the  forces;- 
in  a  moving  frame  of  reference,  tor  example,  on  the  helicopter  rotor.  Suppose  that 
coordinates  measured  in  this  frame  are  defined  by  rt  and  that  the  origin  pf  thtts^  po- X 
ordinates  is  moving  with  velocity  aQM.  Then,  at  ary  instant,  the  q  and  y  coordinate 
systems  ore  connected  via  ~  y  ‘ 


=  y 


-Mat 
~  0 


However,  in  Equation  15  we  are  not  required  to  evaluate  the  integral  at  any  ■  istqnt  i 
Instead,  the  integral  must  be  evaluated  at  the  appropriate  retarded  time  t  -t  -  r/c0; 
i .e. ,  aver 

n  -  y  +  -Mr  -  Mact 


Thus,  in  the  coordinate  transformation  from  fixed  to  moving  axes,  it  is  appropriate.  tP 
use,  as  first  suggested  by  Lighthiii, 


n 


-  y 


+  Mr 


(16) 


This  axis  transformation  also  affects  the  volume  element  of  the  integration  in  (15),  and 
the  integral  must  be  dividod  by  the  Jacobian  of  the  transformation 


J  = 

3n,/9y, 

Srij/dy, 

3n3/^y, 

t 

3n  y/bi 

b1/b2 

3p,/»Y3 

3n2/3y3 

0V»y3 

M;  (x.  ~  y.)/> 

1  |  1  ' 

=  {Ml(xt  - 

y,}+  h^(xj 

=  1 


M. 


where  M 

r 

is  the  component  of  the  convection  Mach  number  in  rhe  direction  of  the  observer. 
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Equation  (15)  thus  becomes 


jL/-r^L_i 

dx.J  I  r(l  -  Mf) 


The  derivative  3/3 x.  is  operating  on  both  an  integral  over  q  and  a  retarded  time 

I  ~ 

operator,  be  hr  of  which  are  functions  of  x.  For  any  function  f(t). 


[f(t)]  =  [f(t  -  r /a0)] 


Thus,  the  partial  derivative  with  respect  to  x.  ,  keeping  q  constant  of  f,  is  given  by 

I 

the  chain  ruie  as 


I*  .  _j _ (  t 

0xi(n)  ®xi(n>  I*1  j* 


h _  J_  3f_ 

C.,  \  °n 

l(q)  0 


X.  -  y. 


(x.  -y.) 

J  J 


From  Equation  (16), 


^ - (x.  -y. )  =  6..  +  M.  7T^- 

«x.,o  J  J  '!  J  oxr/ 


Hence, 


x.  -  y. 

i  '  i 

r(l  -  M  ) 
r 


and  when  this  is  in  Equation  (18), 


— H  ■  & 


X.  -  y. 


1  3f 


,  l  J  3x.,  x  r(l  -  M  )  a  3 1 

i(q)  i(q)  r  0 


Now  the  oroer  of  magnitude  of  the  first  term  on  the  right-hand  side  of  Equation  (19)  is 
1/r  where  r  is  a  typical  distance  from  the  source,  whereas  the  order  of  magnitude  of 
the  second  term  is  f/ac  =  1/X  where  f  is  a  typical  frequency  and  X  is  a  typical 

wavelength  Thus,  when  the  observation  point  is  many  wavelengths  from  the  source,the 
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first  near-field  term  moy  be  ipnorec  .  Hence,  when  (19)  and  .(17)  ore  used,  the  expr**-"- 
4?©n  (or  the  sound  radiation  its  the  far  field  due  to  an  g^t.jroi  ij.^mcv.iog’  fryce  distribution' 

it  ,rv 


r 
-  * 


.1 


r  *  1  2  /  j  a.  r;?  .Mi  at  in.;  -M 

4i8,»  1  «  r  — «- '  r  /  I  •_  --7;— — 

If  desired,  the  time  differentiation  may  be  performed  to  give 


where 


] 

_i_  f 

x.  -  y.  {  3F. 

s  *  t  1 

f.  2M  * 

»  f 

l 

A  rv 

W  J  ^ 

0  fit  -M  I9* 

I  -  Mf  3t 

■i 

vf|  . 

u  f 

•  .(21) 

3M 

r 

x.  —  y.  3  to. 

1  »  < 

■  (22) 

3  r 

r  dr 

is  the  component  of  the  acceleration  (divided  by  q^)  in  the  direction  of  the  observer  . 
Equation  (21}  was  derived  for  the  special  cCse  of  a  point  source.®  Note  that  both 
Equations  (20)  and  (21)  apply  only  if  the  force  distribution  is  not  changing  in  size. 
This  will  be  the  case  for  helicopter  applications. 


Further  discussions  on  the  utilization  of  Equation  (21)  for  computer  calculation  of 
sound  levels  is  given  in  Section  5.0  .  In  the  next  section,  analyticsolutiomfor  the 
noise  radiation  From  a  helicopter  will  be  given  using  methodswhich  ere  an  extension 
of  those  presented  In  Section  3.2. 

3.4  RADIATION  FROM  A  HELICOPTER  ROTOR 


As  shown  in  Section  3  .2,  the  sound  radiation  from  a  fluctuating  force  can  be  found 
from  the  solution  for  a  simple  source  by  differentiation.  In  the  helicopter,  fluctuating 
axial,  circumferential,  androdiai  components  of  force  occur.  Rather  then  calculate 
the  effects  of  each  force  component  separately,  it  is  convenient  to  .ofve  the  problem  for 
the  rotating  simple  source  first  and  then  just  to  differentiate  in  ■'he  appropriate  direction 
to  determine  results  for  the  various  force  components. 

Suppose  that  the  fluctuating  source  is  defined  by  a  Fourier  series 
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X  =  l 


23 


where1  8^5* the -angle  oyoyncHhe' j  otoc'diic  from  the  y*cxis  .and  is  eq«pi  iv  Q  « 

the  anauivJf  veioci.i-y  of  the' rotor.  It  is-  convenient  fcr  the  'present  premiere 
t©  define  oi*o  <i  cpfTipi.dx.r^yj^rjc^i^  by  - 
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-  Since  the  source  j$  denned  in  terns  of  /racing  source  coordinates,  it  is  necessary  to 

-  .transform.  Equation  (i 5)  from  y  to  n  coordinates  using  Equation  (16) .  in  q  manner 
'  idfoticatto  that  'used  to  obtain  Equation  (17),  we  find  that 
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_  -wberetlse  relation  p  -  ’a^p  has  also  been  used. 

Specializing  to  a  point  source  causes  the  n  integration  to  vanish,  for  any  harmonic 
.funytion'the  Fowler  coefficients  are  given,  in  complex  form,  by 


C  =  a  +  i  b  =  |  f  exp  i  n  u  t  d  t 

rt  n  n  n  u 


where  the  integral  is  over  any  period.  Replacing  f  by  the  point  source  version  of 
Equation  (26)  give*  the  complex  magnitude  of  the  sound  harmonic  os 
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n  Q  t  dt 


Now  change  from  observer's  time  f  tc  source  time  t  by  t  - t  -  r  /d,  , 
dt  =  dt{l  -  M  )  ,  giving 


(27) 
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»herc  n  t  hgs  beer  replaced  by  9. 
;  con  oe  seen  from  Figure  14  ;ha? 


r2  -  xJ  r  Y'  ^  -  2  Ys  cos  ;,6  -*)  •  -  •  (2?} 

where  x,y  =  Cartesian  coordinates  *ith  x  along  the  rotor  ox.li 

Y  -  the  radial  distance  of  the  observer  from  the  rotor  axis. 

R  —  the  effective  radius  of  the  point  source 

41  -  a  reference  angle  between  the  y-axts  Qnd  the  observer 

If  the  observer  is  far  from  the  rotor,  so  that  ix2  A  Y2  )  >>  R2  t  then  Equation  (2?)  may 
be  approximated  to  give 
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t  =  f!  -  — —  cos  (6  -  9) 


(30) 


*%  c  .  ... 

where  r{  -  (x2  +  Y2)  is  the  distance  of  the  observer  from  the  rotor  hub. 

Using  fhe  approximation  of  Equation  (30)  and  the  source  tenms  defined  by  Equation 
(24)  in  Equation  (28)  gives 
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which  can  be  written 
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Since  :he  integral  applies  ove  *  erva!  2-n,  it  can  be  expressed  in  Bessel  function 
form,  using  formula  42  in  M- 


Equation  (33)  gives,  in  complex  form,  the  sound  harmonics  from  a  point  simple  source 
describing  a  circular  path.  The  equation  could  be  applied  directly  to  the  calculation 
of  rotating  mass  sources,  provided  that  the  time  differential  is  observed  (see  Equation 
(9)  anc  that  proper  account  is  taken  of  momentum  output  (see  References  8  and  9). 
Howev  ;r,  we  are  interested  in  deriving  the  results  for  the  force  cases,  which  are 
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Only  the  far-field  terms  have  been  retained  in  Equations  (34). 

The  prime  on  the  Bessel  function  in  the  radial  expression  denotes  differentiation; 

M  =  Q  R/oq  is  the  rotational  Mach  number  of  the  point  force.  Negative  signs 
must  bo  applied  in  the  first  two  equations  because  differentials  are  based  on  observer 
coordinates,  whereas  the  differential  in  the  last  equation  is  on  a  source  coordinate. 

Notation  in  Equations  (34)  must  now  be  changed  to  specify  the  forces  acting.  The 
three  components  by  simple  Fourier  series  are  defined  r  s 
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The  thrust,  drag,  and  outward  components  of  force  ore  assumed  to  act  in  the  axial, 
circumferential,  ana  racial  directions  respectively.  In,  conversion  to  thrust,  amines 
sign  rrusr  be  ir.corporoted  because  fne  farce  on  the  ah  is  in  the  negative  x-aireciion. 

Equations  {34}  are  converteo  to  the  required  form  using  Equations  (25).  Note  that 
terms  for  both  pius  and  minus  X  in  the  summations  in  (34)  contribute  to  th«  resuit 
for  any  given  loading  harmonic  ,  Thus,  the  final  resuit  tor  trie  complex  magnitude 
of  the  sound  harmonic  is 
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where  the  ocditlenal  subscript  X  implies  that  only  the  one  X  ;oad*rg  hormpnit  »&•_ 
considered .  The  complete  sound  field  :s  given  by  the  summation  cf  sound  radiated  by 
a!i  harmonics.  •-  ■._■  ? 

3.5  COMMENTS  ON  THE  RESULTS  . 


Equations  (37}  give  the  complete  solution  tor  the  sound  radiated  by  fluctuating  air-  c 
loads  on  a  rotor.  The  results  are  limited  .n  two  ways.  First,  the  equations  do  not 
cpply  close  to  the  rotor,  soy,  within  a  few  wavelengths  or  rater  diameters.  This  <s 
not  important  for  most  practical  noise  problems.  Second,  they  do  not  Include  the 
effects  cf  biade  motion.  This  point  has  beer  studied  using  the  general  computer 
program  based  on  Equation  (21)  and  is  also  studied  analytically  In  Appendix  I.  'One  ' 
extremely  important  effect  not  explicitly  giver,  in  Equations  (37)  is  the  effect  of 
blade  number.  If  B  blades  ore  present,  harmonics  which  are  not  integral  multiples 
of  B  will  cancel  out.  Those  harmonics  which  are  multiples  of  B  will  odd.  Thus, 
the  effect  cf  blade  number  rrcy  be  included  in  Equation,  (36)  ond  (37)  by  replacing 
r.  by  rr  B .  In  this”  case,  the  coefficients  of  the  force  harmonics  must  be  token  as 
the  values  for  the  complete  rotor,  which  ore  8  times  the  values  for  the  individual 
blades . 


It  is  of  particular  interest  to  consider  the  case  X  =  0  in  Equation  (36).  This 
corresponds  to  the  case  of  steady  loading  only,  as  is  assumed  for  a  propeller. 
Using  X  =  0  in  Equation  (36)  gives  substantial  cancellation,  with  the  final  result 
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The  first  term  in  the  above  equation  is  identical  with  the  classical  propeller  noise 
solution  due  to  Gutin,  while  the  second,  the  radial  component  term,  is  the  same 
as  that  deriveJ  In  Reference  9.  The  reduction  of  the  general  solution  to  particular 
cases  previously  obtained  is  a  helpful  test  of  the  mathematics. 

The  effect  of  forward  velocity  is  also  of  considerable  interest.  In  Reference  8  it  was 
pointed  out  how  the  equations  for  constant  velocity  convection  of  the  hub  could  be 
obtained  from  those  for  the  stationary  case  by  replacing  the  term  r?  in  the  stationary 
case  by  r(  f  ■  -  \^r),  where  MQr  is  the  component  of  the  hub  convection  Mach 
number  in  the  direction  of  the.  observer.  In  utilizing  this  transformation,  it  is  import¬ 
ant  to  note  that  it  applies  to  the  retarded  position  of  the  helicopter.  In  other  words, 
the  dimension  rj  used  must  be  taken  as  the  distance  from  the  observer  to  the  position 
of  the  helicopter  when  it  emitted  the  sound.  Relation  of  the  results  to  the  instantane¬ 
ous  position  of  the  helicopter  requires  another  transformation.  Details  are  given  in 
Section  5.1  c  id  in  the  companion  report  .^0  It  may  also  be  noted  that  virtually 
identical  transformations  were  presented  in  Reference  18,  where  it  was  also 
shown  that  the  above  (1  -  Mor)  correction  term  gave  the  Garrick  and  Watkins,^ 
moving  propeller  result  directly  from  that  of  Gutin, ^  for  the  stationary  case. 

Perhaps  it  should  also  be  noted  that  the  results  of  Equation  (36)  are  not  entirely  new . 
During  die  course  of  the  current  study,  a  major  new  book  on  acoustics  by  Morse  and 
Ingard,  was  published  .  Equation  1 1 .3.20  of  that  book  gives  the  sound  radiation 
by  a.  propeller  in  unsteady  flow,  which  corresponds  to  the  helicopter  rotor  case. 

Only  thrust  and  drag  terms  are  considered  in  that  equation.  The  results  are  derived 
in  a  rather  different  manner,  but  they  agree  with  those  derived  here.  No  numerical 
analysis  of  the  results  is  given,  but  several  additional  points  of  interest  are  discussed. 
A  report  by  Arnold  etal,  ^  also  meats  a  similar  case,  the  problem  of  the  "singing" 
propeller  in  underwater  acoustics.  Here,  frequencies  other  than  the  fundamental  were 
allowed  to  occur,  and  it  was  found  that  multiple  frequencies  were  produced  by  fre¬ 
quency  modulation  effects.  Finally,  a  recent  report  by  Lowson,  covers  the  related 
case  of  jet  engine  compressor  noise.  Again,  only  thrust  and  drag  terms  were  con¬ 
sidered,  but  results  were  obtained  by  a  method  different  from  that  given  here  and 
are  in  agreement.  Analytic  expressions  for  overall  acoustic  power  radiation  were 
also  given  which  can  be  of  use  to  the  present  problem. 
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ROTOR  AERODYNAMICS 


If  riSi  noted  <n  Section  2.  I  that  the  feosic  source  of  helicopter  ooiie  is  i!i*  fT.oj-T  Ijfhnj 
rotor  (or  rotors),  end  a  brief  discussion  of  the  general  peobiern  «»  p?ssc“t©dt-  ".'ppoirt 
this  statement.  In  this  section  a  closer  took  is  taker,  at  the  fundamental  oe?adynomi_cjIr 
phenomena  whloh  generate  rotor  noise.  cr.d  suiiubic  techniques  for  evniijatinpdhetr 
magnitudes  ere  examined  .  It  is  not  intended  to  present  a  rigorous  treatment^ptblhs  V. 
subject  sf-rotor..oerpdynwnics?_  this  is  a  complex  subject  in  its  own  right,  but  the  fun¬ 
damentals  ore  covered  ir,  sufficient  detaii  rb  serve  os  an  sniraduvilor.  to  rater  noise 
generation  end  to  postulate  some  realistic  source  terms  for  the  noise  theory.  Emphasis 
is  placed  on  the  acoust'c  implications  of  each  aspect  of  the  problem. 


The  notation  in  this  section  is  somewhat  different  from  that  in  preceding  sections, 
since  an  attempt  has  been  mode  to  use  conventional  helicopter  aerodynamics  symbols 
wherever  possible  (in  particular,  following  References  27  and  28).  Symbols  used  in 
this  section  that  have  different  meonings  elsewhere  are  listed  belcw  . 
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lift  curve  slope 

coning  angle,  rods 

blade  chord,  ft 

effective  root  chord,  ft 

effective  disc  area  correction  «  s2  -  s2 

2  1 

28 

blade  twist  integral  =  f. 

harmonic  number 

blade  radiai  position,  ft 

nondimensionat  radiai  station  -  r/R 

station  defining  inboard  end  of  blade 

station  aefining  effective  outboard  end  of  blade 

blade  taper  integral  =  4  J  2  —  5'  1  dE 

s  c 

disc  area  =  nR*  ,  ft2 

longitudinal  cyclic  pitch  coefficient,  rads 
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;  Jcijrql  cyclic' fitch  'coefficient,  rods 
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bled*  pitch  ,  rads 
cojUct'v^. pitch  (nominal),  rads 
poitective  pitch  jne<u«r*d  at  blade  root,  rods 
total  biod*  twist  roof  to  tip  (washout) 

■»  -  --v, 

dawnwosh  coefficient 
advene*  ratio  =  V/QR 
Wade  angle  of  attack,  rods 
rotor  azimuth  angle,  rods 
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4  • 1  NOISE  GENERATION  BY  AERODYNAMIC  FORCES 


We  have  seen  in  Section  2.1  that  the  major  contributions  to  rotor  noise  are  made  by  the 
fluctuating  lift  and  drag  forces  which  act  on  the  rotor  blades  and  are  constantly  in 
motic  > .  The  basic  equation  of  sound  generation  by  an  aerodynamic  force  in  motion  is 
(Equc.  ion  20) 
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when  the  space  integral  has  been  removed,  since  the  F.  represent  a  point  force.  Thus, 

the  sound  pressure  (or  density)  fluctuation  at  any  point  due  to  the  force  is  a  function 
of  its  position  relative  to  the  observer  (x|  -  y. ) /r,  its  rate  of  change  (3Fj/9t),  and  its 

velocity  and  acceleration  toward  the  observer  (Mf,  9Mr/3t).  The  F;  are  of  course 
the  three  orthogonal  components  of  the  resultant  force  acting  at  the  point  y. .  In 

reality,  *he  forces  acting  on  a  rotor  are  distributed  pressures,  and  it  is  only  for  mathe¬ 
matical  convenience  that  we  perform  an  integration  of  these  pressures  to  arrive  at 
equivalent  point  forces.  If  is  also  for  convenience  that  the  net  aerodynamic  load  on  a 
rotor  blade  section  is  conventionally  broken  down  into  the  two  components  known  as 
lift  end  drag.  It  will  be  seen  that  these  components  are  not  particularly  useful 
for  the  numerical  acoustic  calculations,  and  further  resolution,  into  a  coordinate  system 
which  is  fixed  in  space,  is  preferable. 


However,  Equation  (20)  shows  the  informa¬ 
tion  required  for  a  sound  calculation  based  on 
a  single  source.  A  helicopter  rotor  represents 
a  complex  distribution  of  aerodynamic  pres¬ 
sures  in  spaceyand  it  is  necessary  to  make  use 
oF  a  distribution  of  sources  which  lie  on  the 
rotor  blade  axes.  The  sketch  shows  how  the 
distributed  lift  and  drag  loads  acting  on  the 
blades  are  simulated  by  a  series  of  point 
forces.  V/hat  this  diagram  does  not  show  are 
the  variations  and  motions  of  these  forces  as 
the  Made  rotates. 
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The  hovering  rotor  is  relatively  straightforward  aerodynamical ly  since  the  forces  acting 
on  the  blades  are  effectively  steady.  Provided  the  rotor  is  operating  in  undisturbed 
free  space  ,  the  flow  is  completely  axisymmetrical ,  and  there  are,  in  principle,  no 
perturbations  to  change  the  steady  forces  as  the  blade  rotates.  In  forward  flight,  how¬ 
ever.  when  the  rotor  is  essentially  sideslipping,  providing  both  lift  and  propulsive 
force  for  fde  helicopter,  rhe  flow  over  the  blades  is  asymmetric,  due  to  a  velocity 
differential  over  the  advancing  and  retreating  blades.  Rotor  control  is  obtained  by 
"cyclic  pi  ch"  change  -  which  is  the  name  given  to  the  first  harmonic  variation  applied 
to  the  blare  pitch  angle  as  rotates.  Since  the  relative  air  velocity  over  the  blade 
also  has  a  first  harmonic  variation  and  since  aerodynamic  forces  are  proportional  to 
the  square  of  the  relative  velocity,  we  may  expec*  to  find  at  least  three  harmonics  in 
the  force  fluctuations  acting  on  the  blades.  This  would  be  true  if  the  flow  through 
the  rotor  v  ere  uniform.  However,  due  to  the  proximity  of  the  rotor  to  its  own  vortex 
wake,  which  is  swept  backwards  under  the  rotor  disc,  the  flow  is  far  from  uniform, 
and  velocity  fluctuations  are  induced  which  give  rise  to  very  many  harmonics  of  blade 
tbading.  The  calculation  of  these  higher  harmonic  blade  loads  is  an  extremely  cc'i- 
plex  problem  which,  to  date,  has  been  performed  numerically  only  by  digital  com¬ 
puters  (foi  example.  References  29,  30,  and  31)  and  then  with  only  limited  success. 
Such  calculations  were  beyond  the  scope  of  this  study  and,  as  will  be  shown,  would 
not  be  justified  in  any  case. 

As  a  first  ^tep,  a  simplified  analysis  is  followed  which  helps  to  clarify  the  general 
problem . 

4 . 2  ESTIMATION  OF  BLADE  LOADS  AND  MOTIONS  AS  A  FUNCTION  OF 

FLIGHT  CONFIGURATIONS 

One  of  the  main  prerequisites  to  performing  the  noise  calculation  is  to  compute  the 
lifting  rotor  attitude  in  space,  together  with  the  coning  angle  which  dominates  the 
blade  flapping  motion  terms.  A  procedure  by  which  this  information  can  be  accurately 
derived  follows  ■ 


Force  Balance 
The  start!  ig  point 


in  these  calculations  is  the  flight  configuration  defined  by: 
W  ’  -  aircraft  gross  weight,  lb 
V  =  forward  velocity,  ft/sec 
Cq  -  fuselage  parasite  drag  coefficient 
80  -  blade  section  profile  drag  coefficient. 
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and  the  disc  incidence  is 

.  _  .  _i  /H 

•d  ~ tan  r 


When  the  rotor  thrust  vector  is 
assumed  to  be  normal  to  the  tip  path 
plane  (and  experiment  has  shown  this 
to  be  true  to  within  1°  according  to 
Reference  27)  and  when  the  rotor 
weight  with  respect  to  the  total  is 
ignored,  the  shaft  inclination  is 


where  D  =  fuselage  drag  and  H  =  rotor  drag. 


Defining  a  rotor  thrust  coefficient  Cj  =  T/0.5  p  A  QR2,  a  fuselage  parasite  drag 
coefficient  -  H/0.5  pAV2,  a  rotor  drag  coefficient  =  H/0.5  pAV2, 

Equation  (39)  can  be  rewritten 


J 


Payne^  (p.  202)  gives 


ch0=  + 


1 .7  /  t,  S 


l  wo  3  UJ  6 
+  — - — 


(40) 


wh, 


6  =  6. 


1  + 


2  r  7 

°  CTb 


|f32  +  U2  (7  t;2  +  t,  tj)} 


Induced  Flow 


Wh  -n  the  disc  angle  of  attack  is  known,  actuator  disc  theory  is  used  to  calculate 
the  induced  flow  through  the  rotor.  The  momentum  relationship  is 
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T  -  2pA  v  V2  +  v  2 
r  0  f  0 

A  refinement  is  obtained  by  introducing  the  concept  of  effective  area  which  uses  aero¬ 
dynamic  blade*  root  and  tip  stations  to  account  for  tip  losses  and  the  blade's  ending 
distance  from  the  rotor  hub,  e  =  s2  -  s2  .  Using  this,  the  above  equation  can  be 
written 


X«0  [Xio  +m2|S  =CT/4e 


A  more  realistic  induced  flow  distribution  is  one  that  is  not  uniform  but  increases 
linearly  from  the  front  to  the  rear  of  the  rotor  according  to  Glauert's  classical 
equation 

Xj  =  Xjo  (1  +  K  x  cos  i{i)  (42) 
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where  K  is  given  by, 

(4/3)  (MA) 

|.  K  =  1.2+(MA) 

and  X  is  the  feral  mean  inflow,  p  sin  i^  +  Xj0  . 


Rotor  Control  Settings 

Conventional  blade  element  theory  is  used  to  establish  blade  section  loadings  (see 
Figure  15), 

#  =  -  &  -  Aj  cos  t|»  —  B|  sin  iji 

where  is  fh  *  collective  pitch  setting, and  Aj  and  B|  are  cyclic  pitch  coefficients. 
The  normal  an  tangential  velocity  components  ore,  putting  *  =  r/R; 

'  Up  =  X  0  R 

Uj  =  Q  R  (s  +  u  sin^) 

The  lift  increment  dL  on  the  blade  element  of  width  d  r  is 
dL  =~p  (U2  +  U*)  CL  c  dr 

-ypUjOacdr  (since  Up  <<  U-p)  (43) 
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(since  Up  <<  Uj) 


wrefe 


yr  -  9  &  &  - 


Substitution  for  Uy  end  a  in  this  equation  eventually  yields 

r 

—  -  -r  p  i!R‘  ac  R  M  C  s“  p‘  }  vt>n-  #Ts)  -ux6,-  x  >w  \ 

c;  ,  L1  "  ’  ' 

"  |  ( ‘2  +  \  H2)  A,  *  X-c  K  *2  I  cos  « 

T  |2  F  s  "  Is2  '  4~h2)  Bi  '  1"  Xo|  s'n  ♦ 

n  jp  s  Bj  -  ~  p2  |  cos  2$  -  |h  s  A,  -ri^Ks  \c  |  sin  2* 

+  |j  H2  A,  |  cos  2  if  *  j~  p2  S,  j  sin  3*j 


This  analysis  assumed  zero  coning  angle .  If  the  coning  angle  is  a5  ,  then  the  inflow 
X  is  increased  by  a  factor  p  ac  cos  $  .  Thus, the  coefficient  X;5  K  s  is  increased  to 

(X;o  Ks  tfiac).  This  introduces  the  additional  terms  t  p  s  in  the  cos  p  bracket 

and  +  1,  2  p2  ac  in  the  sin  2$  bracket. 


Tne  blades  aaoot  a  "coning"  angle  to  balance  lift  end  centrifugal  forces  (in  the  case 
of  the  conventio;  si  hinged  rotor),;  and  when  moments  about  the  flopping  hinge  are 
equated,  it  is  round  by  the  spanwise  Integration  of  the  steady  term  of  Equation  (44)  that 

-  k4  m  (t*  -r  j  p2  t2)  -  “  a2  kj  -  (p&,  +  X)  t3 

a„  = - ; - - - f -  (A5) 

r 

Since  this  analysis  is  conducted  with  respect  to  the  rotor  disc  plcne,  i.s.,  the  "plane 
of  no  flopping",  the  first  harmonic  section  loadings  must  integrate,  over  the  length  of 
the  btede,  to  zero.  This  yields  two  further  relationships 


Aj  -  -  ' 


(pt3o  t).  <t4) 

(>  *4'2  *2) 


(AS) 


and 


&;  - 


**3  h*  1st  ’  \^Tb  *  *3  N  T  *-h) 

H'i 


rirwfty,*  th'i  i'^ady  farm  of  Equation  (44;  a  Ini  eg;  tiled  ?c-  yield  the  tohs?  rc.«w  thrusi, 
which,  ofie;  JvbititufTvio  -if  the  blece  thrust  coefficient  Cy,  jives  the  following  exptej- 
stop  tor  the  coitecf;v«  c:rcrr  —  -  -  - 


V 


L  , 

\*«e  '  i  h 

*2.!  (Cfb  •  *3  ■ 

1 

'  *  -* 

/ 

v 

■  )  **  ux  »2  i  2  k-j  *  A  tj.! 

! 

f  1  2  \ 

i  ! 

ivfh2  v 

1  *  2  p*  r2  t. 

(48) 


Equations  (48),  (47),  (45t,  and  (46)  ere  solved  sequentially,  in  that  order,  !or 
&t  /  %  and  A,  , 

A*  noted  previously,  these  quantities  ore  derived  with  respect  to  disc  axes  •  ;?  the 
cyclic  pitch  coefficients  ere  required  with  respect  to  ths  shaft  exes  (which  is  mere 
physically  roearvi:igfut)rthe  following  transformations  -rust  be  applied. 


A>  -Ai  tb,  ana  8-  =  S,  -  a, 

^5  ‘  'S  1  S  1  *  5 


where  dt  s  and  Qj ,  ore  flapping  Coefficients  measured  relative  to  the  shaft  axes  ,  in 
this  cate* 


‘d 


cnc 


4.3  LIMITATIONS  Or  SIMPLE  THEORY 


The  simplified  analysis  of  the  preceding  section  has  been  used  to  ccicuiaie  the  con¬ 
trol  settings  ana  blade  loads  for  the  H-34  helicopter  in  a  variety  of  flight  conditions . 
The  computed  results  have  been  compared  with  -he  experimental  date  horn  Reference 
21,  and  some  of  the  comparisons  ore  shown  in  F'gures  16  trough  16.  Figures  16(a) 
through  16(d)  are  the  variations  of  disc  ancle  of  cttack  i^,  collective  pitch  Ct  ana 


cyclic  pitch  coefficients  A  and  with  advance  rcrio  p. 


The  rr.cst  noticeable 

$  's 

feature  of  these  plots  is  the  apporuntiy  large  scatter  o’  the  experimental  data.  How¬ 
ever,  it  should  ue  noted  thot  these  points  correspond  to  a  large  number  of  flights  of 
the  tame  helicopter  under  different  arrr.ospheric  and  loading  conditions  ana  slightly 
different  rotor  speeds.  No  attempt  has  been  mace  to  normalize  these  results, 
although  such  a  step  would  probably  reduce  the  scatter  somewhat  .  The  main  rea>or 
for  nor  doing  so  is  that  insufficient  information  ■  $  available  :n  Keferen£.o  xl  rc  per¬ 
form  un  adequate  normalization  .  However,  >n  spite  of  th.s  »cct*ef ,  :t  con  be  seen 
that  the  agreement  between  rneory  one  experi  rent  ^  --air  . 


Figures  17  and  18  show  the  sine  and  cosine  components  of  the  first  ten  blade  section 
airland  harmonics  for  the  two  advance  ratios  p  =  0.15  and  p  =  0.30.  The  experimental 
poin  s  were  obtained  by  the  integration  of  a  chordwise  distribution  of  differential  pres¬ 
sure  measurements  and  therefore  represent  the  normal  blade  section  loading .  The  theo¬ 
retical  curves  represent  section  lift  loadings  which  are  not  exactly  aligned  with  the 
normal  force,  although  of  course  the  numerical  differences  in  these  two  forces  are  very 
smaP  . 

In  both  velocity  cases,  the  agreement  between  theory  and  experiment  is  good  for  the 
steady  loading,  fair  for  the  first  harmonic  components,  and  poor  for  the  second  har¬ 
monics.  Little  can  be  said  about  the  third  harmonic  comparisons,  since  there  is  very 
little  correlation.  The  simple  theory  of  course  predicts  no  harmonics  greater  than  the 
third.  This  is  quite  clearly  unrealistic;  since theamplitudes  of  the  harmonics  up  to  the 
tenth  are  seen  to  be  significant  experimentally.  The  main  conclusion  io  be  drawn  from 
this  comparison  is  that  the  simple  theory  based  on  momentum  concepts  is  adequate  for 
the  calculation  of  rotor  performance  parameters  such  as  thrust,  disc  angle  of  attack, 
coning  angle,  and  control  settings,  but  it  is  quite  inadequate  for  the  prediction  of 
the  harmonic  airloads  which  are  so  vital  to  the  noise  problem. 

» 

The  deficiencies  in  the  method  are  fairly  obvious .  It  is  important  to  remember  that  the 
momentum  theory  is  simply  a  crude  analogy  for  the  induced  flow  through  fhe  rotor  and 
that  the  main  justification  for  it  is  that  it  gives  good  agreement  with  experiment  in  per¬ 
formance  calculations.  To  obtain  a  more  realistic  picture  of  the  induced  flow,  wemust 
turn  ro  a  consideration  of  the  rotor  wake  and  of  the  vortex  sheet  theories.  Again,  in  a 
somewhat  crude  fashion,  it  may  be  considered  that  if  the  circulation  about  a  certain 
length  of  blade  changes  instantaneously,  then  an  elemental  vortex,  of  equal  but 
opposite  strength  to  the  increment  in  the  bound  vortex,  will  be  shed  from  the  biade 
with  its  axis  parallel  to  that  of  die  blade.  Further,  if  there  is  a  difference  in  the 
circulation  about  two  adjacent  segments  of  the  blade,  then  a  vortex  will  be  trailed 
from  that  intersection,  with  its  axis  parallel  to  the  flow,  and  with  a  strength  equal  to 
the  difference  between  the  two  bound  vortex  strengths.  Now,  in  practice,  the  liftvaries 
continuously  along  the  blade  radius  as  the  blade  rotates  with  a  corresponding  variation 
in  circulation.  Thus,  the  wake  behind  each  blade  may  be  thought  of  as  a  lattice  struc¬ 
ture  of  shed  and  trailing  vortex  elements.  In  fact,  such  a  model  has  been  included  in  a 
digital  computer  program  for  the  numerical  calculation  of  vortex  wake  effects,^  and  the 
sketch  on  the  following  page  illustrates  the  concept. 

If  this  wake  structure  could  be  accurately  and  realistically  defined,  then  the  induced 
flow  in  the  vicinity  of  each  rotor  blade  could  be  calculated  by  summing  the  velocity 
components  induced  by  each  vortex  eleme.  Although  this  method  represents  an  order 
of  magnitude  of  improvement  over  earlier  attempts  to  account  for  the  woke  flow  (which 
were  necessarily  oversimplified  in  the  absence  of  powerful  computing  equipment),  it  has 
been  frustrated  by  further  factors  which  are  difficult  to  include.  One  of  the  funda¬ 
mental  difficulties  is  that  changes  in  the  bound  vortex  strength  do  not,  in  practice,  occur 


'  ,  ;in  steps-  Alt  changes  are  continuous  and  the  real  wake  behind  the  blade  is  a 

vortex  sheet .  a  continuous  surfoce  of  vorticity.  However,  this  sheet  is  unstable  and 
.welj^nown,  tolls  up. into  at  least  one  vortex  tube,  and  possibly  more,  bach 
..y  “tube  .rqpves -away  from  the  rotor  in  the  form  of  o  helix,  which,  because  of  inter- 
'•  action* feefwteejl  different  vortices,  rapidly  becomes  distorted.  Although  in  principle 
'"Titjs  postibleJo  calculate  these  distortion^  o>--rs!!  problem  becomes  immensely 
•”  2  corr^lex .  Furthermore,  the  hosic  heticai  vortex  shape,  even  in  hover,  is  unstable  for 

-  eitchet  below  0.34 Reference  33). 

-  t  v»  — 

-  :"-ilTl»e  above  discuision  iliysnotesi  just  one  of  the  problems.  Further  difficulties  include 
compressibility  effects  near  the  advancing  blade  tips,  reverse  flow  over  the  inboard 
■-•••  end  ofth* -retreating  blade,  blade  stall  roar  the  retreating  blade  rips,  and  inter¬ 
ference  effects  between  the  rotor  and  the  fuselage  and  between  rotors. 

-A  factor,  which  hot  a  profoundaffect  upon  the  wake  colculaf  ions  descr  Ibed  above ,  and  wh  tch 
also  has  a  direct  influence  on  the  rod ic ted  noise,  is  the  dynamic  response  of  the  rotor 
biodata  Blade  mouritin  '  mechanisms  take  a  variety  cf  forms,  the  most  complex  of  which 
i»  the  flap- log  binge  system .  In  this  system,  the  blade  has  freedom  of  rotation  about  its 
root  both  in  and  out  of  the  rotor  plane.  Further,  the  blade  is  flexible,  and  thus  bends  in 
normal  and  transverse  directions  and  twists  about  its  pitching  axis  .  The  fluctuating 
airload  environment  causes  oscillatory  motions  in  all  these  degrees  of  freedom.  Aero- 
■  dynamically,  these  motions  modify  the  applied  airloads  and  subsequent  vortex  shedding 
.action*/ so  thpt  the  entire  blade  loading/response/wake  system  forms  a  closed-loop  prob¬ 
lem-  Acoustically,  the  blade  response  generates  sound  directly  through  the  motions  of 
the  aerodynamic  noise  sources,  in  addition  to  the  indirect  consequences  of  the  modified 
airloads . 
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4-il  AERODYNAMIC  loading  terms  for  input  to  noise  calculations 


During  the  present  study,  considerable  effort  has  been  devoted  to  the  development  of 
simplified  methods  for  the  estimation  of  realistic  cirload  and  blade  motion  input  terms 
for  t he  acoustic  calculations.  The  basic  aim  of  this  work  was  to  devise  semiempirical 
techniques  to  extend  the  simplified  analytical  approach  of  Section  4.2  to  cover  the 
higher  loading  harmonics.  The  analyses  were  programmed  for  digital  computer  sol 
tion,  using  normal  mode  structural  response  theory  to  calculate  realistic  blade 
motions  in  all  degrees  of  freedom  (rom  the  estimated  loading  distributions. 

For  rwo  reasons  this  approach  had  to  be  abandoned.  First,  a  thorough  study  of  avail¬ 
able  experimental  airload  data^'*^  revealed  no  accurately  predictable  trends. 
Second,  as  the  ccouilic  study  progressed,  it  became  increasingly  clear  that  (1)  a  very 
large  number  of  loading  harmonics  are  required  for  the  calculation  of  a  moderate 
number  of  sound  harmonics,  end  (2/  *he  loading  fluctuations  must  become  increasingly 
random  as  frequency  increases.  Tills  is  to  be  expected  from  o  consideration  of  the 
turbulent  flow  processes  wh'ch  generate  these  higher  frequencies. 

Il  is  clear  that  new  information  is  required  before  firm  c^r.wlusions  can  be  drawn  re¬ 
garding  the  acoustically  important  characteristics  of  the  blade  loadings.  Continuous 
recordings  of  blade  section  pressure  differentials  should  be  acquired  by  a  system  having 
an  audio  frequency  range,  and  be  subjected  to  a  power  spectral  density  analysis, 
■rather  than  a  f  ourier.  analysis,  to  determine  not  only  the  load  harmonic  amplitudes 
but  dso  their  bn,.d  v  idths .  This  subject  is  discussed  further  In  Section  5.2  and  Appen¬ 
dix  if.  Nevertheless,  it  has  boon  possible  to  draw  tentative  conclusions  from  an 
examination  cf  die  first  ton  harmonic  amplitudes  experimentally  measured  and  re¬ 
port'd  in  Refer encos  21  and  34. 

Dat  from  these  references  are  presented  in  figures  19  and  20  which  show  a  selection 
of  85-percent-radius  section  loading  harmonic  amplitudes,  plotted  against  harmonic 
number,  on  logarithmic  scales.  These  plots  suggest  that  the  harmonic  decay  follows  a 
paw«r  law  of  the  form 

FX  *  Fo  X'k  (49) 


where  X  is  the  harmonic  number  and  =  y  L*  +  M*  is  the  section  harmonic 

loading  amplitude  .  F0  is  the  steady  component  of  the  section  loading.  Consequently, 
straight  lines  have  been  fitted  to  all  plots  for  the  available  forward  spead cases  (steady, 
level  flight),  enabling  the  loading  law  exponent  k  to  be  calculated  in  each  case. 
Figure  21  (a)  shows  the  exponent  values  plotted  against  advance  ratio  p.  A  line  has 
been  sketched  through  these  points,  but  in  view  of  the  scatter,  there  seems  to  be  little 
point  in  assuming  any  departure  from  the  value  k  •=  2  for  all  speeds.  The  most  sur¬ 
prising  feature  of  these  result?  is  the  relative  constancy  of  this  loading  law.  Despite 
expectations  of  an  increase  in  the  higher  harmonic  amplitudes  with  forward  speed,  no 
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The!”rm0nic  °"’P!ih'dK  «"•  practically  os  9reot  in  the  hove,  as 
they  c,e  at  the  limit, ng  advance  ratio  of  0.3.  To  study  this  phenomenon  in  mare  de- 

tail,  a  second  set  of  straight  lines  was  fitted  to  the  harmonic  amplitude  plots,  con¬ 
straining  rhe  lines  to  pass  through  the  second  harmonic  point,  which,  even  by  the 
Simple  theory,  is  predicted  to  rise- rapidly  with  forward  speed.  The  alternative  loading 
law  exponents  derived  in  this  manner  are  plotted  against  advance  ratio  in  Figure  21(b). 
More  variation  is  observed  this  time,  and  the  sketched  curve  bears  a  remarkable 
resemblar.ee  to  the  rotor  "power  required"  curve  (compare  with  Figure  16(a)).  The 
significance  of  this  has  not  been  established  and,  because  the  scatter  about  the  fitted 
straight  lines  was  somewhat  increased  by  this  procedure,  the  point  was  not  pursued, 
ft  should  be  noted  that  although  the  bulk  of  the  results  were  derived  from  the  H-34 
data  of  Reference  21,  four  cases  were  available  in  the  UH-1B  data  of  Reference  34. 
The  close  agreement  between  the  two  sets  of  data  is  of  considerable  interest,  since 
the  two  rotor  systems  are  very  different.  The  H-34  has  a  fou  -blade  articulated  rotor 
of  56-ft  diameter,  and  the  UH-1A  has  a  two-blade  teetering  rctor  of  44-ft  diameter. 


Figures  19  and  20  give  results  for  a  single  position  on  the  blade  at  0.85  span .  It  is 
of  interest  to  look  at  the  loading  variation  along  the  span  of  the  blade  which  is 
shown  in  rigure  22 .  Local  sectional  force  coefficients  in  the  various  harmonics  are 
plotted.  The  immediately  obvious  effect  in  Figure  22  is  that  the  higher  harmonic  load¬ 
ing  coefficients  rise  toward  the  tip  of  the  blade.  In  terms  of  power  laws,  the  results 
vary  from  an  inverse  square  law  based  on  second  harmonic  near  the  root  to  an  inverse 
first  power  law  at  the  0.95  span  location.  Based  on  the  steady  loading  coefficient, 
the  three  inboard  stations  correspond  to  about  an  inverse  cube  power  law,  while  the 
three  outboard  stations  correspond  to  un  inverse  square .  Since  the  outboard  stations 
are  much  more  effective  as  producers  of  sound,  this  effect  provides  additional  justi¬ 
fication  for  the  use  of  an  inverse  square  law  based  on  steady  loading  in  the  final  cal¬ 
culation  , 


The  fairly  wide  variation  in  harmonic  levels  should  again  be  noted.  For  instance,  all 
sections  have  second  harmonic  levels  close  to  0.003,  except  one  at  the  0.85  station 
which  shows  an  increase  by  a  factor  3 .  This  increase  in  second  harmonic  level  is  not 
observable  at  the  0.75  or  the  0.9  station,  and  it  is  further  evidence  of  the  highly 
localized  nature  of  the  loading  forces  even  for  the  lew  harmonics.  Similarly,  the 
0.9  station  shows  an  exceptionally  high  level  of  the  fifth  harmonic.  While  it  is  pos¬ 
sible  that  part  of  these  effects  may  be  explained  by  experimental  error,  it  seems 
probable  that  most  of  the  observed  effects  are  real .  For  instance,  the  very  high  level 
of  second  harmonic  is  observed  at  the  0.85  station  for  oil  the  low  forward  velocity 
cases  given  by  Scheiman;^  this  suggests  that  some  repeatable  vortex  pattern  is  the 
cause . 

It  is  also  of  interest  to  look  at  the  measured  loading  effects  for  some  rough  running 
cases  shown  in  Figure  23  in  comparison  with  the  hover  case  in  Figure  22.  Note  that 
the  rough  running  cases  have  higher  levels  of  the  higher  harmonics.  The 
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lougn  running  cases  a  so  generally  have  much  lower  second  harmonic  loadings,  which 
results  in  power  laws  based  on  the  second  harmonics  always  being  at  inverse  first  ,/ower 
law  or  higher.  The  power  laws  based  on  the  steady  loading  show  less  effect,  but  it  is 
thought  rhat  for  these  rough  running  cases  it  will  be  overconservative  tc  use  laws 
based  on  the  steady  loading.  On  physical  grounds,  it  is  clear  thot  for  cases  where 
the  shed  vortex  passe:  very  close  to  the  blade,  there  will  be  an  impulsive  loading 
applied,  which  will  tend  to  produce  all  harmonic  loadings  with  rouphly  equal 
magnitude.  It  is  tentatively  suggested  that  on  inverse  first  power  loading  law  be 
used  for  the  rough  running  cases  which  would  include  blade  slap. 

The  power  laws  discussed  above  appear  to  give  the  amplitude  of  the  loading  harmonics 
with  reasonable  accuracy.  However,  as  will  be  discussed  in  Section  5.2,  the  phase 
of  the  loading  is  equally,  or  perhaps  even  more,  important  in  determining  the  sound 
radiation  characteristics.  During  initial  work  on  this  project, attempts  were  made  to 
predict  the  phase  of  the  various  loading  harmonics  as  a  function  of  span.  This  was 
founc  to  be  impossible;  no  coherent  picture  emerged,  as  is  shown  by  typical  plots, 
again  from  Scheiman's  data,^  in  Figure  24.  Clearly,  the  best  assumption  is  of 
random  phase,  particular  jyar  the  higher  loading  frequencies.  This  introduces  some 
mathematical  complexity,  and  the  approach  which  has  been  used  is  described  in 
-Section  5.2.  ..  .  | 

Effect  of  Blade  Vortex  Interactions 

It  v.  j.  noted  in  the  Introduction  that  a  primary  source  of  fluctuating  airloads  on  the 
blades  is  their  passage  over,  or  even  through,  a  concentrated  vortex  trailing  from 
another  blade.  It  is  worthwhile  at  this  point  to  consider  the  form  of  the  rotor  wake 
geometry,  since  it  does  have  such  an  important  influence  on  the  noise  problem.  As 
discussed  in  Section  4.3,  the  vortex  sheet  shed  by  a  blade  is  unstable  and  rapidly 
rolls  up  into  at  least  one  concentrated  vortex  having  an  apparent  origin  near  the 
blade  tip.  In  fact,  its  precise  origin  varies  as  the  blade  rotates,  and  itdoes lie  a  little 
inboard  of  the  blade  tip.  However  for  th.  present  purposes  it  is  satisfactory  to  assume 
that  each  blade  trails  a  single  vortex  from  its  tip.  Furtf lermore ,  it  will  be  assumed  that 
each  element  cf  *he  vortex  remains  at  the 
precise  point  in  space  where  it  was  shed. 

The  plan  view  of  a  single  vortex  therefore 
will  look  something  like  the  accompanying 
sketch.  The  radial  station  s  of  the  nth 
following  blade  which  intersects  the  vor¬ 
tex  is  given  by  the  simultaneous  equations. 
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where  H*  =  2rn/B  +  tjj  —  , 

2  it  n 

and  — g —  =  the  angle  between  the  two  blades,  where  the 

symbol;  are  defined  in  the  sketch.  These  equations  were  solved  numerically  on  o 
computer  for  a  variety  of  rotor  configurations  and  advance  ratios,  and  the  results  ore 
shown  n  Figure  25,  The  curves  in  noch  case  represent  the  loci  of  biade/vortex 
intersection  points  on  a  single  blade  as  it  rotates.  The  number  on  each  locus  corresponds 
to  the  equency  of  preceding  biades.  For  example,  in  the  top  left-hand  diagram  for 
b  =  4,  j  =  0. 1,  the  number  1  corresponds  to  the  Intersection  point  of  the  vortex  shed 
by  the  blade  which  leads  the  reference  blade  by  90  degrees;  2  is  the  blade  leading 
by  180  degrees,  and  so  on.  Only  the  first  B  vortices  ure  sh  .»vn  for  clarity,  so  that 
the  fast  curve  in  each  case  represents  the  intersection  of  the  blade  with  its  own 
vo.tex  shed  earlier.  Also,  the  vertical  displacement  of  the  vortex  from  the  blade, 
which  of  course  is  not  shown,  will  generally  increase  as  the  number  increases,  so  that 
the  h  igher  numbered  intersections  will  have  a  decreasing  influence  on  the  blade  loads. 

Results  are  shown  for  rotors  having  4,  5,  and  6  blades  moving  at  advance  ratios  of 
0. 1  through  0.5.  However,  the  results  for  2  and  3  blades  can  be  seen  in  the  plots 
for  4  and  6  blades  as  alternate  curves.  This  is  the  reason  for  the  alternate  broken 
lines  in  those  cases. 


The  mcst  obvious  feature  of  these  results  is  the  general  movement  of  the 
intersection  points  to  the  region  of  the  disc  corresponding  to  low  values  of  ty  as 
forward  velocity  increase*.  More  important  than  this  from  the  standpoint  of  loading 
fluctuations,  however,  is  the  decrease  in  the  angles  between  the  blade  and  the  inter¬ 
section  loci.  Since  it  con  be  assumed  that,  near  any  intersection,  local  section 
loadings  will  be  strongly  influenced  by  the  presence  of  the  vortex,  these  loci  indicate 
the  me; ion  of  loading  peaks  (cr  troughs)  along  the  blade.  As  shown  by  the  general 
acoustic  Equation  (21),  aerodynamic  forces  in  motion  generate  sound  which  increases 
rapidly  with  velocity  and  acceleration.  At  low  advance  ratios,  the  curves  show  that 
the  motion  of  these  forces  is  relatively  low,  of  the  order  of  half  a  i-rdiui  in  half  a 
revolution.  This  corresponds  to  a  convection  Mach  number  along  the  liade  of  the 
order  of  0.2.  It  is  clear  that  as  advance  ratio  increases,  this  convection  speed  in¬ 
creases.  Typically,  at  p  =0.5,  the  intersections  move  half  a  radius  in  10  degrees 
of  rotation,  corresponding  to  a  convection  Moch  number  of  the  order  2.0. 
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Th  us,  this  fairly  crude  consideration  ;hows  through  Figure  25  that  even  at  moderate 
advance  ratios,  blade  vortex  intersections  can  travel  along  the  blade  at  supersonic 
trace  speeds.  Such  supersonically  moving  forces  are  efficient  generators  of  sound,  and 
it  is  interesting  to  consider  that  blade  slap,  for  example,  could  be  generated  by  this 
!  ind  of  mechanism.  Although  such  an  approach  has  not  been  taken  in  this  study,  it 
would  be  possible  to  perfornvan  analysis  of  the  sound  generated  by  a  blade/vortex 
intersection  in  arbitrary  motion. 

Physically,  of  course,  a  rapidly  moving  intersection  corresponds  toa  loading  fluctuation 
which  is  practically  !n  f  .use  along  the  entire  blade,  effectively  a  distributed  impul¬ 
sive  sound  source.  Viewed  in  either  manner  it  is  important  to  note  that  such  a  source 
will  be  highly  directional  with  a  marked  lobe  pointing  near  to  the  direction  of  the 
LlacU  axis.  Thus,  blade  slap  generated  by  this  mechanism  may  be  expected  to  peak  in 
azimuth  directions  within  the  first  quadrant  of  rotation. 

21 

Figure  26  shows  a  contour  plot  of  an  airload  distribution  measured  by  Scheiman 
for  a  four-blade  rotor  at  an  advance  ratio  of  0.29.  This  was  plotted  by  computer  using 
Scheiman’s  data  directly.  Superimposed  on  this  p*  ot  are  the  blade/vortex  intersection 
loci  for  B  =■  4,  p  =  0.29.  The  effects  of  the  vortices  ere  reflected  in  thq  loading 
gradients,  remembering  that  increased  angle  of  attack  s  to  be  expected  outboard  of  a 
vortex  location. 
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5.0 


COMPUTATIONAL  METHOD^  AND  ACCURACY 


5 . 1  SUMMARY  OF  COMPUTATIONAL  METHODS 

Two  ehgifxal  computer  programs  have  been  prepared  during  the  course  of  this  stud/  for 
the  CDC  3300  Computer  to  calculate  the  sound  field  cf  a  helicopter  rotor.  The  develop¬ 
ment  of  the  first  of  these,  code  named  HERON  1,  was  in  fact  one  of  the  primary  objec¬ 
tives  of  the  investigation.  This  program  is  completely  comprehensive  in  that  it  calculates 
the  sound  field  generated  simultaneously  by  an  arbitrary  number  of  rotors  of  any  dimensions. 
For  each  of  these  rotors,  the  program  accounts  for  the  aerodynamic  lift  and  drag  section 
loadings,  articulated  and  flexible  blade  motions  in  all  degrees  of  freedom,  rotor  geometry 
including  number  of  blades  and  rotor  diameter,  and  the  attitude  and  motion  of  the  rotor  in 
space ,  From  this  information,  the  program  computes  and  outputs  the  amplitudes  and  fre¬ 
quencies  of  the  fundamental  rotational  noise  component  and  its  harmonics  at  any  specific 
point  in  space.  Both  geometric  and  acoustic  near-field  effects  are  included  so  that  the 
result  is  accurate  at  any  point. 


Progrt  m  HERON  2,  on  the  other  hand,  was  written  specifically  for  the  direct  numerical 
evaluation  of  Equation  (37),  which,  as  discussed  in  Section  3.4,  is  a  closed-form  solution 
for  the  for-field  sound  of  a  rigid-rotor  system  whose  only  admissible  blade, motion  is 
sreadv  coning.  However,  as  will  be  shown  in  Section  6,  this  solution  is  satisfactory  for 
practical  application  in  the  study  of  helicopter  noise,  and  this  second  program  has  proved 
to  be  useful  in  aiding  a  general  understanding  of  an  extremely  complex  problem.  Unlike 
HERON  1,  this  program  was  written  entirely  as  a  special-purpose  mathematical  tool  and 
has  been  revised  and  refined  to  perform  a  wide  variety  of  computations  based  on  Equa¬ 
tion  (37).  Consequently,  no  general-purpose  version  was  developed . 


Program  HERON  'i  is  described  in  detail  in  the  companion  report.^®  A  complete  speci¬ 
fication  including  operating  instructions  is  included.  The  discussion  in  this  section  is 
limited  to  some  of  the  pertinent  fecte'es  of  the  computational  method  together  with  some 
notes  regarding  accuracy. 

The  acoustic  Equation  (20)  is  defined  in  a  coordinate  system  which  is  fixed  in  space 
with  a  distributed  acoustic  source  in  arbitrary  motion  relative  to  a  stationary  medium. 


p  = 


x.  -y. 

i  i 

a  r  (1  -M  ) 
0  r 


9 

3t 


where  q 


is  the  position  vector  of  the  aerodynamic  force  having  components 


F.. 

i 


The  seund  generated  by  a  rotor  blade  in  motion  is  the  result  of  that  motion  and  the  dis¬ 
tributed  aerodynamic  pressure  acting  over  its  entire  surface.  For  numerical  purposes, 
however,  it  is  necessary  to  simulate  the  actual  distributions  by  a  radially  distributed  set 
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of  point  loads,  as  is  ccmmcn  in  rotor  aerodynamics.  Further,  it  is  convenient  to  specify 
these  point  loads  as  having  orthogonal  lift  and  drag  components  normal  to  the  blade 
axis,  ond  acting  at  the  same  point.  When  one  knows  the  orientation  of  the  blade  *n 
space,  os  a  function  of  time,  it  then  becomes  a  relatively  straightforward  matter  to 
resoive  the  lift  and  drag  into  the  three  components  F.  defined  in  Equation  (20).  Strictly 
then  ,  Equation  (20)  should  be  written 


p--E 

4  IT 


xi-y; 


a0  r  (1  -Mr)  9t  \  r  (1  -Mf) 


where  the  tensor  subscript  applies  to  the  three  fixed  coordinate  directions,  and  the  sum¬ 
mation  is  carried  out  over  all  assumed  radial  loading  points  (and  over  all  blades).  In 
addition,  the  near-field  pressure  fluctuations  are  calculated  according  to  Equation  18  of 
Reference  8,  which  is 


P1  = 


F.  (x.  -y.)  (1-M2) 


(1-Mr)2  r2 


- F.  M. 

(1-M.)  '  ' 


To  evaluate  Equations  (51)  and  (52),  the  summations  have  to  be  performed  at  the  appro¬ 
priate  retarded  times  (which  the  brackets  denote).  This  is  simply  saying  that  we 
have  to  calculate  the  sound  pressure  at  the  observer  position  x.,  due  to  all  the  the 

aerodynamic  forces  F.  acting  at  the  point  Y.  on  the  rotor  system  when  they  generated 

the  acoustic  disturbances,  which  all  reach  the  observer  at  the  same  instant.  The 
accompanying  sketch  shows  clearly  that  the  sound  generated  by  a  number  of  points  on 

a  rotor  at  the  same  instant  does  not 
necessarily  reach  an  observer  at  the 

Nr  \  ‘ "  ■  .  """  *"  same  time.  The  arrows  denote  the 

W  - — — “  distance  travelled  by  the  sound 

originating  from  each  point  when 
the  sounc*  reaches  the  observer. 
s''  Figure  27  shows  the  computed  posi- 

I  Nn\Tx^/'  tions  of  the  blade  axes  of  a  four- 

blade  rotor  corresponding  to  the 
appropriate  retarded  time  of  each 

blade  element.  The  sound  generated  by  each  element,  when  it  was  in  the  position 
denoted,  arrived  at  the  specified  observer  position  at  the  same  instant. 

The  observer  positions  and  the  rigid  body  motions  of  the  helicopter  (in  terms  of  its  posi¬ 
tion,  and  linear  and  angular  v  elocities)  are  definedwith  respect  toa  set  of  "ground  axes" 
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which  cc .  fixed  in  space.  For  convenience,  al!  blade  loadings  and  motions  are  speci¬ 
fied  wit,  respect  to  a  set  of  rotor  axes  which  are  parallel  and  normal  to  the  rotor  plane 
of  rotation  and  whose  origin  is  coincident  with  and  moves  with  the  rotor  hub.  The 
.  azimuthcl  variations  of  all  section  airloads  and  motions  are  input  to  the  program  either 
as  time  histones  or  as  Sets  of  Fourier  coefficients.  By  performing  the  appropriate  axis 
transformations,  the  coordinates  of  each  blade  element  and  the  force  components  acting 
upon  it  can  be  defined  with  reference  to  the  ground  axes  in  a  form  suitable  for  direct 
input  to  equations  (51)  and  (52)  as  a  function  of  time.  All  that  is  required  to  pejform 
the  summations  involved  in  those  equations  is  the  correct  retarded  time  for  each  blade 
element. 


This  is  calculated  by  an  iterative  procedure.  For  uniform  linear  motion  of  the  rotor, 
‘the  retarded  time  for  the  hub  itself  can  be  obtained  from  the  equation 


t  -  t'  = 


(x.  -y.)  y.  +  yf](x.  -y.)  y. } 2  +  r2  (a2  -y2) 


(53) 


a 2  -  <,* 

0  Xi 


where  tha  y.  are  the  coordinates  of  the  rotor  hub  at  time  t .  This  gives  a  good  approxi¬ 
mation  to  the  retarded  time  for  the  most  inboard  blade  station,  which  is  then  found  by 
iteration .  The  converged  value  is  used  as  a  starting  value  for  the  next  blade  station, 
and  so  on,  until  all  have  been  determined. 


;The  sound  pressure  at  each  observer  position  is  calculated  as  a  time  history  at  a  series 
of  equal  dme  intervals  over  a  total  period  equal  to  the  rotor  rotational  period  divided 
by  the  number  of  blades.  This  is  the  period  of  the  fundamental  sound  harmonic. 

Starting  et  time  t  =  0,  the  retarded  times  are  calculated  for  all  blade  loading  stations 
on  all  blcdes,  and  Equations  (51)  and  (52)  are  evaluated  to  give  the  sound  pressure 
p(0).  The  process  is  then  repeated  for  t  =  At,  2  At...  ,  and  so  on,  until  the  complete 
time  histery  is  obtained.  This  is  then  Fourier  analyzed  to  give  the  amplitudes  of  each 
sound  harmonic . 

Numericol  Accuracy 

The  implications  of  simulating  the  distributed  aerodynamic  loadings  by  a  set  of  discrete 
point  forces  are  fully  discussed  in  Section  5.2.  From  that  point  on,  the  computed  results 
are  exact  within  the  numerical  accuracy  of  the  computer  For  sound  level  calculations, 
this  is  sor  ev/hat  limiting  due  to  the  very  large  dynamic  range  of  the  hearing  mechanism 
pf  around  140  dB  in  sound  intensity.  Due  to  the  limitations  of  the  computer  word 
length  in  ^he  machine  used,  the  useful  dynamic  range  is  only  around  90  dB.  However, 
for  most  p-actical  purposes,  this  is  more  than  adequate. 

Another  question  of  accuracy  arises  in  the  specification  of  the  airload  data.  The  num¬ 
ber  of  sound  harmonics  computed  by  the  program  is  approximately  m/2B  where  m  is 
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the  number  of  azimuth  stations  used  to  define  the  azimuth  variation  of  loading  and 
motion  and  B  is  the  number  of  blades  in  the  rotor.  Since  m  values  are  sufficient 
to  accommodate  approximately  m/2  loading  harmonics,  the  program  effectively  cal¬ 
culates  l/B  times  as  many  sound  harmonics  as  loading  harmonics.  It  will  be  seen  in 
Section  6  that  a  certain  number  of  loading  harmonics  can  lead  to  the  accurate  com¬ 
putation  or  somewhat  less  than  l/B  times  as  many  sound  hcrmonics;  therefore,  the  possi¬ 
bility  of  error  :n  the  highest  harmonics  must  always  be  borne  in  mind. 

Although  the  program  contains  provisions  for  helicopter  pitch  and  roll  angular  veloci¬ 
ties,  it  is  assumed  that  these  do  not  result  in  any  change  of  flight  direction  during  the 
entire  ange  of  retarded  time,  which  is  approximately  equal  to  l/B  times  e  totor 
rotational  period.  This  assumption  avoids  very  considerable  computational  Complexity 
and  is  believed  to  involve  very  small  errors.  1, 

5.2  ACCURACY  CONSIDERATIONS  -  EFFECTS  OF  RANDOM  LOADINGS 

One  of  Hie  principal  requirements  in  any  numerical  study  is,  of  course,  thaf  accuracy 
be  obtained  in  the  answers.  In  the  present  approach,  theavailability  of  two, independent 
methods  for  calculating  the  noise  has  acted  os  a  powerful  check  against  computational 
errors,  but  several  extremely  important  and  rather  subtle. points  affect  eithet  method 
equally .  In  fact,  it  will  be  shown  that  the  la.  k  of  detailed  experimental  knowledge  of 
the  higher  harmonics  of  helicopter  noise  can  ejsily  lead  to  the  introduction  of  errors. 
These  errors  result  from  the  nature  of  the  assumptions  made  to  cover  the  lactc  of  experi¬ 
mental  data  and  are  by  no  means  obvious. 

It  was  shown  in  the  previous  section  how  actual  blade  lift  and  drag  loading  distribu¬ 
tions  are  represented,  for  the  purpose  of  numerical  calculation  of  the  sound’field,  by 
a  finite  number  of  discrete  point  forces.  At  this  point  the  question  arises:  how,  pre¬ 
cisely.  is  this  simulation  performed  in  terms  of  numbers  and  distributions  of  the  point 
loads? 

If  we  knew  the  actual  pressure  distributions  on  the  blade  surfaces  as  a  function  of  time 
or  blade  position,  with  sufficient  accuracy,  it  is  clear  that  the  most  accurate 
solution  would  be  to  use  a  very  large  number  of  loading  points,  both  radially  and 
around  the  azimuth  (or  in  the  case  of  harmonic  representation,  the  maximum  number 
of  Fourier  coefficients  for  euch  radial  position).  It  is  equally  clear,  however,  that 
(a)  we  cannot  define  the  actual  blade  loadings  with  sufficient  accuracy  to  justify 
more  than  a  certain  number  of  radial  and  azimuthal  intervals  (note  that  the  use  of 
1 -degree  and  2-degree  azimuth  intervals  in  previous  investigations  cannot  improve  the 
accuracy  over  the  15-degree  intervals  for  which  the  loadings  were  originally  specified; 
the  higher  harmonics  yielded  by  this  technique  are  in  error),and  (b)  it  is  necessary  to 
minimize  the  number  of  points  for  computational  expediency.  In  the  preseqt  program 
the  azimuthal  interval  defines  how  many  sound  harmonics  we  can  calculate;  whereas 
the  radial  intervals  govern  the  accuracy  with  which  they  are  calculated.  Therefore, 
for  a  given  number  of  sound  harmonics,  the  required  azimuth  Interval  (or  number 
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of  locding  harmonics)  is  fixed,  and  it  only  remains  to  optimize  the  number  and  distribu¬ 
tion  of  radial  loading  points  in  order  to  obtain  adequate  answers.  Let  us  consider,  then, 
some  of  the  relationships  between  the  loading  distributions  and  the  sound  which  they 
generate . 


Toking  one  pair  of  results  from  Equations  (37)  (specifically  for  n  -  X  even),  and 
dividing  the  cosine  harmonic  by  the  sine  harmonic  we  obtain  the  following  result 


(54) 


n  is  the  sound  harmonic  number  and  X  is  the  loading  Harmonic  number  (where  the  rota¬ 
tional  frequency  is  the  fundamental).  The  coefficients  and  b^  are  the  in-phase 

and  quadrature  amplitudes  of  the  nth  sound  haimonic  that  are  attributable  to  the 
Xth  loading  harmonic  .  The  X  th  loading  harmonic  is  represented  by  the  coefficients 
a  ,  L\..,  etc.,  which  correspond  to  the  harmonics  of  a  point  force  acting  somewhere 

on  the  blade.  The  radial  location  of  this  force  affects  only  the  value  of  M,  the  rota¬ 
tional  Mach  number,  which  is  directly  proportional  to  radius.  Equation  (54)  gives 
the  phasing  of  the  sound  harmonic  since  tan  <J>nXR  -  (anx/bnx'R  '  where  lhe  SUff,X 
R  applies  to  the  force  acting  at  radius  R. 

We  ee  from  this  equation  that  if  only  one  force  is  acting  (having  thrust,  drag  and 
radial  components),  the  phase  of  the  sound  harmonic  is  essentially  independent  of  the 
point  of  action  of  the  load,  apart  from  a  minor  variation  due  to  the  effect  of  M  on 
the  drag  term .  Similarly,  if  many  loads,  acting  on  the  blade  at  different  points,  are 
phase  with  each  other  (i.e.,  a/b^  =  constant),  men  their  contr.but.ons  to  ony 
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sound  Harmonic  ore  cis©  in  phase,,  and  w*:  have  s*7npty 

a 


rX  •  &  °r,X*  ''"  %\R.  ■*  °nXR.  “VV,';  7v>~ 

lOiAw  i  i  •».•.  ,  ..  - _ 


1  'TO *  A !  ^ i</w*  i 


;(55) 


tn  "‘-.V  VO-  -'  V'^rtO 


V.X 


where  R.  are  the  loading  points,  and  the  ratio  <r” 
i  a 


TOTAL  '  Lu 


pXR. 


nX  .TOTAL  . 


■  for  -ali  i . 


•b 


nXR. 


However,  if  the  various  toacs  are  not  in  phase  with  each  other,  or  if  the  ratio? 

°XT;°XO;qXC  and  bXT:hXDibXC  0fenetcr'  ♦ont/theh  the  phases 

vary  with  i .  '*'"■■■■  /■’ 

The  consequence  of  this  is  that  *f  the  harmonic  loading -is  in  ptiase  along  ,+he  blade,  it 
can  be  quite  accurately  represented  by  o  single  loading  point  if  the  correct  forces  (fa- 
example,  the  integrated  btode  thrust,  drag,  and  outward  component  harmonics)  are 
applied  at  the  qpirect  point.  This  is  precisely  what  hnj  been  Oyne  in  accepted  propel¬ 
ler  noise  theory,  •  and  it  allows  considerable  simplification  of  the  confutations.  It 
is  important  to  note,  though,  that  the  effective  loading  point  iso  function  of  X  ahd  .n, 
and  thus  it  takes  on  many  different  vdlues  .In  g-camplete  calculation  -for  oil  X  and  n’-- 
However,  if  the  iooding  is  not  in  phese,  then  an  effective  single  food  would  hove  to'be 
specified  in  terms  of  position  and  phase  fa  each  field  point,  which  is  ooviously  imprac 
tied .  .  . ' 

'  •  -w  -  - 

V.' 

Unfortunately,  os  we  have  seen  in  the  case  of  the  helicopter,  it  is  on(y  the  steady 
loading  which  remains  in  phase  oieng  the  blade  (by  definition),  the  phase  of  the  har¬ 
monics  showing  increasingly  greater  variations  as  frequency  increases,  'As  a  general  - 
role,  it  seems  safe  to  soy  that  the  higher  the  Iooding  harmonic,  the  mote  rodioJly  dis¬ 
tributed  the  point  loads  must  be  fa  accurate  sound  calculation.  It  is  equally  certain 
that,  under  the  present  state  of  the  art,  we  cannot  define  the  loading  harmonics 
with  sufficient  resolution  to  satisfy  the  requirements  of  the  acoustic  theory,  since  the 
final  results  are  extremely  sensitive  to  loading  phase,  as  will  be  seen  in  Section  6. 


Fortunately ,  due  to  the  nature  of  rotor  aerodynamic  loads,  we  are  oble  to  simplify  the 
problem  in  a  manner  which  enables  us  to  make  reasonable  estimates  of  the  sound  field 
up  to  very  high  harmonics  of  the  blode  passage  frequency.  The  basis  for  this  simpli¬ 
fication  is  that  the  flow  processes  which  cause  the  fluctuating  airloads  become 
increasingly  random  as  frequency  increases.  This  Is  appcrer.r  from  experimental  date 
on  the  first  ten  loading  harmonics  (Section  4.4)  and  is  clearly  true  at  frequencies  of 


the  wake  turbulence.  Experimentally,  this  randomness shows  itself  os  a  iock  of  repet¬ 
ition  of  results,  either  in  successis  e  experiments  or  in  successive  felode  revolutions, 
;^cx#xcyjip]e  i  i£_dr>  el  ectnC.o.L  signal,  which  is  proportional  to  the  blade  section 
"loodthg  of  any  station,- were  passed  through  a  very  narrow  hirer  to  exiiuui  sorae  par- 
ti.puAQr  hofmoolc  bevel,  -the  output  woo  to.  be  seen  to  nuuiuuie  in  amplitude  and  phase . 


To  uodersfond  me  implications  ot  a  random  loading  variation  along  the  Place, 
r*k» -integration. ot_tr»e-curwa.ahui*fj Sst Jus. ikc-ah  Jseiqr.v  v.fsich  {(«  be  thought 


VejpreVentlog'ths  .acoustic" contributions  from  each  element  of  the  rotor  blade. 


k“.IAit  l>ay«  the  some  root  mean  square  value  '(H)  over  the  range  of  interest .  Consider 
'p  first  the  simple  straighthUne  case  (A) .  Suppose  that  this  is  divided  into  (soy)  100  parts 
#och  of  iength'i^/lQOl  Then  each  parr  will  contribute  one  100th  of  the  total  integral. 
Now  in  sound  calculation^, we  must take  the  square  of  the  result.  If  we  odd  each  con- 
tribgtion  and  then  square,  the  result  will  be  ;  but  if  we  square  each  contribution 
and  then  add,  the  result  wi  II  fee  R2  bP/lOO,  In  other  words,  the  value  of  the  sum  of  the 
squares  pan  b*»  mode  qs  small  as  we  like  by  taking  more  and  more  elements. 


Next  qanslder  the  curve  B,  which  is  a  perfect  sine  wave.  Divide  this  into  100 
pieces,  and  qd«U  The  result  is  identically  «ro;  but  if  we  square  each  contribution 
and  then  qdd,  thereswlt  is  again  R*  H^/IOO.  Thus  in  this  case,  the  result  of  taking  the 
sum  of  the  jqycres  is  always  yeater  than  taking  the  square  of  the  sum .  The  difference 
lies  in  the  relative  phasing  along  the  length  of  the  curve.  The  straight 
line  was  all  in, phase,  so  each  segment  odded  to  its  nei^rbor.  For  the  sine  wave,  alt 
contributions  cancelled  exactly.  Curve  C  shows  a  further  case  where  the  sine  wove 
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ls  inc^ple-e.  Clear  I  v,  the  sum  of  the  square*  will  again  be  closely  eauai  to-R7-^/^ 

-while  .V  a  direct  summation,,  the  whole  contribution  O*  curve  «.  nui  wgpngrn**, 

the  last  wiconcelleb  pert  angle.  --  *  -  -- ---  ; — 

,  <  ■  .u..  ~r-~  r*-nh!e m.  Suqdoso- a  load 

Tk.««  cornu  cere; ions  sie-unr  siivnwv  * - -  ••  c: 

harmonic  i*  taken  to  be  in  phase  over  the  whole  blade  span .  because  or  atmp^ 
time  effects,  radiation  trotr.  different  parts  or  the  span  will  reoon  ^e  oOserver^^^ 
"ferent  "times"]  Thus';  the  in-phase  'ioodir.g-ir.jw:  azzsxsi 

In  other  wards,  curve,  like  Acre  transformed  into  curves  Uke  S  or  C  vas  ^ 

retaraed  time  effects.  Fortunately,  however,  in  the  integrated  result 
there  is  no  phasing  dependence  on  radius  due  to  the  acoustics  Bu,  the  loodmg  phase 
effects  are  still  of  extreme  importance.  It  is  well  Known  ,see  Section  2.  _  mt* 

blade  loadings  can  be  verv  localized,  being  the  resuit  of  vortex  interaction  e  fee*, 
so  that  it  is  extremely  unlikely  that  they  are  in  phase  over  the  span  -  The  mart  con¬ 
venient  assumption  is  that  the  higher  harmonic  loads  ore  random  in  phase °}^S 
blade  as  was  discusses  above .  Jr.  this  case,  the  sum  of  the  squares  calculation  applies,. 
However,  here  another  problem  arises  .  It  can  be  seen  that  the  sum  of  the  squares  con  be 
made  as  small  as  desired,  simply  by  choosing  a  sufficiently  large  number  o 
olong  the  span.  Again,  this  is  obviously  incorrect.  The  fallacy  behind  this  hm.t  .s 
that,  as  the  divisions  get  smeller  and  smeller,  it  becomes  less  and  less  accurate  to 
assume  that  the  phases  of  successive  divisions  ore  random .  The  phases  of  nearby  .vi¬ 
sions  are  approximately  equal,  so  that  they  must  be  calculated  via  the  square  of  the 

sum  rather  than  the  sum  of  the  squares.  The  key  question  is:  how  many  successive 

divisions  con  be  regarded  as  being  in  phase? 

The  answer  to  this  question  lies  in  The  space  correlation  of  the  fluctuating  loads  along 

Hie  blode .  Consider  the  in-phose  and  out-of-phase  combinations  to  the  nth  sound 

harmonic  ^  o  and  b  bn  ,  b  . .  due  to  K  Wn, 

points  chosen  to  act  at  various  blade  stations  cn3  representing  the  distribution  of  rise 
X  th  loading  harmonic  .  If  these  K  loads  are  constant  in  both  amplitude  and  phase, 
so  are  the  a.,  b  .  (i  =  1 ,  2,  3, . . .  K) ,  end  we  can  writ,  that  sound  pressure  level 


n  i 


b  .  (i  =  1,  2,3,...  K) 

ni 


of  the  nth  sound  harmonic  as 


.2  _ 


n. 


K 

E 


Now  if  the  K  loads  are  random,  i  .e. ,  not  constant  in  amplitude  and  phase,  the  above 
result  is  not  true,  and  we  must  write 
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>if,; 


-K  K  f 


^  --- 
n 


.  ^  {  r..  w  a  t-  j,,  t  b  \ 

\  !j  n.  sj  *.  n.  J 


■a-cdKhfeiSfc  31  v  vffe  thr  c.vfe’siiort  cOvffir.isr.t;  which  dsfinsthe  decree  af 

—  -  .;-.,tj..-.  ij  ■  ' 

•randomness  between  <»  and  a  or  b  ond  i.  ,  Thesu  coefficients  have  values 
—  •‘•v  —  n*  n,  *n;  n- 

■■■>'  .  i  .  j 

;*  ranging  between  ^erri-and  unity,  where  the  former  Implies  complete  randomness  and 

7 1th* ijgrtar/at  follows  from  equating  Equations  (56)  and  (57),  implies  complete 
^v'icarTeictfqn.  If  fhe-K  loads  are  completely  random,  r. .  =  C  (except  when  t  =  j, 
*in§e.r;,.  Js  qlyyoy*  unity)  and  then  _  *■* 


■Ek-< 


.:r7ln7gjeneral ,  the  coefficient*  v*i  If  fie  somewhere  between  the  two  extremes  ond  will  de- 
'cfeos*. «"7'|i  /-  }j  increases*  In  order  to  make  use  of  Equation  (57),  it  is  necessary 

to  deflne.the  T. .  ond  .7,.  .  Lin  doing  so  we  (poke  the  following  assumptions: 

•”  yk- ,  r-y  ..  .  ‘  .  %. 

7  .  "V  -fy  -That  T.  -  T.. 

^  -j  ■  _ 

L 2).  'l;Thdrt  7  decrease*  exponentially  with  the  separation  |  (expressed  as  a 
;  -  -fraction  of  the  redius)  between  the  ith  and  jfh  looding  stations. 

_  '  3)  That  the  "cotrelctian  length11,  defined  here  as  the  distance  over  which  the 

.7  7‘  '  'correlation  coefficient  foils  to  some  fixed  value,  is  inversely  proportional  to 

harmonic  number  X. 

Using  these  assumptions,  we  arrive  at  the  result 

7..  =  T..  -  «"°X«  (55 

'  .  ‘J  ‘J 


where  _o  is  a  constant  to  be  determined.  The  variation  of  7..  with  £  and  X  is 
'sketched  on  the  following  page.  ‘-I 

In  order  to  use  Equations  ( 57)  ond  (59)  in  c  numerical  procedure,  a  further  calcula¬ 
tion  is  required .  If  the  correlation  coetticienf  does  not  differ  substantially  between 
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successive  loading  points,  then  Equation  (57)  probably  gives  sufficient  accuracy..  if 
the  coefficient  does  change  ilyr.Jfrcsxt!/,  -  -  . 

however,  the  ronuvut  chcnsw  that  $cew  „  *  ■■„. 

over  c  sir.gie  loading  segnve-it  may  be  t 

cApecivC  ic  affect  the  results,.  since  the  '  ^  .  *  ~  *  j 

effective  load  is  somewhat  iower  than  LXv\x"V%^.  *  _ I~J 

the  assumed  value.  This  is  clarified  be-  5  VOO'v.  ~  1 


Considei  o  single  rodia!  segment  of  the 
rotor  blade  between  the  stations  S;  one 

S  .  if  the  section  loading  is  well  corre- 
2 

iated  along  its  length  then  the  harmonic 
amplitude  of  the  total  force  acting  upon 
it  is  simply 


iF>  ~  { 


FI  (slds 

A 


\;hs'r-  F^(s)  is  dF^/ds,  the  amplitude 

of  the  X  th  harmonic  section  iooding.  If 
F^  (s)  is  not  well  correlated,  then  we  must 
resort  to  the  method  illustrated  by  Equa¬ 
tion  (57)  to  determine  a  vaiue  of  an 
"effective"  load  increment  f  acting  on 

EX 

the  segment.  In  fact,  we  can  write 


w 


f S2  f  S2 

=  J  J  r  12  FX  ^Fx  Wds  * 


Now  assuming  the  same  form  for  the  correlation  coefficient 


=  e-0  ^  ^  ,  where  £  = 


and  assuming  a  "triangular  loading"  pottern  (i.e.,  the  local  loading  is  proportional  to 
radial  distance) for  the  Xth  harmonic  amplitude  so  thef 


Fx  <s>  *  2FX  •  = 


we  CO;'  c  , press  the  following  relationship  for  the  effective  load  Fp  . 


:  -  4F*  P  f  f‘ 

■x  x  A.  i  A. 


—  -aX(?  -  s)  ,  f  ?  _  -aX($  -?)  , 

s  •  s  ■  e  ds+l  s  •  s  •  e  as 


ds 


(62) 


Evaluation  of  this  integral  gives 
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;X  aX 
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aX 


S2  +  S2  - 

'  2  a2  X2 


(63) 


It  is  of  interest  to  put  =  0  and  =  1,  for  rr.en  we  obtain  the  effective  loading 
.for  the  entire  blade 


r- 


4FX 


^X  aX 

Note  <Ho.'  for  a  X  »  1 


’  2 

1  2 

6  -  e“°XNi 

3 

aX  a3  X3 

V  c  / 

1 

3  F2 

1 

3a  X 

(64) 


(65, 


Combining  Equation  (65)  with  Equation  (49)  and  using  the  recommended  value  of 
2  for  k,  we  have  the  approximate  result 
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(where  Tg 


tots’  blode  fh-rv.  ii. 


Thof  is,  If  the  tooas  ore  ccTt-cUy  fetorg  in  mOgnitude  as  X  %  then  the  effective  power-’' 

■  •■ -■'a4  Uu*  i  i  r-  /?x.ri  i  r*  tn  b  j»fv.T»>yi»n*rAw't  i-T"-**-- 


Equations  (63)  and  (66)  have  been  used  fc  study  the  effects  of  blade  loading  da- ' 
tribofionsy  and  the  results  are  discussed  in  Section  6. 


A.  r-:on<2  equivalent  error  wobiera  occurs  in  the  calculation  of  the  sound  from  "he 
summation  of  oil  the  loading  harmonics,  tleferer,-*  io  Equation. (36)  or  {37)  shows* 
that  the  contributions  to  any  sound  harmonic  of  every  second  loading  harmonic  wilt  • 
have  c  ;>gn  change.  Now,  by  the  time  the  loadings  reach  the  twentieth  or  so ,  their 
magnitude  changes  little  with  change  in  harmonic  number,  'thus,  ail  the  sound  from, 
the  loading  harmonic  would  be  cancelled  out  by  the  sound  from  the  next  if. they  both 
hod  the  same  phase  angle.  Clearly,  it  is  extremely  unlikely  that  such  harmonics 
would  have  the  same  phase,  but,  as  shw.i  In  Section  4r  beyond  about  the  second 
harmonic  we  have  no  real  Idea  of  what  the  phase  is.  Again,  the  only  logical  choice 
is  to  assume  that  the  harmonics  have  random  phase,  so  that  the  sum  of  the  squares  of 
successive  harmonics  must  be  taken  rather  than  the  square  of  the  sum.  If  realistic 
answers  are  to  be  expected.  On  the  other  hand,  if  the  phase  is  weft  kr.c.  n  (a*  for  - 
instance  on  the  first  harmonic  airiocd),  then  it  con  be  included,  and  the  second  expres¬ 
sion  can  be  used.  It  may  be  noted  thct  the  definition  of  the  loading  harmonics  via  a 
roof  mean  square  value  automatically  takes  care  of  these  problems  In  the  azimuthal 
integration.  These  points  were  not  considered  by  previous  investigators ,  *^r  12  ^,j4 

is  another  source  of  the  low  levels  they  find  in  the  highe'  noise  harmonics,  Ali  the 
points  mode  above  via  physical  argument  can  be  duplicated  mathematically .  A 
mathematical  analysis  is  giver,  in  Appendix  II. 


As  a  final  point,  it  is  worthwhile  to  consider  the  chordwise  integration.  The  genera! 
conclusions  obcve  obout  correlation  functions  for  the  high  harmonics  do  apply  in  a 
modified  form,  but  the  abject  here  is  to  study  the  direct  effects  of  a  chordwise  distri¬ 
bution  of  loading.  Schieget  ef  ai,  foiiwed  Gutin  **  in  using  a  rectangular  dis¬ 
tribution,  and  Loewy  and  Sutton  used  a  straight-line  approximation  which  gave 
a  complex  algebraic  result.  In  fact,  these  modifications  are  both  unnecessary  and 
inaccurate.  Figure  28,  token  from  Reference  35,  shows  the  Fourier  coefficient*  of 
various  chordwise  loading  shapes,  a  typical  reai  distribution  (from  Reference  36),  and 
assumed  rectangular  and  sinusoidal  shapes.  Several  points  may  be  observed.  First, 
the  lock  of  certain  harmonics  is  peculiar  to  the  artificial  shapes  and  does  not  occur 
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forth  real  case,  5econd,  the  real  loading  gives  higher  levels  of  harmonics  than 
do  tb<  artificial  shapes.  This  is  basically  due  tc  the  strong  first  peak  of  the  real 
loading.  Third,  solidities  typical  of  a  helicopter  have  much  higher  harmonic  con¬ 
tent  than  equivalent  propeller  cases.  Fourth,  the  delta  function  approximation 
(aqua-  harmonic  content)  is  very  acceptable  for  helicopter  loodings  and  is  certainly 
more  accurate  than  rectangular  distributions.  {The  results  shown  in  Figure  35  are 
not,  of  course,  noise  harmonics.  The  effective  values  must  be  taken  through  the  full 
computation  of  Fauation  (36)  to  give  the  final  noise  radiation.)  !t  is  clear  that 
detai’,  of  the  mean  chordwise  loadings  distribution  are  unlikely  to  be  important. 

This  result  was  found  in  the  computations  of  Loewy  and  Sutton. ^0,11  Thus,  it  is  sug¬ 
gested  that  point  chordwise  loadings  (delta  functions)  be  used  in  helicopter  noise 
calco  a1>ons,  and  this  has  been  followed  throughout  the  present  work. 

It  is  quite  probable  (hot  the  use  of  the  point  loading  approximation  is  conserva¬ 
tive,  In  that  it  will  generally  give  higher  harmonic  noise  levels  than  any  other  chord- 
wise  loading  distribution.  The  overall  sectional  lift  coefficients  used  in  the  present 
work  ire  derived  from  Scheiman’s^  V  data,  as  discussed  in  Section  4.4.  The  data  are 
the  result  of  averaged  pressure  readings  over  the  blaue.  Examination  of  raw  pressure 
data  fas,  for  instance,  sketched  in  Cox  and  Lynn's  report^),  inevitably  shows  many 
local  fluctuating  pressure  peaks.  These  pressure  peaks  will  radiate  noise  efficiently 
at  the  high  frequencies.  Alternatively,  the  peaks  would  be  reflected  in  an  increase 
of  the  Fourier  amplitude  coefficients  with  an  increase  in  frequency  in  Figure  28. 
Watkins  and  DurKng37  showed  a  similar  effect,  in  that  any  part  of  the  blade  with  a 
Si'all  wading  pattern  will  also  possess  high  levels  of  the  high  Fourier  coefficients.  It 
mighl  also  be  noted  that  the  smoothing  operations  used  by  Scheiman7'  could  be  a 
sourc  of  considerable  error  in  estimating  the  levels  of  the  harmonic  airloads.  This 
poin^  certainly  justifies  further  study. 
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6.0 


RESULTS  AND  DISCUSSION 


6 . 1  PARAMETER  STUDY 

The  .!  sbal  requirement  in  the  anclysis  of  the  results  Is  to  determine  the  general  features 
of  the  ncse  radiation  from  the  helicopter.  This  has  been  achieved  by  a  thorough  study 
of  the  results  using  Equation  (37)  of  this  report  (Computer  Program  HERON  2) .  The  more 
detail  -:d  study  of  near  field  and  blade  motion  effects  was  accomplished  by  using  the 
exact  computer  program  (HERON  1)  via  Equation  (20)  of  this  report. 


Acous  ic  Effectiveness  of  Individual  Loading  Harmonics 

First  cf  all,  the  effectiveness  of  the  various  loading  harmonics  as  sound  sources  will  be 
evaluated .  Figures  29  and  30  give  the  sound  radiated  in  various  harmonics  of  a  four- 
blade  otor  by  various  loading  harmonics.  The  magnitude  of  each  loading  harmonic 
(including  the  steady)  has  been  assumed  to  be  the  same,  and  in  each  cose  the  load  has 
been  educed  to  a  single  point  fluctuation  force  acting  at  0.8  of  the  blade  span.  The 
rotatb  nal  Mach  number  of  this  point  is  0.5,  so  that  t!;?  graph  corresponds  to  a  tip  Mach 
numbc  of  0.625.  The  magnitude  of  the  thrust  component  has  been  assumed  to  be  equal 
to  10  ‘irr.es  that  of  both  the  drag  and  outward  components  of  force,  which  are  assumed 
to  be  quol .  This  proportionality  is  typical  of  a  helicopter  rotor.  Figure  29  gives 
curve  corresponding  to  o  field  point  10  degrees  below  the  rotor  disc,  and  Figure  30 
gives  urves  for  a  field  point  10  degrees  from  the  rotor  axis.  The  same  general  features 
may  b  •  observed  in  each  of  the  harmonics  plotted .  Only  o  limited  range  of  loading 
harmc  ,ics  contributes  to  each  noise  harmonic.  Take,  for  ir.rtance,  the  fourth  sound 
harmonic  (m  =  4)  in  Figure  29.  Loading  fiarmonics  below  the  eighth  can  be  seen  to 
produce  little  noise.  Between  the  eighth  and  twenty-fourth,  the  sound  produced  varies 
but  is  of  roughly  the  same  order  of  magnitude,  while  beyond  the  twenty-fourth  loading 
harmonic  the  sound  radiation  falls  away  rapidly.  Thus,  it  may  be  concluded  that  on  a 
real  helicopter  rotor,  when  all  loading  harmonics  can  contribute  to  the  observed  noise, 
loading  harmonics  be'ween  the  eighth  and  twenty-fourth  must  be  included  to  obtain  an 
accurcie  calculation  of  the  fourth  harmonic  .  This  conclusion  explains  why  Schlegel 
et  al  .  calculate  such  low  levels  for  the  fourth  harmonic  in  their  report  12  (see  also 
Figure  12),  since  they  include  only  loading  harmonics  up  to  the  tenth. 


This  e  rect  .nay  be  understood  by  reference  to  Equation  (36).  There  are  two  basic  terms 

in  this  equation,  the  J  .  and  the  J  .  Fc,  the  values  of  the  argument  of  the  Bessel 
'  n-X  n+A  57 


function  typical  of  the  helicopter  problem,  the  J  .  terms  wilt  be  quite  insignificant 

•'A 

compo.ed  to  the  J  .  terms  and  ccn  up  ignored  (except  oerhaps  for  the  low  harmonic 
n-A 


noise)  .  Now  n  -  mB  where  m  is  the  harmonic  order  and  B  is  the  number  of  blades. 
Thus,  us  the  loading  harmonic  X  increases,  the  order  of  the  Bessel  function  n-X  will 
decrease,  eventually  going  negative.  For  example,  for  m  ~  4,  B  =  4;  n-X  =  6  for 
X  -  10  (the  tenth  loading  harmonic);  and  n-X  -  -6  for  X  =  22 .  The  absolute  value  of  a 
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Bessel  functi  on  with  a  negative  order  is  equal  to  the  value  of  the  Bessel  function  with  a 
positive  ord-  r.  Thus,  we  expect  the  results  for  the  X  =  10  and  22  cases  to  be  the  same 
for  the  fovrtr  harmonic  case,  inspection  of  Figure  29  will  show  this  to  be  so.  Indeed 
it  will  be  observed  thot  each  curve  in  Figure  29  is  symmetrical  obout  X  =  mB,  which 
corresponds  to  X  =  16  for  the  fourth  harmonic. 

Figure  29  she  wed  the  sound  pressure  levels  calculated  for  on  elevation  of  10  degrees 
below  the  re  ar  disc  plane.  Figure  31  illustrates  the  variation  of  the  sound  level  at 
other  elevations  and  shows  how,  as  the  point  of  observation  is  moved  away  from  the 
disc  plane  toward  the  rotor  axis,  the  number  of  loading  harmonics  contributing  to  a 
specific  ace.  stic  harmonic  is  reduced.  Indeed,  in  the  limit,  immediately  under  the 
rotor  disc,  only  the  single  loading  harmonic,  X  =  mB  contributes.  The  effect  may  again 
be  understood  by  reference  to  Equation  (36) .  The  argument  of  the  Bessel  function  terms 
includes  the  'actor  y/r.  Thus,  moving  away  from  the  rotor  disc  reduces  the  magnitude 
of  the  argume  nt  and,  hence,  the  range  of  effectiveness  of  the  loading  harmonics.  In 
the  . limit  when  y/r  ~  0,  then  only  the  JQ  term  has  a  finite  value,  equal  to  unity.  The 
increase  in  amplitude  of  the  loading  harmonics  on  moving  toward  the  rotor  axis  may  also 
be  Understood  by  reference  to  Equation  (36).  The  thrust  term,  which  dominates  the 
results,  is  multiplied  by  x/r.  Moving  toward  the  rotor  axis  gives  increasing  values  of 
x/r  and  accounts  for  about  a  15-dB  increase  between  the  10-degree  and  80-degree 
cases.  The  i  smaining  increment  comes  from  the  increasing  peak  magnitude  of  the  lower 
order  Bessel  functions.  The  fact  that  sound  radiation  immediately  under  the  helicopter 
tetor  is  dependent  on  a  very  limited  range  of  loading  harmonics  may  be  of  some  poten¬ 
tial  significance  in  naval  applications.  However,  it  should  also  be  noted  that,  imme¬ 
diately  under  the  helicopter,  the  refractive  effects  of  the  downwash  can  also  be 
important. 

The  key  requirement  now  is  to  be  able  to  predict  which  loading  harmonics  are  necessary 
for  the  acousric  calculation  of  any  giver,  sound  harmonic.  From  considerations  of  the 
basic  feature':  of  Bessel  functions,  it  is  possible  to  show^  that  the  range  of  interest  of 
loading  harmonics  is  roughly 

m  B  (1  -  M)  <  X  <  m  B  (?  +M)  (67) 

where  M  is  the  rotational  Mach  number.  The  accuracy  of  this  equation  may  be  checked 
in  Figure  29  where  M  =0.5.  The  equation  will  be  found  to  be  conservative  when 
applied  to  Figures  30  and  31  for  points  removed  from  the  rotor  disc.  In  fact,  here  the 
formula  mB(l  -  My/r)  <  X  <  mB  (1  +My/r)  can  be  used  if  desired.  However,  for  general- 
purpose  calcs,  lotions.  Equation  (67)  is  appropriate .  Figure  3?  shows  the  effect  of  rota¬ 
tional  Mach  number  on  the  noise  radiation.  It  will  be  observed  that  the  range  of 
loading  harmonics  which  contributes  to  the  noise  is  substantially  increased  as  Mach 
number  Increases.  Again  the  effects  are  as  given  by  Equation  (67).  Figure  33  gives  a 
plot  of  Equation  (67)  and  may  be  used  to  determine  the  range  of  loading  harmonics 
necessary  for  accurate  calculation  of  any  given  noise  harmonic. 
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Figure  34  shows  the  contribution  of  the  three  force  components  to  the  overall  sound 
pressur The  proportions  of  components  used  in  this  and  all  the  other  calculations  of 
the  pa  j-meter  study  are  thrust:drag:outward  component  -  10:1:1  .  It  will  be  observed 
that  in  the  inefficient  region  of  radiation  (well  down  the  shoulders  of  Figure  29),  the 
effects  of  all  three  components  are  about  equal .  However,  over  the  efficient  radiation 
range,  the  thrust  term  dominates.  The  substantial  movement  of  the  plots  oetween  suc¬ 
cessive  harmonic  levels' corresponds  to  moving  past  lobes  in  the  sound  patterns,  as  will 
be  discussed  in  more  detail  later.  Note  how  the  thrust  and  drag  terms  show  similar 
fluctuc.  Hons,  but  in  "antiohase"  to  the  outward  component  term.  Thus,  t"<»  curve  for 
overall  level  which  combines  the  three  components  shows  smaller  fluctuations.  The 
dreg  term  is  very  inefficient  in  the  central  region.  This  is  the  effect  of  the  (n-X) 
multiplier  in  Equation  (36) .  For  most  practical  purposes,  it  would  appear  to  be  a  good 
approximation  to  consider  the  thrust  term  only.  The  other  terms  contribute  significantly 
only  directly  in  the  plane  of  the  rotor  disc  and  in  the  inefficient  radiation  region;  for 
instanc  e,  the  steady  loading  contributions  at  low  rotational  Mach  numbers. 

Combi  rati  on  of  Loading  Harmonics 

All  the  data  presented  so  far  apply  to  the  effects  of  a  single  loading  harmonic  .  In 
reality,  the  observed  noise  at  any  point  results  from  the  sum  of  the  effects  of  all  har¬ 
monics  As  was  discussed  in  Section  5.2,  the  sums  of  the  squares  of  the  contributions 
of  the  individual  loading  harmonics  must  be  taken  for  a  realistic  result.  The  key  effect 
in  this  summation  is  the  comparative  magnitude  of  the  loading  harmonics.  This  problem 
was  discussed  in  Section  4,  and  it  was  shown  how  a  harmonic  inverse  power  law  approxi¬ 
mation  seemed  to  be  the  most  accurate.  In  order  to  study  the  effects  of  various  har¬ 
monic  .oading  laws.  Figure  35  has  been  prepared.  The  figure  corresponds  to  a  summa¬ 
tion  of  the  results  like  those  presented  in  Figure  29  with  an  appiopriate  weighting 
according  to  the  harmonic  loading  inverse  power  law  applied.  The  first  60  loading 
harmonics  arc  summed  for  each  point.  The  results  apply  10  degrees  below  the  disc  of  a 
four-blade  roto: . 


Figure  35a  gives  the  results  of  the  zeroth  pov/c-  law  case,  which  corresponds  to  the 
direct  summation  of  all  the  harmonic  levels  with  equal  weighting.  Note  how  the  noise 
is  predicted  to  go  up  with  harmonic  order,  actually  rising  at  6  dB  per  doubling  of  order. 
This  moy  be  compared  with  Loewy  end  Sutton's  rf  suits,  where  their  input  of  whet 
was  es-  enticlly  a  zeroth  harmonic  power  loading  luw  gave  rapidly  decreasing  sound  har¬ 
monics  .  Although  this  must  be  partially  due  to  the  limited  range  of  input  harmonics 
(0-20) .  it  is  thought  that  their  results  also  strongly  suggest  an  undetected  computational 
error.  Note  that,  because  of  the  limited  number  (60)  of  loading  harmonics, used,  the 
oresen  calculation  loses  accuracy  cbove  a  sufficiently  high  sound  harmonic.  Referring 
to  Fig.  re  33,  it  may  be  predicted  that  accuracy  will  be  lost  beyond  about  the  eighth, 
ninth  <;.ncl  tenth  harmonics  for  M  =  1  .0,  0.75,  and  0.5,  respectively.  Figure  35a  shov's 
how  this  is  indeed  true.  The  effect  of  Mach  number  is  also  shown  l,i  Figure  35a.  It  will 
be  observed  that  the  eKect  of  Mach  number  is  small  for  this  zeroth  power  case.  This  is 
consistent  with  the  effects  noted  ! n  Figure  32,  where  increase  in  Mach  number  gavt  an 
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toc'eosed  tons*  °f  loading  -horriwiici,  but  not  a  significantly  increased  Amplitude  of 

by  c  he -monte  ~ 


.'v:Fi.5utei.i,35e.aJ.so;giif«.?h4|  results  of  assuming  only  a  steady  loading.  This  corresponds 
to-jrfiij^lo^c  propeller  noise  qialculfttifioj  J  Note  hew-'  the  increase  in  Mach  number 
f»os  o  very,  prprspuncoj  u>>  the  ncisE  rcd'cri"  by  the  steady  feeding.  Tn«  sc  me 

effect  may  be  observed  by  studying  the  intercepts  of  the  various  curves  in  Figure  32 
yrith  the  vertlsai  ox’s  -  in  feset" Figure  35o  also  shows  how  the  steady  loading  is  por- 
hcu’twJy  high  fix  the  ,ta  s.1 ,0  co$e .  ...This  effect  may  be  understood  by  reference  to 
Equation*  *  (37) .  Both  a  ^  qnac  Jn+X  +erm  are  present.  In  .general,  as  was  men- 

tioned  above,  the.J^^  term  ismucb  !arger  than  the  j  .erm.  Howe/er,  for  the 

X=?0  COM/  the  “terms  are  identically  equal  and  are  thus  additive .  Tht  lower  first  hcr- 
mpqic  Joqding  points  In  Figure  '32  result  from  the  particular  phasing  used  in  the  cal¬ 
culations  and  are  net  c  general  result. 


Figures  35b  to  ■  35d  give  the  results  of  assuming  ether  loading  harmonic  power 
Tdsys.T  In'  each  cose  the  hooding  ho-monics  have  been  derived  from  the  zeroth  by  using 
the  appropriate  inverse  power  law,  as  shown,  (The  first  and  zeroth  loading  harmonics 
.fl.heqyjTOrve equal  magnitude .)- In  each  figure,  the  result  for  steady  loading  only  is  oiso 
rshown,'  This  djey  be  considered  os  the  result  for  the  minus  infinity  power  Saw  for  looding 
harmonics.  The  inverse  first  power  law  case  is  shown  in  Figure  35b,  For  M  *  0.5 
and  0.75,  th*  result  is  very  closely  a  constant  sound  level  for  all  harmonics.  This  would 
correspond  to  a  very  impulsive  banging  type  of  nqige.  For  M=  1.0,  the  sound  rises  with 
frequency .  How  ever  it  will  be  observed  tbet  the  effect  of  the  steody  iood'ng  alone 
gives  virtually  all  iht  'bserved  effects.  Thus,  it  appears  That  the  observed  levels  in 
Figure  35b  at  M>.  1 ,0  are  basically  due  to  steody  iooding  effects,  with  fluctuating 
forces  ploying  a  minor  rule.  Figure  35c  shows  the  inverse  square  law  loading  results. 
Here  'the  curves  drop  off  more  rapidly  with  increase  in  sound  harmonic.  It  rnOy  be 
observed  that  change  of  rotational  Mach  number  is  starting  to  hove  a  more  significant 
effect.  Also,  it  can  be  seer,  that  the  M  =  0.75  case  Is  governed  by  the  steody  looding 
for  the  first  five  harmonics,  and  only  above  this  is  there  any  significant  effect  of  the 
fluctuating  loads.  It  may  be  noted  that,  because  of  the  markedly  reduced  effect  of  the 
higher  loading  harmonics  in  this  inverse  square  loading  law,  the  inaccuracies  introduced 
by  the  limitation  to 60  loading  harmonics  disappear,  and  Figure  35b  is  probably 
accurate  out  to  the  sixteenth  sound  harmonic.  The  inverse  third  power  iow  resuits  given 
in  Figure  35d  show  ell  .htse  effects  to  an  even  ryealer  degree,  with  the  steady 
loading  dominating  the  M  =  0.75  case  our  ro  adout  the  twelfth  harmonic  end  the 
M  =  0.5  case  out  to  the  fourth . 


All  these  results  applied  to  the  ccse  10  degrees  below  the  rotor  disc,  which  is  one  of 
the  prir.».ipcl  directions  of  interest  from  the  detectability  point  ot  view.  The  same 
general  effects  occur  in  oil  directions  u;tdt  the  helicopter,  but  .natural iy  detail  changes 
do  occur.  Figure  3b  shows  the  effects  10  degrees  from  the  rotor  axis,  in  this  position, 
the  erfect  of  the  steody  loading  is  extremely  small  at  aii  Mach  numbers,  indeed,  the 
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second  and  higher  horn  ionics  are  too  small  .to  be  plotted  on.  hbe  ‘•Tb.d*#i 

Ul'r  .IOa  4.-a  *Ua  '*'X;  j  j  1 ’'Ijp  'iCfin  f  tfl£  •  T'’ 
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spend  closely  to  9  straight  line  an  *ne  iog^iog  plot  .  FigvfP  36;?>ay-bs  regarded  to  some  r 
extent  os  giving  the  basic  effects  of  The  fluctuating -forces, -uncamoiica  ted -ovutbeeffec  tv 
of  the  steady  feeding  which  had  Ts  be  jneldaBS-ia  Figure  3o.  -  .  ._•  •_ 
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retical  argument.  Reference  to  figure  29  shows  that,  ovec  the  .range  :of  jwjximymjjpisp^- 
production,  higher  sound  harmonics  produce  mote  noise  thcn:tqw«r<r>P? ,  regce :rr 

18  it  was  shaun  how  the  acoustic  efficiency  was  gctuaily-proporfianahto  m8M  .  The  v ’• 
peak  noise  levels  for  the  various  ham.oriic»  in -Figures  29  through  32  will  belfoqpd  to  -•;• 
follow  this  lew  closely.  However,  there  is  no  basic  trend  due  to  loading  hgrmoni.c  on 
any  given  frequency  observable  in  Figure  29,  apart  from  the  limit  in  the  Oufflber  tijohr'.-.. 
contribute.  As  suggested  by  fquotion  (67) ,  the  loading  harmonics  contribute  over  the  -  > 
range  mB  {1  - M)  <  X  <  mB  (1  +  M) .  furthermore,  sound  intensity  is  proportional  to  ■••-' 
loading  squared,  and,  thus,  if  the  loading  harmonics  follow  arvinverse  k*b  pewer  low, 
the  arguments  above  suggest  that  the  acoustic  intensity  follows  the  taw  given  by 


That  is. 


p2  ~  mBM  j 

-W-M) 


p2  —  (mB)I-J* 


mB(1+M)  - 

X  aX 


(’  +M'l'"zk  -  (1  -M)J~2k  } 


Thus,  the  predicted  sound  power  Saw  variation  with  hannpnic  number  is  * 2 ,  0,  -2,  -4 
for  loodmg  power  laws  of  0,  i,  2,  3,  respectively.  These  laws  agree  virtually  exactly 
with  the  results  shown  in  Figure  36.  The  expression  for  dependence  on  Mach  number 
shown  above  is  not  universal  and  essentially  applies  oniy  in  the  plone  of  the  disc . 

The  law  reduces  to  an  M2  law  for  the  zeroth  power  loading,  and  this  can  also  be  shown 
to  hold  near  the  rotor  axis  for  ail  looding  laws.  This  M2  law  can  be  observed  in 
Figuie  36  -  (Note  that  this  is  based  on  constant  thrust.)  The  results  shown  in  Figure 
35  also  show  the  same  trends.  At  a  sufficiently  high  sound  harmonic  number 
the  curves  follow  the  harmonic  law  stated  above  .  At  a  lower  sound  harmonic  dependent 
on  the  loading  law  and  rotational  Mach  number,  the  effects  of  these  laws  are  overcome 
by  the  contribution  of  the  steody  looding.  The  effect  of  rotational  Mach  number  in 
Figure  35  can  also  seam  to  follow  broadly  the  law  given  in  Equation  (68),  beyond 
the  range  where  the  steady  loading  effects  are  significant.  Thus,  the  equation  does 
give  an  approximate  indication  of  the  effects  of  rotational  Mach  number  on  the  higher 
harmonics  of  the  noise. 

However,  the  most  important  use  of  these  laws  given  in  Equation  (68)  is  in  prediction 
of  the  higher  sound  harmonics,  whicl  cannot  be  readily  calculated  on  a  computer . 
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Using  Equation  (68),  predicted  sound  levels  can  be  extiapolated  with  confidence  to 
the  highest  he  monies.  Furthermore,  this  suggests  that  it  will  be  quite  accurate  to 
determine  the  oasic  blade  loading  power  laws  from  the  acoustic  spectrum.  In  either 
case,  the  con  dation  length  discussed  in  Section  5  must  be  included.  Nevertheless, 
from  these  considerations  it  does  appear  feasible  to  predict  the  amplitude  of  the  higher 
harmonic  blae  ;  ioadings  from  a  measurement  of  the  radiated  sound  signal  . 

It  is  of  particuiar  interest  to  compare  the  predictions  made  above  with  experimental 
data  trends,  f  igures  5  through  7  ga \>e  i  .suits  from  Stuckey  and  Goddard's  work,^  and 
were  discussed  in  Section  2.1  .  It  wus  shown  how  at  low  tip  Mach  numbers  the  first 
harmonic  'eve  s  were  lower  than  the  fluctuating  loads,  while  for  high  tip  Mach  numbers 
the  contribution  at  the  lower  sound  harmonics  was  greatly  increased.  Essentially  the 
same  effects  a  e  shown  here.  Figure  35  clearly  shows  the  increased  significance  of  the 
steady  loading  os  tip  Mach  number  increases.  Because  of  this,  the  first  few  harmonics 
can  be  seen  tc  rise  at  a  far  faster  rate  with  Mach  number  than  the  high  harmonics.  As 
mentioned  in  Section  2.1,  theory  suggests  that  the  fi.st  harmonic  levels  rise  according 
to  a  M^®  law,  and  this  was  shown  to  be  broadly  consistent  with  the  data  trends  of 
Figure  6.  Equation  (68)  suggests  an  M2  law  for  ihe  high  harmonic;,  since  Goddard 
and  Stuckey's  data  correspond  to  a  fairly  rough  running  case.  Again  this  M2  law  was 
objerved  in  the  data  of  Figure  7.  Thus,  the  basic  trends  established  by  the  parameter 
study  are  reflected  in  the  experimental  data. 

The  effect  of  blade  number  B  may  also  be  studied  via  Equation  (68) .  As  the  loading 
inverse  power  law  varies  through  k  =  0,  1,  2,  3,  the  suggested  overall  acoustic  power 
treno  is  through  B2  ,  B° ,  B~2  ,  B"4.  Thus,  the  equation  suggests  that  during  blade  slap 
conditions  when  k  ~  0  or  1 ,  the  sound  power  will  go  up  with  an  increase  in  blade 
number,  while  under  more  normal  flight  conditions,  where  k  ~2,  the  sound  power  will 
go  dawn  with  blade  number.  Although  these  arguments  do  not  apply  to  the  first  few  har¬ 
monics  of  the  round,  these  harmonics  are  not  too  important  subjectively.  Thus  it  ap¬ 
pears  that  an  Increase  in  blade  number  will  reduce  the  noise  radiation  under  normal 
flight  conditions.  This  is  consistent  with  the  considerably  lower  "vortex"  noise  levels 
predicted  by  Schlegel  for  a  five-  and  six-blade  rotor,  compared  to  Davidson  and 
Hargest's  results^  for  a  three-blade  rotor.  However,  the  equations  also  suggest  that  if 
a  multiblade  rotor  gets  into  blade  slap  or  rough  running  conditions,  much  more  noise 
will  be  radiated.  This,  in  turn,  suggests  that  multiblade  rotors  will  not  be  advantageous 
on  tandem  he!  copters.  Little  data  on  this  point  are  available. 

Forward  Speed  Effects 

The  effect  of  forward  velocity  was  discussed  in  Section  3.5.  It  is  clearly  of  consider¬ 
able  interest  r->  study  the  effects  of  this  velocity  on  the  radiated  sound,  and  Figures  37 
through  39  gi v?  some  computed  results.  The  calculations  were  made  using  the  basic 
transformation  discussed  in  Section  3.5.  Those  transformations  were  based  or:  the  re¬ 
tarded  time  pc.ition  of  the  helicopter;  that  is,  the  position  of  the  helicopter  when  it 
emitted  the  se  nd .  The  results  shown  in  Figures  37  through  39  are  based  on  the  actual 
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pos.  i on  of  the  helicopter  when  the  sound  is  observer.  This  requires  an  additional 
fra  reformation,  which  essentially  is  that  followed  by  Garrick  and  Watkins^  and 
Loe  .y  and  Sutton.  ^ 

Figures  37  through  39  again  apply  to  the  basic  case  of  a  single  point  loading  with 
rote  ional  Mach  number  M,  for  a  four-blade  helicopter  rotor,  with  the  observer  10 
degrees  below  the  rotor  disc.  For  zero  forward  speed,  the  sound  field  is  completely 
syrn  metrical  around  the  rotor  (at  least  under  the  random  phase  assumptions  made  here.) 
However,  for  the  case  of  forward  speed,  more  sound  is  radiated  forward  than  aft.  Thus 
the  esults  include  both  the  forward  and  aft  radiation  cases.  The  key  factor  in  the 
thec  y  is  (1  -  M  )  where  MQr  is  the  component  of  the  hub  convection  Mach  number  in 
the  direction  of  tf>e  observer.  Thus,  when  the  helicopter  is  flying  toward  or  away  from 
the  observer,  the  maximum  effect  is  observed;  when  the  helicopter  is  flying  at  right 
angles  to  the  observers'  line  of  sight,  the  convection  Mach  number  component  is 
negligible.  In  Figures  37  through  39,  both  forward  and  aft  radiation  cases  are  given, 
together  with  the  zero  velocity  condition  which,  thus,  also  corresponds  to  the  side- 
way;  radiation.  As  before,  the  effect  of  steady  loading  only  is  given  in  each  case  as 
a  re  erence. 

In  e  ich  case,  the  forward  radiation  is  higher  than  the  aft.  For  the  zeroth  ioading  law 
(Figures  37a  through  39a),  little  change  is  noticeable,  but  for  the  steady  loading  alone, 
the  'orward  velocity  has  a  s"b$tantial  effect.  For  the  0. 125  forward  Mach  number  in 
Figure  37,  the  sound  due  ro  the  steady  loading  alone  is  increased  by  around  20  dB  in 
the  orward  direction  compared  to  the  aft.  Plots  of  the  other  loading  law  cases  in 
Figure  37  show  that  in  each  the  sound  is  increased  in  the  forward  and  reduced  in  the 
aft  directions.  Piets  for  other  rotational  and  forward  Mach  number  cases  in  Figures 
37  tnrough  39  show  the  same  general  effects.  The  higher  harmonics  due  to  steady 
loaoing  alone  are  increased  by  over  50  dB  forward  compared  to  aft  in  Figure  38. 

Figures  37  through  39  are  the  results  for  the  integrated  loadings  under  various  power 
law  .  It  is  of  interest  to  study  the  individual  harmonic  effects,  and  a  typical  plot  is 
shown  in  Figure  40.  It  can  be  seen  that  for  the  forward  cases  a  larger  ronge  of 
harmonics  contributes  than  for  the  aft  radiation  cases.  Figure  40  also  shows  how  the 
effc:t  of  velocity  on  the  first  order  (steady)  loading  harmonic  is  considerably  greater 
the;  on,  say  the  sixteenth. 

The  effects  of  forward  velocity  noted  on  both  Figures  37  through  39  and  on  Figure  40 
are  ery  similar  to  the  effects  observed  simply  due  to  change  in  rotational  Mach  number 
(cor  pare  with  Figures  35  and  32) .  This  suggests  that  it  may  be  possible  to  estimate  the 
effe  ;rs  of  forward  velocity  simply  by  choosing  an  effective  rotational  Mach  number  for 
the  orv/ard  speed  case.  That  this  is  so  can  be  seen  from  a  study  of  the  basic  Equation 
(36)  as  modified  by  the  (1  -  Mor)  term  suggested  in  Section  3.5.  The  argument  of  the 
Bess..*  I  function  terms  is  r,My/i(l  -  MQr) .  Now  n  =  mB  i;  the  sound  harmonic  number 

and  is  unaffected  by  forward  speed.  Nea'  the  plane  of  (he  disc,  both  y  and  r  are 
increased,  by  a  roughly  equal  amount,  due  to  the  effects  of  the  retarded  position  of 
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i  he  helicopter .  Their  effects  therefore  cancel .  Thus,  if  an  effective  rotational 
Mach  number  Me  is  token  such  that  Me  =  M/1  -  MQr,  then  the  effects  of  forwa-d  speed 
on  the  Besel  function  term  ore  simulated.  This  Bessel  function  term  is  the  most  impor¬ 
tant  single  term  in  the  results.  It  governs  the  overall  harmonic  efficiencies  which  are 
shown  grcjjhicaliy  in  Figures  29  and  40.  This  effective  Mach  number  approximation 
takes  care  of  this  term.  Further  study  of  the  equations  shows  that  both  the  drag  and  the 
outward  force  component  terms  are  correctly  approximated  by  this  method.  However, 
the  th  st  term  is  in  error  by  a  factor  (1  -M  )  and  will  be  too  high  by  this  amount  in 
the  effect  ve  Mach  number  approximation .  °iJnfortunateiy,  as  has  already  been  shown, 
the  thrust  term  dominates  the  sound  field.  However,the  error  is  quite  small  for  typical 
helicopter  cases,  of  the  order  of  I  dB  near  the  plane  of  the  disc  for  a  forward  Mach 
number  of  0.125,  and  about  2  dB  for  0.25.  It  may  be  concluded  that  this  effective  ro¬ 
tational  Mach  number  approximation  is  a  useful  tool  for  the  study  of  forward  velocity 
effects . 


Figure  41  shows  the  effectiveness  of  the  approximation.  A  forward  Mach  number  case 
of  0.25  on  top  of  an  actual  rotational  Mach  number  of  0.75  is  compared  to  the  case  of 
rotational  Mach  number  of  1,0.  Since  the  field  point  is  near  the  plane  of  the  disc, 

M  0.75/1-0.75  1 .0.  Figure  41  shows  good  agreement  between  the  more  exact 

i  © 

and  the  approximated  solution. 


It  should  be  noted  that  the  M  is  the  component  of  the  convection  Mach  number  M 

or  o 

in  the  direction  of  the  observer.  It  is  thus  positive  in  the  forward  direction, having  a 

maximum  in  the  plane  of  the  disc  equal  to  +M^.  Moving  the  observation  point  around 

the  helicopter  will  reduce  M  ,  and  it  is  zero  at  the  side  of,  or  immediately  beneath, 

or 

the  rotor.  In  the  af<  direction, M'  becomes  negative,  reading  a  minimum  of  -M  in 

or  o 

the  plant  of  the  disc.  Thus,the  presen.  results  for  forward  velocity,  and  the  approximate 
rule,  exp  lain  one  of  the  frequently  observed  effects  in  helicopter  noise  -  that  the  heli¬ 
copter  hep  a  harsh  popping  sound  on  approach,  but  a  very  moderate  thumping  sound  as 
it  recede  .  This  can  be  attributed  directly  to  the  increased  effectiveness  of  the  higher 
harmonic  in  front  of  the  helicopter  compared  to  behind.  It  will  be  observed  that 
'Figure  3E.  predicts  a  IQ— dB  difference  if  an  inverse  square  loading  law  is  assumed.  Un¬ 
fortunately,  no  published  data  for  direct  comparison  with  this  prediction  are  available. 


•  Directior  ality  Fatterns 


The  next  major  point  to  be  discussed  is  directionality.  Clearly, the  direction  in  which 
the  souna  radiates  is  of  key  importance  in  determining  its  significance.  Some  of  the 
effects  have  already  been  noted,  for  instance  in  Figure  31,  but  the  problem  justifies 
a  closer  study.  Note,  first  of  a!!,  that,  under  the  randomizing  approximations  used  in 
this  repost,  no  variation  in  sound  pressure  around  the  azimuth  occurs.  Little  variation 
is  actually  observed  in  practice.  The  small  differences  in  sound  pressure  measured 
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could  eosily  be  due  to  «*perbnefl»ol  w ror  ^  the  effect  'of ;  th*  ra i I  rotor ,  to  monthly  Vv. 
lust  two  possible  effect .‘'Aiimu'ha'.l  ^ri£tionj' wp>.  ld  be' predicted  •«  <?  definite  pitOM  v  - 
fp,  ^  various  loading  njdnnpn.'cs  Were.' Introduced ;-i nto<  t  Ke calcylot ifin , ,o<r*d  t<*uot;or» . .  _ 
(371  couldfee  w«dtw».  dyf-uno J^heraryaom^nbi^q^ 

.  _  *■*  •  .  -  ,  *i'i  ■"i"--:- ■  j  ^ l  -  _ _'j  ,x;-. i -Jk mV tu.  « 

cow,  iu ch  vonot ionv or e -nor  ?w*s 


cummafr » r  . 

/.  v’.f -v.  •*;.;J: , ■•. 


sourd  .f^r^tare-oredlcfed-on'oossinq.obpve  Of  -  - 

below  the  helicopter ’disc .  -figure  42  shows  tbe  directipnoUty  poftet.>  <vc  to  various 

;nHlv iAiol  force  ®n 

observer  a-C'>;*!a1it  tOC4UJ  trocr,^rse.rTuc,  ano^noycuvrrt/^'^^--™;- --j-  -.7,;;:  v  ...:■-  ■_ 

ticopl  Attach  number  effects  enter  only  tbi;oogb  tbiXeqoency  peromet^r  mfcM.^I 

d&  line  on  a,’*.  plqtsco^ospo*?*  to  the. ov<*«pe>ound  po^ 'rbdkJ^.w . ; 

curves  give  the-correctldn 

ing  low  to  give  the  sound  in  Any 

contours  of  equal  ground.oosse .  As' the  helicopter. flies  over  *nt  Crpu/id,:tbe^nd  ':.r 
radiated  fore  or  aft  bps  to  travel  much  further  Wore  str  ik ing  thegraind  than  the  . sown  ,  - 
radiated  immediately  downward ,  Over-this  additional;  distance,  the  sound.  w.i < j- aheh"  ;„' " 

uaie  according  fo’the  inverse"*quore  law  (see  equation  \2)}..  Th«ido,t+.ed  I^rses "dh  .f i^gur«v  , 
42  allow  for  this  spherical  spreading,  -.0«um|agthatth«  rdh*dJl9 

ground.  Thus,  sound  levels  lying  or,  theAanve  dphed' . ines-w  ill  be  oteeryea  a»; the  *«ne  _ : 
level  at  the  ground.  The  .plats  in  figure  .42  are  given,  rn  ‘he.firtt  quadrant  out  are,  v 
entirely  symmetrical,  so  that  sound  . radiated  up  equal*  sound  radiated  d>vn  .n  _  ■  •  • 

magnitude.  .  .  .  .  .  ■  >. .  • -v,  w-.-.w.5.  v-t** ti'.-'x'.'. 


In  each  of  the  plats,  the  same  general  effect*  oa.n  be  obsitypd.  .-f  av  .large,vatues  of 
mB  -  X  and  law  frequencies  (mB.M),  the  sound  radiation  pattern  is  suhstont tolly  outward, 
in  the  plane  of  the  disc.  This  is,  of  course,  particularly  d.iw4ya_ntageous,.but  rt.shoutd 
be  noted  that  the  overall  sound  efficiency  in. these  case*  is. poor,  with  sound  levels  well 
down  the  shoulders  of  the  curves  in  figure  29 .  On  the  other  hand,  these  cases  assen- 
tiolly  correspond  to  the  effect  .of  low  harmonics  of  the  looping,  which  are  generally  <* 
greater  magnitude.  As  frequency  is  increased.  Of  os  rdB  -X  is  decreased,  the.peo.H  of  the 
directionality  pattern  moves  around  toward  the  «c«i.pf  rotation and  thq. patterns  hove 
a  stronger  and  stronger  lobe  structure.  All  the  cases  with  a  lobad  sfrvcture.con-esppfsrf 
to  efficient  acoustic  radiation,  with  levels  across  the  top  of  the  curves  m  figure  ?9.H 
may  be  rioted  that  on  increase  in  rotational  Mach  number  M  corresponds  .to  a-nK»v«pent 
from  left  to  right  in  each  pattern  matrix  given.  Howevef,en  increase  in  harmonic .oeder 
m  or  blade  number  B  corresponds  to  o  movement  both  across  and  down  (at  least  for 
rhe  lower  looding  harmonics  X),  sincemB  occurs  as  a  parameter  an  both  axes. 


Individual  comparison  of  the  plots  in  Figure  42  shows  that  all  the  thrust  terms  g*  U>.*«0 
in  the  plane  of  rotor  disc.  This  is  the  effect  of  the  x/r  term  on  the  thrust  m  Equotran  «$) / 
and  =  0  at  the  rotor  disc .  Note  ofso  that  all  harmonics  of  all  force  components  go  to 
zero  at  the  rotor  axis  except  for  one  case .  That  is  the  thrust  component  for  mB  -X  -  0,  ‘ 
which  has  a  maximum  there.  Thus, only  the  single  loading  mode  X  -  mB  contributes 
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Immediately  under  the  helicopter.  This  was  also  observed  in  Figure  31  .  The  drag  terms 
shown  n  tgure  42  heve  the  opposite  effect.  All  mB  -X  ;r  0  cases  are  virtually  zero. 

This  can  also  be  seen  in  Figure  34.  This  occurs  because  of  the  mB  -  X  multiplication 
factor  ir  '.e  drag  term  of  Equation  (36) .  Note  that  there  is  also  an  mB  +  X  term  which 
does  not  ;  o  to  zero,  but  this  is  so  small  as  to  be  negligible  in  virtually  all  cases.  Also, 
itr  reiativ  magnitude  is  a  function  of  mB  and  X  separately  rather  than  the  single  term 
mB  -  X  .  7  .us,  no  case  corresponding  to  mB  -  X  -  0  has  been  plotted  on  Figure  42  for 
the  drag  component.  The  mB  -  X  multipl ication  factor  has  another  effect.  When  X  is 
greater  th  n  mB,  it  goes  negative.  The  amplitude  of  the  negative  cases  is  the  same  as 
the  equivalent  positive  case,  but  the  phase  has  been  shifted  180  degrees  with  respect 
to  the  thrust.  Hence,  the  drag  terms  add  to  the  thrust  terms  below  the  rotor  for  positive 
mB  -  X  one  sub*ract  for  negative.  Since  the  thrust  is  dominant,  the  effect  is  of  little 
practical  importance,  but  this  does  explain  the  small  asymmetric^  near  the  plane  of  the 
rotor  disc  n  Figure  43,  to  be  discussed  later.  The  same  general  effects  can  be  observed 
on  the  outward  components  of  force  terms.  These  terms  become  significant  slightly 
nearer  the  rotor  axis  thar.  the  thrust  and  drag  terms  for  any  given  harmonic.  Note  that 
all  the  zeros  in  the  thrust  and  drag  terms  are  matched  by  maxima  in  the  outward 
terms,  ana  vice  versa.  This  is  due  to  the  outward  terms  containing  the  differential  of 

the  Bessel  unction  term  used  in  the  thrust  and  drag  expressions. 

>  1 

Figure  43  uves  the  result  of  the  summation  of  the  three  terms  for  each  loading  harmonic 
oh  the  assumption  that  the  thrust  to  drag  to  outward  components  are  in  the  ratio  10:1:1, 
consistent  with  the  previous  results.  Figure  43  must  be  plotted  in  two  quadrants  because 
of  the  asymmetry  in  the  final  result.  As  discussed  above,  thrust  and  drag  add  in  the 
downward  and  subtract  in  the  upward  quadrant  for  mB  -  X  >  1 .  For  mB  -  X  <  l,fhe  posi¬ 
tion  is  rev*  sed.  The  figures  should  be  inverted  for  this  case.  The  thrust  dominates  the 
results  for  he  total  case,  as  can  be  seen  by  comparison  of  Figures  42  and  43.  The  drag 
and  outward  component  terms  are  only  significant  in  the  immediate  vicinity  of  the 
plane  of  th  *  disc  .  Figure  34  shows  that  even  10  degrees  from  the  plane  of  the  disc 
the  thrust  was  dominant.  An  additional  effect  is  that  the  overall  plots  do  not  go  to  zero 
between 'lobes  as  in  the  case  of  the  individual  components.  Since  the  outward  compo¬ 
nents  are  c  t  of  phase  they  do  contribute  there.  These  valleys  of  the  lobes  are  shown 
only  appro:  imatelv  in  Figure  43. 

The  lobed  suund  patterns  given  in  Figures  42  and  43  are  basically  idealizations  and  are 
unlikely  to  occur  in  practice  for  two  reasons.  First,  the  random  effects  discussed  before 
will  rarely  jIIow  such  an  ordered  pattern  to  occur.  Second,  and  far  more  important, 
the  observed  sound  level  is,  of  course,  the  sum  of  the  contributions  from  all  loading  har¬ 
monics.  The  effects  of  any  one  harmonic  will  usually  be  lest  in  the  overall  pattern. 
Figure  44  g  ves  a  typical  overall  pattern  for  various  sound  harmonics.  It  is  based  on  an 
assumed  2.5  inverse  power  law  for  the  loading  harmonics,  with  realistic  inputs  which 
will  be  discussed  in  more  detail  later. 
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Fi£Kr*  44 c  o  pcior  plot  of  me  so ur.o  radiation  .fj.eld'.in. hover*  .■  As  expected.  the  ■ 

detailed  lobe  structure? 'of  •fi’gL«-«;42- p«<J '43^do  r>ot  opp#or.  in  Ffj^r®  44  at  any  har-  -'••■- 
monte’.  However,  o  gross  di  rec  Nona  Ifry  pattern  f  »:-evtdenK('‘-The  »Os*vcl  fj eld  is  sywwwt-  • 
ric  on  both  siqes-or  t’rs  ios.be  . 

yoc  j  rs  i?r*  i  ‘-tmym  %/%.  ,u-t  >  ww^y  ■*_«:  '*  *!►  ‘  y>v^  »  ■* :  >  ^-».-  -  j  r  — _v  ——  —  .  r  -  —  -  w  —  -— - 

Oo^  {pi>  e$  lOvyGf^  *Hc  /OtOT  .<is£C~f©r_-*^ C'Ki/Q^VCf ;  monies , «  ^ .: ; 

texi  CC  \  ly  ££*£  *0  ff>e  rnrr eH-.fion 

minimum ,  r.:.«'»oyno‘:fTeto‘ T» 

fairiy  uniform, anfl.  In the  high?'  harmonics  ircnJy  varies  a  few &.&. front 

_  y ;  ,^._  I  >.  t_  fTt  ^  <4  «■»  a  e  A/i  rf  ^  irmry  ~  H  f^e 

domizatton  and  harmonic  i coding  Jaws  ossymed/-  On  areolhelicoprer  ro^^jt  /noit  b®'.'; 
expected  That  porticutor  combinations  of  loading  baxmonics  w  iUiccsor  Cependinfl  ®n  ,tt!*-- 
f tight  condition.  This  will  be. reflected  In  o  much  stronger 

Figures  42  and  43.  'Furthermore,  particular  phasings  of -these  loading  comhiodtipns  \«  -  , 

can  give  azimytM  variations  in  'sound  level .  However  fi^/re 

basic  trends.  *. '"  ’  . --■  -  vr 


Figures  44b  and  44c  give  results  for  forward  flight  cases.  Sound  is  given  relative  'toy*?* 
the  position  of  the  he!  i  copter  when  the  sound  is  heard,  as  before .  Both  £  gyres  show 
how  forward  flight  causes  a  genera!  swelling  forward  of  "the  sound  field . ;  ".Details  of 
tiie  minimum  also  vary  somawnat  with 'forward  speed,  but  triis  :sof  compixotively  little 
significance.  In  both  forward  flight  coses  the  maximum  effect  is  observed  in  the  plane 
of  the  disc,  and  also  for  the  second  ond  third  harmonics  plotted,  for  a  forward  Moch 
number  of  0. 125  (Figure  44b),  the  difference  between  the  forward  ond  aft  radiation 
amounts  to  about  5  dB  for  both  the  first  honramc  and  for  the  higher  harmonics  but 
up  to  10  dfi  for  the  seccr.a  and  third.  At  o  forward  Mach  number  of  0-25  (Figure  44c)  - 
the  corresponding  figures  ore  50  dB  and  20  dB^espectively .  ~For  harmonics  higher 
than  the  tenth  (not  shown  hsre),the  differences  are  slightly  less.  ;  -  r 


Necr-Ficld  Effects 

A  limited  number  of  cases  have  been  investigated  in  order  to  estimate  the  order  of  mag¬ 
nitude  of  near-field  effects.  Figure  45  shows  the  polar  distribution  of  the  fiat  harmonic 
sound  radiated  by  a  hovering  four-blade  rotor  represented  by  single  lift  and  drag  forces 
acting  at  the  90-percent -radius  point.  The  radiation  is  in  a  plan#  normal  to  the  rotpr 
disc  .  The  actual  levels  have  been  normalized  to  an  arbitrary  reference  since. the  two 
cases  were  computed  at  distances  of  1  diameter  ond  25  diometers  from  the  rotor. 
Spherical  spreading  corrections  have  been  applied  in  order  to  illustrate  the  direct 
effects  of  the  near-field  terms.  In  fact, the  far-field  results  {25  diameters)  Ogreed.with 
those  calculated  using  Gutin's  equation  to  within  a  fraction  of  a  dB  at  all  potrvt*,whJch 
verifies  the  accuracy  of  the  program.  At  1  diameter  it  can  be  seen  that  the  biggest  ‘ 
influence  of  the  near-field  effects  occwn  above  the  rotor,  although  the  absolute -levels 
ore  still  slightly  less  than  those  which  occur  below  the  rotor.  It  seems  that  in  addition 
to  smoothing  out  the  lobes  the  neor-fleld  pressures  are  more  evenly  distributed  about 


fhe  rotor  plane  than  is  the  radiated  sound.  In  this  case, the  maximum  amplification  in 
the  near  lefd  occurs  10  degrees  above  the  disc  and  is  equal  to  27  dB.  This  corresponds 
to  the  petition  where  the  thrust  and  drag  components  cancel  in  the  radiation  field. 
However,  below  the  disc  the  maximum  difference  is  approximately  4  dB.  This  smooth¬ 
ing  out  or  the  dios  in  the  radiated  sound  by  near-field  effects  was  also  shown  by  Loewy 
and  Sutton. ^  However,  it  should  be  noted  that  the  far-field  lobe  structure  will 
theoretically  occur  sufficiently  far  from  the  rotor.  In  practice, the  lobes  may  not  be 
observed  because  of  atmospheric  erfects,  but  any  observed  smoothing  in  the  far-field 
will  not  be  due  to  near -field  effects. 

Possibly  of  more  interest  than  these  particular  results  for  a  distance  of  1  diameter  are 
those  presented  in  Figure  46,  which  illustrate  the  decay  of  the  near-field  effects  with 
increasing  distance  for  the  first  three  harmonics.  Three  curves  ore  shown  in  each  case. 
The  first  s  the  radiation  field, which  can  be  seen  to  decay  at  6  dB  per  distance 
doubling.  The  second  includes  geometric  near-field  effects;  I .e .,  the  effects  caused 
by  the  special  distribution  of  the  source  at  distances  of  the  order  of  the  extent  of  the 
source.  The  third  includes  both  geometric  and  acoustic  near-field  effects,  where  the 
latter  arc  given  by  Equation  (52).  It  should  be  noted  that  the  acoustic  near  field 
effects  a.  tually  mitigate  some  of  the  effects  of  source  geometry  for  the  first  two  har¬ 
monics,  reducing  the  amplitude  of  the  first  by  some  3  dB  at  a  distance  of  1  diameter. 
The  most  significant  finding  is  that  the  near-field  is  of  little  importance  at  distances 
greater  than  2  diameters  from  the  rotor.  This  is  in  disagreement  with  the  suggestions 
of  Loewy  and  Sutton,  ^  that  the  near- field  may  b?  significant  at  distances  as  great 
as  100  diameters. 

Blade  Motion  Effects 


The  dire'-  '  effects  of  blade  vibration  on  the  noise  field  were  examined  by  computing 
the  sound  field  of  a  rotor  with  steady  loading  but  with  20  realistic  harmonics  of  flap¬ 
ping.  The  analysis  presented  in  Appendix  I  indicates  that  the  maximum  influence  of 
blade  mo:  ion  will  be  found  in  the  close  vicinity  of  the  rotor,  so  the  sound  was  calcu¬ 
lated  at  e  distance  of  1  diameter  at  10-degiee  intervals  from  directly  above  to 
directly  below  the  rotor.  The  only  noticeable  effect  of  the  blade  excitations  was 
found  within  10  degrees  of  the  rotor  axis  where  the  amplitude  of  the  first  and  second 
sound  her  -ionics  was  increased  by  approximately  2  dB.  However,  since  the  absolute 

level  was  close  to  zero  dB  (re:  0.0002  dyne/cm2),  this  is  of  no  practical  significance. 
Elsewhere ,  the  effect  of  the  flapping  motion  was  negligible,  of  the  order  of  0.  1  dB  in 
all  harmonics  up  to  the  fourth.  Consequently,  it  is  concluded  that  blade  motion  has  no 
acoustic  effect  of  practical  importance,  and  no  further  studies  of  motion  effects  were 
made . 

Accuracy  of  Computation  Methods 

Some  of  the  precautions  necessary  to  ensure  accuracy  were  discussed  in  Section  5.0. 
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In  addition, fhe  accuracy  of  the  computer  programs  has  been  verified  by  comparison 
with  previous  work  and  with  each  other.  The  numerical  method  (HERON  1)  and  the 
analytical  method  (HERON  2)  agreed  with  each  other  and  with  Gutin's  equation 
for  the  case  of  steady  blade  loading.  The  actual  comparisons  are  shown  in  Figure  45. 
The  slight  differences  between  the  computed  figures  and  the  Gutin  curve  can  probably 
be  attributed  to  errors  in  the  latter  which  was  hand  calculated,  using  tables  of  Bessel 
functions.  The  two  programs  agreed  with  each  other  in  all  cases  where  valid  compari¬ 
son  could  be  made,  namely,  in  the  far-field. 

The  problems  of  simulating  a  continuous  aerodynamic  pressure  distribution  by  a  number 
of  point  loads  foi  the  purposes  of  acoustic  calculations  were  discussed  in  Section  5.2. 
In  order  to  examine  the  magnitudes  of  possible  errors,  a  series  of  cases  were  computed 
which  clearly  illustrated  the  importance  of  the  phase  relationships  between  the  various 
loading  harmonics.  Using  realistic  harmonic  load  data,  with  definite  phase  relation¬ 
ships  between  the  harmonics,  the  distributed  loads  were  represented  by  1,  2,  5,  10 
and  20  radial  loading  points.  It  was  found,  as  may  be  expected,  that  th'e  results 
varied  somewhat  and  that  the  differences  increased  with  acoustic  harmonic  number. 

In  general, there  was  little  difference  in  the  results  for  10  or  20  load  points  (less 
than  a  dB  or  so  for  frequencies  up  to  the  10th  sound  harmonic);  a  slight  difference 
between  10  and  5  (typically  up  to  2  dB);  but  very  great  differences  among  5,  2,  and 
1  .  At  the  higher  frequencies, for  example,  differences  of  more  than  20  dB  were 
found. 


In  $  action  5.2  it  was  shown  that  the  large  differences  are  to  be  expected  tor  definite 
load  phasing  patterns,  but  that  these  differences  are  illusory.  It  was  also  suggested 
that  the  assumption  of  randomized  ohasing  of  the  spanwise  loads  removed  the  phase 
sensitivity  of  the  acoustic  field.  Accordingly^  computer  program  was  modified  to 
calculate  the  effective  load  distributions  defined  by  Equation  (55)  and  to  use  the 
space  correlation  concepts  reflected  in  Equations  (57) and  (59)  to  calculate  the 
sound  generated  by  a  variety  of  loading  points.  A  thrust  of  10,000  lbs  was  assumed 
with  a  ratio  thrust: drag: outward  component  of  10:1:1,  as  before.  The  value  chosen 
for  the  correlation  parameter  a  was  3/3, so  that  the  asymptotic  form  for 'the  overall 
effective  loading  (see  Equation  (55))  was  Fg  =  F^/V*,  where  is  the  total 


X-hcrmonic  force  amplitude  acting  on  the  blade.  The  particular  value  of  the  correla¬ 
tion  parameter  chosen  above  was  somewhat  arbitrary, although  an  examination  of  the 
blade  loading  distributions  shown  in  Figures  17  and  18  suggests  a  number  of  this  order. 
A  more  realistic  value  could  only  be  obtained  through  experimental  measurement  of 
correlation  patterns  over  a  blade  in  flight  or  through  a  comprehensive  comparison  of 
theoretical  and  experimental  sound  measurements.  A  second  inverse  power  law  was 
chosen  for  the  harmonic  loading  amplitude,  in  accordance  with  Equations  (50)  and 
(53)  .  The  results  of  this  test  are  summarized  in  Figure  47,  where  the  sound  harmonics 
calculated  for  1,  5,  and  10  loading  points  are  compared.  Although  the  first  harmonic 
levels  for  the  10  loading  point  case  were  not  obtained,it  is  clear  that  the  differences 
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are  small,  the  maximum  discrepancy  being  less  rhan  2  dB,  However,  the  most  striking 
feature  of  these  results  is  that  the  errors  now  occur  at  low  frequencies,  as  opposed  to 
the  high  frequencies  encountered  previously.  At  the  higher  frequencies  all  three  dis¬ 
tributions  lead  to  the  same  results.  All  the  present  calculations  have  used  this  random 
phase  assumption  with  a  single  point  loading.  Figure  47  shows  that  errors  introduced 
by  this  are  small . 

6.2  THE  THEORY  IN  PRACTICE 

The  abc/e  parameter  study  has  given  a  considerable  amount  of  information  on  the  basic 
trends  which  can  be  observed  in  the  results.  It  was  also  shown  that  the  general  trends 
predicted  by  the  theory  were  in  agreement  with  the  trends  observed  in  available  experi¬ 
mental  data.  It  is  thus  worthwhile  to  compare  theory  and  experiment  more  closely  and 
to  attempt  to  develop  a  prediction  technique  for  helicopter  rotor  noise  calculation. 

Comparison  of  Theory  and  Experiment 

As. was  noted  in  Section  2.1,  very  little  reliable  experimental  data  on  helicopter  noise 
are  avc'lable.  Consequently  it  has  not  been  possible  to  perform  a  detailed  comparison 
of  the  theoretical  and  experimental  results.  Results  of  a  comparison  in  a  hover  case 
are  shown  in  Figure  12a,  which  gives  data  (and  also  theory)  from  Schlegel  et  al.^  It 
can  be  seen  that  fairly  good  agreement  is  achieved,  at  least  for  the  first  three  har¬ 
monics.  The  fourth  harmonic  is  a  little  low.  The  theoretical  results  shown  in  Figure 
12a  are  the  result  of  the  60  loading  harmonics,  randomized  inverse  2.5  law  discussed 
in  Section  6.1 .  It  will  be  observed  that  these  levels  are  significantly  higher  than 
Schlegel's  theory.  This  is  due  simply  to  the  inclusion  of  a  higher  number  of  loading 
harmonics  in  the  present  calculations.  Figure  33  shows  that  for  the  present  case 
(mB  =  16,  M  =  0.5)  to  up  to  24  loading  harmonics  are  required  for  accurate  calcula¬ 
tions.  Schlegel  et  al  used  Scheiman's  data^  which  are  limited  to  the  first  10 
loading  harmonics. 

A  further  comparison  of  the  present  calculations  and  Schlegel's  results  is  given 
in  Figure  48.  This  applies  to  the  same  S-58  (H-34)  helicopter  in  level  flight  at  a 
velocity  of  40  kts.  This  case  was  chosen  to  correspond  to  a  flight  condition  reported 
by  Scheiman^l  for  which  loading  data  were  available.  For  this  comparisorvScheiman's 
data  were  also  used  in  the  present  program.  For  this  reason, accuracy  cannot  be 
expected  in  the  fourth  harmonic  case.  Agreement  between  theory  and  experiment  ii 
fair  for  the  first  and  second  harmonic,  but  it  falls  off  rapidly  at  the  third  and  fourth  as 
expected.  It  will  be  observed  that  only  fair  agreement  with  Schlegel's  theoretical 
results  is  obtained.  The  maximum  differences  occur  at  maximum  helicopter  range, 
that  is, necr  the  plane  of  the  rotor  disc,  and  reach  above  10  dB  for  the  fourth  harmonic. 
These  differences  are  due  to  slight  differences  in  theory.  The  present  program  included 
both  th*3  radial  force  components  which  arise  through  the  coning  angle  and  forward 
speed  effects.  These  corrections  are  expected  to  cause  differences  near  the  disc 
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P  ane.  It  may  be  noted  that  other  check  calculations  run  during  the  present  study  for 
simple  cases  showed  good  agreement  between  Sehlegel's  theory  and  the  present  one. 

Figure  48  also  shows  the  differences  in  the  sound  radiated  to  port  and  starboard  of  the 
helicopter.  This  difference  is  a  function  of  the  phasing  of  the  blade  toads  only  and 
again  shows  the  sensitivity  of  the  acoustic  field  to  the  loading  phase.  Differences  of 
c  er  20  dB  may  be  observed.  Whether  such  differences  would  be  measured  in  practice 
is  debatable.  Better  definition  of  the  harmonic  loads  would  undoubtedly  remove  most 
o  the  bigger  differences  from  the  results.  , 

It  was  shown  in  Section  5.2  that  ho.monic  blade  loads  may  be  assumed  to  vary  in¬ 
versely  as  the  2.5th  power  of  han  ionic  number.  Using  this  approximation  and  a  single 
leading  point  (as  justified  by  the  findings  repor'ed  in  this  Section  6. 1),  the  rotational 
noise  spectrum  for  ;!.e  B?ll  UH-1  helicopter  has  beer  calculated  for  comparison  with 
available  measurements.^' '  'he  randomized  phase  assumption  was  used  (Section  5.2), 
and  60  loading  harmonics  were  included.  The  comparison  is  shown  in  Figure  49.  The 
experimental  data  were  obtained  on  three  separate  occasions  under  very  different  con¬ 
ditions  and  show  remarkable  consistency.  Because  of  uncertainties  regarding  the  over¬ 
all  levels,  they  have  been  normalized  on  the  basis  of  power  in  the  third  and  higher 
harmonics.  This  step  reduces  the  probability  of  error  due  to  the  !uw-frequency 
response  of  the  microphones  and  tape  r*»/;crcers  which  is  certainly  poor  at  the  funda¬ 
mental  frequency  (around  12  Hz) .  Although  for  this  reason  nothing  can  be  said  about 
overall  levels;  the  agreement,  insofar  as  spectral  shape  is  concerned,  is  good  up  to 
the  thirtieth  harmonic. 


Figures  50  and  51  show  some  further  experimental  data  obtained  by  Wyle  Laboratories 
for  a  twin-rotor  helicopter,  the  CH-47B  Chinook.  No  harmonic  blade  loading  data  are 
available  for  this  aircraft, so  that  no  computations  could  be  made  for  comparison.  How¬ 
ever,  in  both  figures,  the  slope  of  the  fitted  line  is  approximately  -20  dB  per  decade  at 
an  inverse  second  power  law.  When  the  simplified  relationship  between  the  loading 
and  acoustic  power  laws  given  in  Section  6.1  is  used,  namely 


P 


2 


mB 


(2 -2k) 


it  can  be  assumed  that  for  the  Chinook  the  loading  power  exponent  is  1 .5, as  opposed 
tc  2.0  for  the  UH-1  and  H-34  (ignoring  the  correlation  effect  which  adds  a  further 
0.5).  This  rise  in  higher  harmonic  loading  levels  is  almost  certainly  due  to  the  over¬ 
lap  of  the  two  rotors  which  causes  one  to  pass  through  the  wake  of  the  other.  The 
undesirability  of  overlapping  rotors  from  the  standpoint  of  noise  is  thus  fairly  clear. 

Two  further  points  regarding  Figures  50  and  51  should  be  mentioned*  The  first  is  that 
Figure  50  compares  internal  and  external  noise, and, a  I  though  detail  differences  are 
sicjnificant,the  trends  are  obviously  similar  The  second  is  that  one  case  plotted  in 
Figure  51  is  the  analysis  of  a  recording  of  moderate  blade  slap.  The  main  difference 
fct  r  veen  these  points  and  the  second  set  of  data  (which  v/as  recorded  a  few  seconds 
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later  when  no  .blade  slap  was  observed)  is  Ir>  Hie  region  of  the  tenth  harmonic,  u; 
1.0C ili,wiv<ire  the  ievsJs  ere  sese  6  dB  higher  in  the-presence  of  slap. 


AJihoygh.  tfe#  comparisons  ot  theory  and  exp'rimen*  are  somewhat  limited  they  uu  tend 
to  support  the  bdsic -theoretical  results.  More  importanfiy,rhey  show  that  the  lock  of 
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Prediction  Methods  for  ho  rot  Noise 


This  Study,  tike  previous  workon  the  subject,  has  shown  that  methods  for  the  caicuie- 
fion©f  helicopter  rotor  noise  require  extensive  numerical  calculations  which  con  be 
^performed  only  on  a  digital  computer.  Furthermore,  the  accuracy  which  can  be 
.achieved  ,is  limited .. 

-T .^th«u4}^vt-of  the  comparisons  of  theory  and  experiment  described  above,!*  is  con- 
jded  that,  the  most  realistic  approach  of  the  present  time  is  to  make  use  of  the  assump¬ 
tion  pf-rar^d.om  diode  Ipadjng  phase.  This  gives  results  which  are  at  least  os  accurate 
jqs  rnehhi  ds  which  make  use  of  known  phase  relationships,  if  not  more  so.  It  also  offers 
the  advantage  ft»at ,  in  conjunction  with  the  power  law  assumption  for  estimation  of 
bernyonic  loading  levels,  it  can  be  used  to  calculate  high  sound  harmonics.  In  addi- 
■•tiprwit  hq$  been  shown,  that  although  the  forward  velocity  of  the  helicopter  has  an 
Important  effect  oo  t.ie  radiated  rotor  noise,  it  can  be  accounted  for  with  reasonable 
“  accuracy  by  ujing  the  concept  of  an  Heffective  Mach  number  “  . 

Consequently, 0  set  of  tentative  design  cnarts  has  been  prepared  which  enables  the 
joutv,  field  of  a  rotor  to  bft  calculated  for  any  conditions  of  steady  flight.  These 
charts  art  presented  os  part  of  Appendix  ill,  where  the  detailed  instructions  for  their 
use  rijay  also  be  found.  A  review  of  the  underlying  theory  is  included.  With  careful 
use  of  these  charts,  a  few  simple  hand  calculations  will  yield  any  reasonable  number 
of  rotational  noise  harmonics,  at  any  point  in  the  for  field  of  the  rotor,  to  within  2dB 
of  the  accuracy  obtained  by  compu»ei .  The  experimental  comparison  discussed  above 
suggests  that,  although  the  design  charts  may  be  in  error  for  the  overall  levels, they 
should  give  the  parameter  trends  quite  accurately.  The  chorts  should,  therefore,  be 
useful  tools  for  design  trode  -off  studi  es . 


6.3  METHODS  FOR  REDUCING  THE  NOISE  OUTPUT 


The  osults  of  the  last  section  show  fair  agreement  between  the  present  theory  and  avail 
able  data.  It  is  therefore  of  interest  to  suggest  methods  for  reducing  the  sound  output 
of  the  helicopter  rotor.  By  no  means  have  all  the  points  made  in  the  theory  been 
verified,  but  it  is  thought  that  the  basic  trends  should  be  correct. 

The  r.ost  useful  data  from  the  noise  control  viewpoint  were  presented  in  Figures  2  and 
3.  If  was  shown  how  the  sound  radiated  by  the  helicopter  rose  substantially  at  both 
high  and  low  values  of  collective  pitch,  due  to  still  and  wake  interactions,respectively 
It  therefore  appears  worthwhile  to  attempt  to  define  optimum  collective  pitch  settings 
for  minimum  noise.  As  far  as  is  known,  little  work  has  been  accomplished  on  this. 

Note  particularly  that  some  of  the  results  shown  on  Figures  6  and  7  are  as  much  as 
10  dB  lower  than  expected.  It  seems  quite  possible  that  this  is  a  real  effect  due  to 
collective  pitch  variations. 

The  oasic  mechanism  underlying  the  effect  of  collective  pitch  is,  the  displacement 
of  ti  e  shed  vortex  wake  as  far  as  possible  beneath  the  oncoming  blade,  so  that  har¬ 
monic  airloads  are  substantially  reduced.  Any  effect  which  achieves  greater  wake 
displacement,  less  vortex  strength,  or  less  blade  load  response  to  the  fluctuating 
aerodynamic  input  will  reduce  the  noise.  Work  on  increasing  the  size  of  the  vortex 
and  thus  reducing  its  peak  velocity  was  reported  by  Sternfeld  et  al  .,^8  and  this  has 
considerable  potential  for  noise  reduction  purposes. 

It  should  be  noted  that  collective  pitch  is  the  important  parameter  for  optimization 
based  on  the  operation  of  a  particular  rotor.  For  a  given  rotor  producing  a  given 
amount  of  ihrust,  the  downward  velocity  of  the  wakes  is  essentially  constant,so  that 
the  vertical  distance  between  a  blade  and  the  vortex  trailing  from  the  tip  of  the 
previous  blade  is  increased  by  reducing  the  tip  speed.  To  do  this,  of  course,  collec¬ 
tive  pitch  must  be  increased.  Beyond  a  certain  limit  the  blade  stalls, causing  the 
increase  in  noise  which  has  been  noted.  The  tip  speed  thus  affects  noise  in  two  ways: 
through  the  direct  effects  of  Mach  numbers  and  through  the  blade/wake  spacing. 

Another  step  which  would  appear  to  be  beneficial  is  to  increase  the  number  of  blades. 
The  trailing  vortex  from  each  blade  would  be  of  less  strength;  but,  more  importantly,  if 
the  blades  were  sufficiently  close  together,  the  wake  would  not  have  had  time  to  roll 
up.  Results  on  wings'^  suggest  thai  a  distance  of  the  order  of  one-half  the  span  is 
reqe  ired  for  vortex  rollup  to  occur.  It  appears  that  if  each  blade  was  close  enough  to 
the  preceding  one,  rollup  effects  would  be  averted.  Indeed  it  can  be  seen  that  for 
very  closely  spaced  blades,  the  wakes  would  be  affected  more  by  their  neighbors  than 
by  ihemselves,  and  rollup  might  be  substantially  delayed.  Possibly  this  is  a  partial 
exp  ■□nation  of  Schlegel  e*  al's  results,  ^  (Equation  (4)),  which  implied  7.5  dB  less 
"vortex"  noise  radiation  from  five-  and  six-blade  rotors  than  Davidson  and  Hargest's 
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j  r«ulri\tfe"qj,t.Hrc?-biade  rotor  ^Equation  {2)) .  furthermore,  as  noted  in  Section  2,1,, 
— .'Aj.Q6s.onq  'Tcsutti^  sagged  that  a  three -blade  rotor  produces  3  aB  less  noise  then 

.. i l-^a.'.t^ro.hojigKie. prje-^iTIve-tr ends  for  Wade  .number  suggested  in  the  ineoiy  in  Section  6.  5 
— - rArc-typicoijY.o  .B^-jaw,  and  this  is  certainly  consistent  with  the,  above  results.  The 
;->j^t^^Stn;fr4«»u«ncY^wh.ich'W';Juid',coaempgr>«  the  use  oLaiv;-*  fcladesjs-X'ld  have  a 
— -double  effect^  At  short  distances  from  the  rotor  the  effect  would  be  undesiraote  In 
— —ttfai.tne_hqrmCt'r-^i  y.-Qu la ibe Shifted  to  frequencies  where  the  hearing  is  more  sensitive. 

.. -T^rthtf ^Othaf^Ond-.  atmospheric  att-Payction  of  sound  irxteases  with  frequency,  so  that  . 
thAfregy.eocy  increase  cOuld  result  in  a  further  decrease  in  detectability. 

V .PisV-gfav *«$•*.(/  vrjdesi.aWa  feature  from  the  wake  interaction  ooint  of  view  is  the  use 
--■of  mdlt|pl.e_-rotoru;heiicopters ,  Here,  one  rotor  con  interact  with  the  woke  of  the  other, 

‘  orvd  figure  51  (.discussed  on  poge  72)  shows  how  the  higher  harmonics  are  increased  in 
,  rg  .tandem  rotor  case.  From  a.  general  viewpoint,  it  might  be  possible  to  design  a  quieter 
-  ’  tandemkotar  Helicopter  =if  the  vortex  interactions,  which  are  the  major  noise  producers, 
:"v day td.be  minimized,  while  at  the  same  time  allowing  sufficient  rotor  separation  to 
,prp<mitTthe.wqke  to  diffuse.  However,  this  seems  to  pose  a  difficult  design  problem.  If 
■  may  alia  be  noted  that  interaptions  of  the  main  rotor  with  the  fuselage  and  the  toil 
•’rbhy  con  also  be  undesirable  acowsticqily  because  of  the  local  blade  loading  incre- 
merSts/thjty produce.  Indeed,  it  should  be  possible  also  to  reduce  the  noise  radiated 
by  the  fail  ratpr  by  minimizing  intercchon  effects  due  to  the  main  rotor.  It  is  fairiy 
__  straightforward  to  minimize  these  Interactions  by  design  if  desired. 

An  alternative  way  of  minimizing  wake  interactions  is  to  use  high-lift  airfoil  sections 
on  the  blades  Such  rotors  will  tend  to  direct  the  woke  well  away  from  S’icceeding 
glades.,  thus 'rtuucing  higher  harmonic  airloads  and  sound  suppor. .  This  will  generally 
-.^.-result  in  higher  disc  loadings  ar.d  increased  frequencies,  so  that  sour. H  </utput  at  the 
r  -.fimdbmentgt  frequency  becomes  more  important.  Reduction,  of  tip  speed  when  using 
.•/  such  devices  is  therefore- indicated . 

rThe  theory  also  indicates  that  noise  output  is  proportional  to  the  product  of  thrust  and 
disc  loading.  Since  the  thrust  is  not  a  design  variable,  the  rotor  diameters  should  be 
as  great  as  possible  for  minimum  noise. 

It  hos  also  been  shown  that  the  sound  radiation  theoretically  rises  as  the  square  of  the 
tip  speed.  Unfortunately,  even  a  reduction  from  7CC  to  500  fps  in  tip  speed  will  give 
only  a  reduction  of  about  3  dB  according  to  this  trend.  On  the  other  hand,  it  has 
been  demonstrated  that  the  sound  output  at  the  lowest  frequency  is  proportioned  to  a 
high  power  of  the  tip  velocity.  Although  this  First  harmonic  sound  output  is  only  pf 
limited  significance  in  the  far-fieldh  can  be  more  importan1  in  the  internal  noise 
field,  because  of  potential  increases  in  subjective  response  at  high  levels  due  tc 
various  forms  of  coupling  through  structural  vibrations.  Particularly  for  high  blade 
numbers  it  is  probably  well  worthwhile  to  adopt  a  minimal  tip  speed  for  noise  control . 
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A  further  possible  method  for  reducing  the  noise  i-i  to  tncreose  the  blade  chore  .  Re¬ 
sults  derived  in  Reference  1 8  showed  that,  for  a  fixed -frequency  input,  the  integrated: 
aerodynamic  toad  was  proper flcaajjp  the  inverse  square  root  of  the  chord  length; 

Thus,  sound  output  would  be  predicted  a s  inversely  prepen  iunal  «- 
increase  of  chord  at  fixed  thrust  and  coiiecrive  pitch  --o  Uici  cs  c*ccr«e«i*c  *n  r“r 

with  its  associated  acoustic  benefits.  Perhaps  any  increase  m  rotor  blade  area  'would 
be  more  ett>ct«wiy  employed  us  un  extra  v’pdc .  Me  vert  I  ess.-  It  is  thought  that  a  ~ 
mare  detailed  study  of  blade  chord  effects  would  be  valuable,  — 


One  last  feature  of  potent iai  use  is  the  'airectionasity  crswo^ferisi ics .  A  suite,".'  -  ~~ 
definite  minimum  just  above  the  plane  of  the  disc  is  predicted.  It  may  be  possible  to. 
aesign  or  fly  c  helicopter  so  that  this  minimum  occurs  at  the  field  position  where 
minimum  noise  is  desired.  It  should  be  particularly  noted  however  that  this.  minimum 
has  not  been  confirmed  exper5-,entally . 


It  v,-ould  seem  well  worthwhile  to  perform  detailed  (possibly  scale  model)  experiments  . 
fo  study  some  of  these  noise  control  methods.  The  experiments  should  be  deliberately 
designed  to  cover  cases  outside  normal  operating  ranges  so  that  the  trends  can  be  well 
defined.  Such  experiments  could  be  of  considerable  value  in  reducing  helicopter  rotor 
noise  radiation . 


The  major  design  requirements  for  minimum  noise  can  be  summarized  as  follows: 

•  Low  tip  speed 

•  Large  number  of  blades 

•  Low  disc  loading 

•  Large  blade  chord 

e  Minimum  interference  with  rotor  flow 

•  Any  features  which  will  reduce  rhe  high  frequency  airload  fluctuations. 


■v  : 

£ 

*  ^ 
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7.Q  CONCLUSIONS 

'"The -prew nt  ■  verb  has  uncoverea  many  new  ♦actors  of  :;£r»pa>iu«ce  in  the  study  of  heti- 
.epprer  rotor  noise .  in  order  to  categorize  tbs  resvit;,  the  co.no iujioni  wifi  be  grouped 
vn&sr  -severe!  subheadings . 

.  'It.:./  ~  BASIC  McCriANiv**^ 

g)  farrfield  noise  radiation  from  tb*  helicopter  h  dofftir.sfec  by  main  roloj 
n  -  .'• — —  '-rufite  . 

b),  jMpip  rotor  noise,  especially  at  the  higher  frequencies,  is  due  tafhe  high 

■  ’  --.v ;•  harmonics  of  the  fluctuating  airloads  Imposed  upon  it, 

■  “  *  jc)  Blade  slap  should  not  be  regarded  as  essentially  different  from  other  forms 

of  rotor  .noae . 

4)  Discrete  frequency.peaks  dye  tq  main  rotor  rotational  noise  are  observable 
:  to  400  Hz , 

-h-  -  * 

p)-  ‘H*Vprt«xu  noise  as  previously  defined  includes  o  substantial  contribution 
r.v  -u  -V  -drom  rotational  noise. 


|  '  _.IL,  EXPERIMENTAL  DATA  REVIEW 

^ •  v  : « — - — — : - 

{  Ii\  V  -  V-  -  ^  ^ 

.  ;  .  ...  a)  The  most  striking  feature  of  available  data  is  its  wide  scatter,  perhaps 

■J  , ,r  •;  t  lo ds . 

'  b)  Very  little  experimental  data  on  helicopter  noise  ore  yet  published, 

J  '  *  ...  p)  High  harmonics  of  rotationol  noise  obey  a  tip  velocity  squared  law  for 

f  -  consent  thrust. 

v  .v  _  d}  Low  bormanics  of  rotational  noise  obey  a  tip  velocity  to  the  sixth  or  eighth 

~-ii  pewer  lew  (constant  thrust)  . 

>’  e)  Spectral  content  of  helicopter  noise  is  strongly  dependent  on  tip  speed. 

f  f)  Noise  output  in  all  harmonics  obeys  a  thrust  power  law  between  1  and  2. 

”  _  g)  Noise  levels  increase  as  the  rotor  approaches  stalled  operation  (high  col- 

.  lectivt  pitch) . 

! 

.  „  h)  Noise  levels  increase  at  law  collective  pitch  due  to  wake  interaction 

\  effects . 

i)  An  optimum  collective  pitch  for  minimum  noise  exists. 
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Ill .  ACOUSTIC  THEORY 


Excellent  agreement  has  been  obtained  between  analytical  and rui/neocd.l 
sci'/bonj  to  the  fu*<J-3me«tc}  acoys+>c  equation  »o  regions  whet e  fh*  onaiy  t* 
IvUi  cct> ©*i motions  cue  /olio* 

IV.  AERODYNAMIC  LOADS  STUDIES 

a)  Simple  rotor  theory  is  adequate  for  predicting  performance  aafa  only 
since  predictions  of  secona  and  higher  harmonic  loading  levels  are  badly 
in  error. 

b)  Available  dote  show  a  negligible  effect  of  forward  flight  on  higher  har¬ 
monic  aerodynamic  loads.  Substantia!  levels  occur  even  in  hover, 

c)  Higher  harmonic  loads  are  extremely  localized  on  the  blade  Sport/  bat 
they  are  generally  more  intense  toward  the  blade  rips. 

d)  It  is  impossible  to  predict  the  phase  of  the  higher  harmonic  loads. 

e)  The  airloads  follow  an  approximate  loading  harmonic  power  law,  which 
is  very  approximately  an  inverse  square,  based  on  the  steody  loading  for 
two-  and  four-blade  single  rotors. 

f)  In  rough  running  cases,  higher  levels  of  the  higher  harmonics  occur, 
suggesting  an  inverse  first  power  low. 

V.  COMPUTATION  AND  ACCURACY 

a)  The  present  results  are  substantially  selr-checking  due  to  the  use  of  two 
independent  computer  programs. 

b)  Results  obtained  agree  with  previous  results  obtained  for  steady  looding 
input.  Spcnwise  loading  introduce*  small  errors. 

c)  For  limited  loading  harmonic  inputs,  computation  time  has  been  sub¬ 
stantially  reduced,  to  the  order  of  10  seconds  a  field  point  for  the  general- 
purpose  program  based  on  Equation  (20),  and  to  about  one-tenth  of  a 
second  a  field  point  for  a  special  program  based  on  Equation  (37)  . 

d)  The  phase  of  the  loading  is  equally  as,  oi  move,  important  than  the 
amplitude.  Lack  of  knowledge  of  the  phase  of  the  aerodynamic  loads  con 
be  reflected  in  substantial  underestimation  of  the  sound  produced. 

e)  A  random  phase  assumption  has  been  used  in  the  present  computations,  both 
over  the  span  of  the  blade  and  between  loading  harmonics,  fhis  introduces 
ti  e  necessity  for  definition  cf  a  correlation  length  not  given  by  available 
data,  but  suggested  theoretically  to  be  inversely  proportional  to  loading 
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— ->-b«i'Hs.<»kC  r  ■;<  .am*  Using  ihi*_  assumption  the  sound  field  coe  fee  com- 
Tpyfsri  wflr*  -e'-ipnofcie  accuracy  usmg  a  singie  spanwue  !vc d-V.g  point . 

f) gsjufispJlop  of  $t  point  chardwise  loading  it  probably  conservative . 

Sf  SULTS  Of -PARAMfc  ifcx  bTuDY  . 

V-'-"-11"  v  - - - .  ■— ■• —  .. 

-d}~&  -limi-ted  •rOTT^e:of  loading  'hsrrnaa«cf  contributes  te  any  one  sound 
"harmonic*  -.  - 

ai“  ■  ***-c?b-*»  v--*-u_  -•  “ 

b)  -Theironge  is  centered  an  the  loading  harmonic  X  “  mS,  where  m  Is  the 
sound,  bprrnoaic  and  B  is  the  number  of  blades . 

p)  The  harmonic  X  — m§  tithe  only  harmonic  contributing  to  noise  directly 
on -the  rotor  axis.. 

<Sy  -  in  the  plane  of  the  rotor  disc,  locdiog  harmonics  up  to  nsB  (1  +  M)  con¬ 
i'  tribute  .  •  for  instance,  for  the  tenth  sound  harmonic  of  a  four-biude  rotor 
_T  ,.  with  loading  harmonics  up  to  the  sixtieth  contribute  . 

e)  Neither  .theory  nor  experiment  gives  aata  on  these  high  loading  liarmonics . 

*,■  -  ■  .»  -'Wj  — -V-  ,  ** 

J) .Fflr  typical  helicopters,  the  thrust  fluctuations  dominate  the  noise  field 
except  close  to  the  plane  of. the  disc ._ 

g)  .  For  the  loixiihg  power  law  typical  of  helicopter  rough  running,  an  acoustic 
'  spectrum  identical  with  blade  slap  is  produced, 

h)  For  fhe  loading  power  law  typical  of  normal  operation,  the  acoustic 
spectrum  closely  resembling  normal  helicopter  noise  is  found . 

i)  /  The  steady  loading  dominates  the  levels  of  the  first  few  harmonics . 

j)  Typical  venation  of  the  higher  harmonics  is  as  tip  velocity  squared, 

k)  Typical  variation  of  the  first  harmonic  is  as  tip  velocity  to  the  2B  power, 
where  6  is  the  number  of  blades , 

l)  Typical  variation  of  all  harmonic  levels  is  as  thrust  times  disc  loading. 

m)  Higher  levels  of  sound  are  predicted  in  the  forward  direction  compared  to 
aft  for  forward  flight  cases . 

n)  Forward  flight  coses  ore  well  predicted  applying  on  effective  velocity  rule 
to  the  hover  results. 

a)  The  overall  sound  radiation  pattern  has  a  maximum  at  a  small  angle  below 
the  rotor  disc  and  a  minimum  slightly  cboye  . 

p)  Near-field  effects  are  negligible  more  than  2  diameters  from  the  rotor. 

q)  Blade  motion  does  not  produce  a  significant  sound  field. 


wn 


C  YOCOTAiCNiTA  )  DOCrM/^TT/^M 


a)  Computed  and  measured  ! eveU  of  the  first  fqyr  bonnonics^  ogree^f airly  w«!! 


w  a  hover  cate  . 


(  \  /-  *  i ..  •  •  •  t  .  -  i  i  •  .  ,  _ -j_»_  r .  j _ -r  . ..  v _ i  r:  •  j. 

4.7  uicjnei  <cv-to  uj  c  isyrr-ci  iww.a^Iv  kS?*  -V 


c)  Spectra!  shapes  ore  predicted  we!!  for  q  hover  case  out  to  the  thirtieth 
harmonic ,  -  — •  •  ••  •  . 

a)  Prediction  methods  bosea  or,  the  Theory  should  be  occ-ureie  for  predicting 

♦ ./»  rv  hLr  k.  t  4  rtin  ‘  a  *-»■»-.  w  I  ,-T  n>t  II  I  a>  ;aIv  "*  - 

-  —  *  wv '  -  «f*» »  trg  in  on  *e»  1 1 »  vm&'W  »-*  m  *«,(*■  « — -  _•**■-*“  — 

e)  Design  charts  for  general  use  in  helicopter  rotor  noise  prediction  are 
presented  in  Appendix  III. 


vm.  NOISE  CONTROL 

a)  An  optim..,:  collective  pitch  exists  for  minimum  noise,  which  requires 
detailed  experimental  study. 

b)  Minimum  rotor  noise  output  requires  that  the  rotor  interaction  effects  be 
minimized. 

c)  Configuration  chonges  should  ensure  minimum  interactions  due  to  the 
fuselage  ond  tail  rotor. 

d)  Multiple- rotor  systems  are  unfavorable  acoustically . 

e)  Increase  of  blade  number  seems  ta  offer  the  best  possibility  of  noise 
reduction . 

f)  Use  of  hi gh- lift  blade  sections  should  reduce  high  harmonic  noise. 

g)  Reduction  of  tip  velocity  will  have  significant  effects  only  on  the  funda¬ 
mental,  but  this  mcy  be  important,  particularly  in  conjunction  with  (e) 
and  (f)  above . 
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8.Q,  .  .,v  RECOMMENDATIONS 

Qft.-.tne  baits^pf  theresults  of  xhls  irtv.esii.goti  on,  if  is  recommended t 

:;r 7 J5)- A  systemofi c  experimental  *h>c-y  ha  performed  to  provide  accurate  and 
-•‘-T  ‘detailed  acoustic  data  on  helhropters  of  ptl  types, 


b)  Measurements  s>e  mao*  of  mode;  ond/or  full  scale  rotor  air!  ocas  using  cubic- 
*  frequency  range  irair^rnt-:tatiorw  possibly  la  conjunction  .with  acoustic 
-  —  -  -  - *^,^*^* ^ i  ~?iTTc  t-*c* .*  j. v jjyts  *s  Mn&iy^Cw  inrc  pow’fir  criu  C4,C4S",powt«/ 

-  ~~-*p#ctrp 

l"r  '  ."  s,  The  noise  control  measures  suggested  in  Section  6.3  be  evaluated  through 
'*■'  '  ..  '  a  .model  or  full  scale  study. 

1  V  '  *  4-  -  ' 

~:d)  further- theoretical  work;be  performed  to  specify  helicopter  noise  output  in 
■_. .  ”  terms  of  the  craw-power  spectrum  of  the  loading  inputs  (see  Appendix  Jl) 

✓  -  e)  A  theoretical  study  be  mode  of  the  noise  generated  by  the  motion  of  blade/ 

'  .  vprt ex' i attraction  points  along  the  blade  span, 

/■"  _f)  . The  theoretical  prediction  methods  presented  be  subjected  to  detailed  experi- 
..*•  r-%>  .  jmeata  [“verification,  both,  to  establish  confidence  in  their  utility  and,  possibly, 
-J-'-  to  .suggest  empirical  correction  factors  which  may  be  applied. 
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APPENDIX  I 

ANALYTIC  EXPRESSION  FOP  BLADE  MOTION  EFFECT 


If  would  cieoriy  be  of  interest  fc  attempt  fo  describe  analytically  the  sound  £i?„i_d  wbU 
fit  *  refer  is  flapping.  Let  us  of  tempt  to  do  this  following  the  methods  of  Sec.fi  op  3.4„ 
The  key  change  :5  ;h=  vui-jc  of  j  used  in  me  ,-etcraed  time  integral-;- .*«£•*«>  bled.; 

motion  case,  suppose  the  motion  is  given  by 


x.  - 

i 


yt  !)  sin  9  -  R  cos  S,  -rj  cos©  -  X  jin  8-.. 


taad  q  are  functions  of  6.  £  represents  the  out-of-plane  flapping  motion,and  rj 

••eptesenrs  the  'n-plcne  motion.  Radio!  motions  could  be  included  if  desired. 


7hus,  r 2  -  (x  -  £  }2  m  (y  +  q  cos  6}2  -  2  (y  +  q  sin  9)  R  cos©  +  q2  cos2  6 

t-  2r(  R  cos  9  sin©  +  R2  (6?) 


Following  through  the  same  approximation  as  before,  ignoring  squares  ana  produuis  of 
£  ,  n,  and  R  in  comparison  with  x  and  y,  gives 


r  *  r 

l 


y  R  cos  8 


-  vn 


sin  ©} 


(70) 


Substituting  in  the  expression  for  the  sound  radiation  gives  then  (identifying  y  with  Y 

Ond  8  with  9  -  p) 


w  /  £  ■“ 


4na2  -  .  . , 

O  0  A=  -CD  1 


A)8 


n  r , 


n  Ci  Y  R  . 

- -  cos  (9  -  f) 


0  0  ) 
n  Q  x  £  n  Q  Y  r>  sin  (9  -  $) 


0  t 


a  r, 

C  l 


d© 


(71) 


If  £  ond  rj  are  arbitrary  functions  of  0,  then  evaluation  of  this  integral  is  possible 
only  on  the  computer,  For  some  special  cases  an  analytic  solution  is  possible,  and 
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one  cose  ot"*ame~  interest  i  s  wnen  £  =  c(  cos  to  —  q),  r|  ~C.  This  corresponds  to  u 

ftr.sf  Harmo-ic  flopping  motion ,  The  £  <>nd  Y  terms  in  the  exponential  then  combine . 

then  iimj3t.y'tbg.Kgi'yen  b^’EguotEon  (33),  except  that  the  argument  of  th* 
fte^iifyhcfiorf  ~is 


,a»nce^  (/^  .pqh/be'sofeiy-aispfriect-to  &£.srnoii/  ir  con  be  seen  :ha;  the  sugar  disiorticr, 

erf  the  spuhdlfrhth'pro^cjdi.m^hjs.  cose  h  of  little  practical  consequence  {except  close 
"  to;  the  he  >.•  axis)  .  ,;  • 


UofortvnqrteJy,  use  oj  o  Fourier  series  to- describe  to  flapping  motion  does  not  lead  to  any 
simpler  result  in  thg^presfnt dose,  since  the  various  components  cannot  be  separated. 

It  J*  prpbdb'te  thpt  a  solution  can  be  .obtained  In  terms  oh  Hypergeometric  functions,  but 
thi i  apiproqch'wa^ not  tqken  in  the  present  work.  Thus,  analysis  of  flapping  and  the 
_bkjde^mptidnrfl.^fes:tj  has  be*  n  tiny  ted  to  the  computation  using  the  genera!  program. 


APPENDIX  II 

NOISE  RADIATION  BY  RANDOM  LOADS 


The  basic  objective  of  the  present  report  has  been  to  calculate  the  noise  radiated  by 
the  fluctuating  loads  on  the  rotor.  The  assumption  in  Section  3.4  that  the  loads 
could  be  expanded  as  a  Fourier  series  limited  the  cnalysis  to  the  loadings  which 
occurred  at  some  harmonic  of  the  rotational  speed.  However,  it  is  clear  that  random 
loads  can  also  act  on  the  blades  at  any  frequency, so  that  ?t  is  important  to  be 
able  to  calculate  the  noise  radiation  from  such  sources.  Furthermore,  it  was  shown  in 
Section  5  how  that  even  the  harmonic  airloads  had  to  be  assumed  to  be  random  in  phase 
to  obtain  meaningful  results.  It  is  clear  that  a  more  detailed  analysis  of  the  random 
case  is  necessary.  This  appendix  presents  a  preliminary  study  of  the  problem.  It 
should  be  noted  that  the  analysis  presented  is  not  rigorous  and  several  assumptions 
which  could  possibly  be  of  importance  have  been  made;  but,  although  only  an  out¬ 
line  proof  is  given,  j-he  results  are  of  fairly  clear  physical  significance,  and  can 
be  used  directly  in  calculation  of  noise  radiation  from  experimental  data. 

First  of  all,  since  random  loads  can  occur  at  any  frequency  it  is  necessary  to  obtain 
an  expression  for  the  spectrum  function  of  the  sound  radiation  at  any  frequency .  The 
analysis  follows  that  of  Reference  18,  and  the  results  of  Reference  24  may  also  be 
stud:ed  for  comparison.  Ffowcs  Williams  and  Hawkings  41  have  recently  presented 
arguments  essentially  similar  to  those  given  here. 

We  wish  first  of  all  to  evaluate 


G(u) 


p  exp  i  u  t  d  t 


whe  e  G  is  a  spectrum  function.  The  limits  on  the  time  integral  are  left  undefined 
at  present,  p  is  the  sound  pressure  observed  at  a  point  x  and  time  t,  and  is  given 
by  Equation  (20)  as 


p(x,  t) 


±  |  Fi  ) 
9  t  j  4  n  r  ( 1  -  M^)  J 


(73) 


Equcfion  (72)  appl  ies  fora  single  point  input  on  the  blade .  When  the  variables  are  changed 
bad  to  source  time  t  =t-r/a  ,  Equations  (72)  and  (73)  give 


G(u)  = 


/ 


-  y. 


( 


) 


aor 


3 'r  |4nr(l  -  M  )| 


exp  i  o(t  +  r/a.)  d  t 
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twd  integrating  by  parrs,  ignoring  the  near  field, gives  \see  aiso  Reference  18  > 


^  ^  , ,  V  K 

/_in_ iL-i 


itti  •  “  s/a  a/  8t 


?,f5  rotor  cose,  **e  can  aoply  ike  vSvOl  formula*  for  the  force  components  end 

-retarded  -Oppr.wr  ;mot.lKC.  ' .  .  rv.  yr-rc 

-  G(«l  =  -r  XJ£  sir.  Cl-r  -x£rc«0-r! 

’V  *  '  J  4*Vt(fi  rl  fi  ) 

k.  -t  exp  j i u co*fiT^|dT  (/5) 

■  '  -  {  \  °c  3ori  /) 

-Notation  is  the  same  as  that  used  in  Section  3,4. 

'Nov*  pyt  T  =  T  exp>(-  t  v  t)  and  similarly  for  the  other  force  components,  that  is 

"we. gssume  that  the  forces  are  varying  at  a  circular  frequency  not  necessarily 
eaual  to  a  multiple  of  the  rotational  frequency.  Also  from  formula  44  tr. 
Mclqchlar.'s  booh/* 

'  -  *>■  V 

"*~V;  •  exp  i  2  cos  9  =  /]  inJ^U)  cosn8 

—  -  n  =  -  oo 


Thus,  (75)  becomes 


/{xT  jr-v  /-V  „  ( 

“■>  ; _ 1  +  - -  sin  Q  t - cosQt 

4»Vi  ri  -I  ri  ) 


l__, 


expi  (u-  v)  t 


r 

-c:  _  j*r\  '  f 


cos  n  Qr  dr 


brow  put 

-sin  s2  t 

.11  . 

*  ..  .  ...  j 

*  2  ii  -  *ap  |"  *  0 -T-< 

'  *  "  " mri 

—  srcv  eAP  * 
4i  ‘ 

and 

cos  Q 

JL1  Di 

-  Tj^Pi 

i-  tJ  +  exp  x  ■*  i  1 1 

. 

_ —  ’ .  _ 

on<j  us*  orthogonality  relations  on  the  integration.'  For  the  thjyxt,  tenij;fo* /if|^.n:ce£ 

o/viy  «-Ti£  ‘"0ug3*“CiCj  *  v  7  "j  ^  m 


»«  j^v  ,-  ,r>  ,  yUv  /•  Mrvri  ,  _  c.hi^'TV-'i 

C(U/  =  4^T)~('°  2TT  \^'} 

0  l  \  l  ‘  ■■'T  •■•?* -  •<?'-. 


..  -•-  V„->*£.V 


Ic 

2 


V  /,  ,nr|  , 

—  {(.-;*  J  . 

r  \  .  nrl 


-  {-an*,j 


V) 


n.£  J . 


.  W  .j"  .  -%.' 


where  the  argument  of  alt  the  Besiei  functions  ii  My§/qpr( 


Now26 

Jn-W  " 

and 

Thus,  Equation  (76)  can  be  written  as 


G(u)  = 


■iea. 


ifciv .  jo  M, 7iriv+  ;yCv ,  /£i-vi 

r,  j\ri  u  w  )  n  sce  ri  /  r  fl\°cri/i 


where  G(u)  exists  only  at  u  =  v  ±  rs  CJ 

Thus,  input  of  an  arbitrary  frequency  v  on  the  blade  gives  rise  to  many  frequencies  in 
the  sound  field,  displaced  on  integral  number  of  rotational  frequencies  from  the  input 
frequency.  The  effect  is  the  same  as  is  observed  in  frequency  modulated  radio  signals. 
It  occurs  here  because  of  the  frequency  modulation  due  to  the  varying  Doppler 
frequency  shift  as  the  blode  rotate;  toward  and  oway  from  the  observer.  Equation 
(77)  is  consistent  with  Equation  (36)  derived  for  the  integer  larding  harmonic 


case . 


If  w <&  now  suppose  that  a  complete.  frequency  spectrum  F(v),  the  same  for  each  com¬ 
ponent,  is  present,  then  clearly  ail  loading  frequencies  v  removed  nQ  from  the 
acoustic  frequency  of  interest  that  u  will  contribute.  Thus,  if  the  overall  con¬ 
tribution  at  u  is  S(u),  we  may  write 


+  oo 

S  (to)  =  ^  F  (u  -  n  Q )  G  (u)  (78) 

n  =  -oo 


where  the  T,  D  and  C  terms  in  G  now  contain  only  magnitude  and  direction,  and 
not  frequency. 


It  now  remains  to  extend  the  result  to  the  complete  blade.  The  observed  pressure  from 
fne  total  blade  p^  is  given  by 


PT  <~'f)=y*  Pk/H/Odb 


where  p  is  the  single  point  pressure  previously  determined  and  q 
fixed  in  ti  e  moving  blade.  Thus  also  the  overall  spectrum  function 


is  a  coordinate 
is  given  by 


By  the  usual  rules  for  establishing  power  spectra  from  spectral  functions  (for  example 
Reference  22  ),  the  power  spectral  density  of  the  sound  is  given  by 


P(x,u) 


Lt 

T—  oo 


T  f  f  S*  (x  ,  n,u)  S(x  ,  r,’,(d)  dr,  dn' 

i  J  J  ~  ~  ~  ~  ~  ~ 

n  n' 


(79) 


If  we  write  q'  -  q  +  £  so  that  £  is  a  vector  separation,  then  Equations  (78)  and 
(79)  may  be  combined  to  give 
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■co 


co 


p. 


k-y  -  T-  m  T //  £  £ 

ri  4  *n^*g&  -.  -*-:'‘v  ■  * 


-*j  ftvTY  t  rjjJf  irr'p  irpj'.i  •trt~*»t7  .r;T!j"t’f  ?vv7fr?v»!  *%>*  H- 

vary  iittfe  over  the  typical  dimension  <  ,  then  Equation  {80)  becomes 


■oo 


F(x,u)  -  II E  Pg(n  ,^,W'fr.Q) 

1  4  n  —  -  CO 


0_  {x  ,  n  j 


4  4d  n. 


{8.1} 


Here  I  G|  is  the  absolute  value  of  the  acoustic  transfer  function  given  by  (72]  -and 
can  be  e  function  of  observer  position  x  ,  point  on  the  blade  n  ,  and  input  frequency- 
is  the  cross -power  spectra!  density  on  the  blade  of  frequency  u  -n  Q  and  spacing 

4  -  Pg  has  been  permitted  to  vary  with  rj  »  allowing  far  inhomogeneity .  Pg  con  be 

measured  experimentally  using  a  transducer  array.  In  fact, the  4  integration  con  be 
performed  to  give 


P(> 


. 

>u)  =/  T, 


Au-nQlTl)  V"'"'1"01 


G  (x  ,  n) 

ft 


d  n 


(82) 


where  Pg  is  now  the  spectral  density,  and  A  ^  is  the  correlation  orea  which 
can  be  a  Function  of  frequency  ar.d  blade  position. 

Note  again  that  the  observed  sound  at  any  given  frequency  is  the  result  of  random 
loading  contributions  at  a  wide  range  of  frequencies.-  each  removed  nQ  from  the 

1  1 2 

acoustic  frequency  of  interest.  In, fact  the  transfer  function  |G|  is  virtually 
identical  to  that  used  for  the  harmonic  loadings  in  the  text  of  the  report,  and  it  will 
thus  possess  all  the  features  discussed  in  Section  6,1.  In  particular,  it  will  actasan 
efficient  acoustic  transfer  mechanism  over  the  range  -oM/Q<  n<  uM;  Q;  i.e., 
over  the  range  of  loading  spectrum  spaced  ever  uM  on  either  side  of  u,  Thus,  for  a 
typical  rotational  Mach  number,  M  -0.5,  the  broad-band  contribution  at  any  given 
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jSQOjncy -cOfro.'-frQtti  q|  loading  spectral  den5iti«  wiiht*  a*  octave  bend  Ccnt&re* 

*1-  .-^i:-*—  e  :  '  -  - 


-T  Ka  "  »T\  A  ♦  W/%ri~  n/-  **  u*w  lA^r*.  A  f  ^  il*  .A  ..  W-,  l/v»  —  —4  -»M  ,  .■  «  _  J  *.  _ .  _  i *  .  . 

nr- vv-^; r.. .*? rs. Y’stv^’ *7  W5^-  jitb  r-viivx.  <syvij 

:i?^^^~5^1>^«*wntioiiy'“CO!T«j»ndi  to  -Equctipn  (82)  -Note  that  the  overall  level 
•i>*S9  ?*>,<>■  fr'X  .?h^.;Sap  thjt  SSdF/iS*  8?  th*  dPpinbjjH.ng  ba-rosci-c  loads,  artbihoj  the 

inpur. -specr>un>— ’ ri— ^.cwvt-;pi;^w  J^y  -m  cOff«»wf*OA  vm^c  .  ££lci«  o? 
occyr^fnjfhe  presenf'wcfk » ‘  If  is  of  :nferest  to  rote  that  for  the  harrr-C'nc  input  case, 
.d.  definite  p+^se  relation.  will  ,e^3it  between  the  forpe  input  on  each  blade,  so  that  the 
harmonic*  cancel  put  {ex  frequencies  that  are  not  multiples  of  the  blade  passage  fre- 
quency.vqee  Section  3,5),  However,  no  such  definite  pfiase  relation  will  occur  far 
the  fanddrij  loads,  so  that  oH  vgtuej  of  n  must  be  token  in  (82),  and  the  effect  of 
-6  blpdes  will  st.mpjy  be  to  multiply  alt  output  by  B,  It  should  also  be  notea  that  the 
.. effective  neglect  -of  retarded  tjtnein  Equations  (81)  end  (82)  could  be  &f  signifi- 
'ccnce'in  some  circumstances . 
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DESIGN  CHARTS  r OR  The  ESTtMATJOr*  uj* 

UC!  f  «A  T 1  in/Sk'Vi  «.  i  »  r-r-\ 

*  t  *  *_  l  v —  I  ;  C.A  AV_Ti.»*.UU  'IMV  ' 


Background 


The  equation  for  me  magnitude  of  the  mtn  sound  hannonic  of  c  B-bloded  cat  or  in  ,rh>« 


Kcv 


iV#*f  <r^-r»trA^  O. 


03  / 

Ev  MiM  .,  CXD  nM 

*  .  t- ET"  sin®cXT“t - £—  J-  ->■  -5—  cos  9  c 


A -0 


t  * 


AC  ) 


(S3) 


where  Kuo  constant 


r  is  the  distance  from  ihe  rotor  center  J 

’•  defining  the  field  point  • 

o  is  the  angle  from  the  rotor  disc  planet 

R  is  tne  radius  of  action  of  the  blade  forces 

CAT'  CXD'  an<ScXC  Cf*  ^rustjdrog,  o.nd  radtai  force  harmonic  coefficients 


n  -  mB  =  harmonic  number  x  number  of  blades 

M  is  the  rotationcl  Mach  number  =  QR/o„ 

'  */  ■* 

Jj ,  J2,ond  Jj  are  complex  collections  of  Bessel  function*  of  argument  nMc<*8 

Cn  Is  the  amplitude  of  the  nth 
sound  harmonic  observed  at  the 
field  point  0  when  the  heli¬ 
copter  is  hovering  at  point  x 
(see  sketch)  a  distance  r  away. 

Now  if  the  helicopter  is  moving  along 
the  x-axis  at  a  Mach  number  Mp,  the 

sound  that  reaches  the  observer  at  0  when  the  rotor  is  at  x  wos  actually  generated 
at  some  previous  instant.  At  this  time  the  rotor  was  at  x  so  ‘hot  the  sound  actually 
traveled  a  distance  r* . 


When  the  necessory  transformations  ore  applied  ro  account  for  the  forward  speed  Educa¬ 
tion  (83)  is  written 


•The  rotation  in  this  seefior  differs  from  that  in  the  main  body  of  the  report. 
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*  2-  *•  s 

w,  r  11 < 


1  i'n  3' 


*yV"  5  -  *“;  Ji- 


‘  D  ;  .  nW  cos  9 

r**  «/»  **  r-'^  ’. >  '~l. . is  s 


'2  r  .(!•«;;  xc  i  i 


-whoa  gtc  S^sHuricfisyw  invciveu  nvr-'  .ha its  the  argument  ft  M  cos,  O',  (t  -  Mp,  and 
M*  a  the  component  a/f'Mf  along  the  line  r‘,  but  r!{l  -M,>  a  r  >o  Yhet 


t  -  M 


odd  £«juotIon(84}HmpUf-es  to 


T  ■=£  sm  9  end 


■  1  -  m;  1  -  Mr  . 


E  Kf 

x =0  ! 


1  inM  a  < 

rpr ,m  9  Vr°i  - 


nM  cos  0 

--  u~.\0 


•\C  w3  j 


~Nk?w  »uppo*e  insfeod  of  performing  the  correct  velocity  transformations,  w«  merely 
substitute  -  xyp  -Mr)  into  Equal  ion  (83)  to  give  the  correct  forward  speed  vclue  of 
the  &e**ei  functions-  We  obtain  ^ 


•  <*>  .  .  c 

r  m  V*  v  ^  t  nM  sin  9  .  XD  » 

-  ^  t  CXf  "  R  2 


nM  cos  8  *  I 

*t-*g  -i 


Note  that  Equation  (86)  is  identical  to  Equation {85} with  the  exception  of  the  additional 
( J  w Q)  in  the  denominator  of  the  thrust  terra. 

-  Now  since  Mr  *  cos  8,  th»v(I  *  Mj.)  factor  is  most  significant  os  3  —■  G,  i  >e ,,  os 
.sin  9  —  0.  As  d  increases  and  the  thrust  terns  becomes  important,  Mf  decreases. 
Consequently,  If  may  be  expected  that  the  error  in  rhe  thrust  term  may  not  be  too 
significant.  To  illustrate  this  point,  two  cases  were  computed  .  The  first  was 
computed  using  the  correct  fonvQrd  speed  transformations,  end  the  second  was 
computed  .merely  using  a  modified  Mach  number  {'  !^/)  -  Mp  cos  6).  The  results 
ore  shown  in  Figure  52  for  M  =  0.5,  Mp  —  0.125,  for  the  first,  second,  fourth, 
and  eighth  harmonics  of  c  typical  four-blade  fotor ,  The  random  phase  assumption 
discussed  in  Section  5.2  of  the  report  was  used  for  the  additionally  extreme  cose 
of  random  phase  among  the  thrust,  drog,  and  rciiai  force  components .  (Note 
that  this  leads  to  on  acoustic  symmetry  about  the  rotor  plane,  since  the  phese  dif¬ 
ference  between  the  thrust  end  drag  components  is  eliminated.) 
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If  can  be  seen  in  Figure  52  fhaf  the  maximum  error  in  the  equivalent  Mach  Number 
cere  is  less  than  2  dB.  In  thu  plane  of  the  rotor  there  is  no  error  since  the  thrust  term 
is  zero.  The  average  error  is  considerably  less  than  1  dB.  Thus,  it  is  concluded  that 
the  concept  of  an  "Effective  Mach  Number"  takes  adequate  account  of  forward  speed 
effects,  at  least  for  typical  helicopter  velocities. 


It  has  been  shown  in  Sections  4.5  and  6  of  the  report  that,  using  a  randomized 
phase  assumption,  and  based  on  experimental  airload  data,  the  thrust  loading  harmonic 
amplitude  can  be  represented  for  the  purpose  of  sound  calculation  by  the  relation¬ 
ship 


Furthermore,  using  these  assumptions,  a  single  loading  point  gives  accurate  results. 
If  we  assume  that  the  ratio  thrust  :  drag  :  radial  component  =  10:1 :1 ,  which  is 
typical  of  helicopter  rotors,  we  can  write  Equation (84) in  the  form 


.  oo 

Cn  *  ^  j(10  nM  sin  0)  J.'  -  J,'  +  (nM  cos  0)  J.' 

X=0 Rr*  5  [ 


(87) 


Firally,  the  sound  intensity  of  the  nth  harmonic  can  be  represented  by  a  relationship 
of  the  form 


r*  A 


f(mB,  M,  0) 


(88) 


where  A  is  the  rotor  disc  area  (nRJ),  and  the  function  f  can  be  inferred  from 
Equation  (87). 


2 

The  sound  pressure  level  corresponding  to  Cn  has  been  computed  for  a  wide  range  of 

mB,  M,  and  0  for  a  rotor  having  nominal  dimensions  (T  -  10,000  lb,  radius  =  20  ft, 

r  =  1000  ft,  and  with  the  loads  acting  at  the  80-percent  radius  point).  The  results 

arc  presented  in  polar  chart  form  in  Figure  53.  The  charts  thus  effectively  give 

the  form  of  the  function  f  in  Equation  (88)  above,  and  they  allow  straightforward 

calculation  of  C2  for  anv  specified  condition. 
n 

Parameters  Required 

The  following  parameters  are  required  for  use  in  the  noise  calculations  using  the 
design  charts. 
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x,  y,  z 


A 

Q 

V 

ao 

Jd 

m 


T 

R 


Field  point  coordinates  relative  to  helicopter  measured  in 
feet  with  x  measured  positive  in  rhe  direction  of  motion 
(parallel  to  ground  in  hover)/  y  measured  sideways  in  the 
plane  of  the  disc,  z  measured  downwards  from  helicopter. 
(Results  for  +y  equal  results  for  -  y.) 

Disc  area,  ft2  (or  T/A  =  disc  loading  in  I b/ ft2) 

Rator  angular  velocity,  rads/sec  (Q  =  rpm  x  2ir/60) 

Flight  velocity,  ft/sec 

Speed  of  sound  in  free  air,  ft/sec 

Disc  incidence  (angle  between  disc  and  x-axis),  deg 

Sound  harmonic  (equals  1  for  fundamental,  2  for  second 
harmonic,  etc.) 

Number  of  blades 

Thrust,  lb 
Rotor  radius,  ft 


Instruction'  for  Use  of  Design  Charts 

To  calculate  the  rotational  noise  spectrum  occurring  instantaneously  at  any  point  r,  0 
relative  to  the  rotor  center  and  its  direction  of  motion 


1)  Calculate  range  r  -  \x2  +  y2  +  z2 

2)  Calculate  the  rotational  Mach  Number  M 

M  =  0.8  — 


3)  Calculate  the  flight  Mach  number 

Mp  =  V/aQ 

4)  Calculate  the  angle  8!  between  the  flight  direction  and  the  line  joining 

he  rotor  and  the  field  point 


0'  -  cos' 1  (x/r) 
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5)  Calculate  the  Effective  Rotational  Mach  Number 


ME  l  -  Mp  cos  0' 


6)  Calculate  the  angle  0  between  the  rotor  plane  and  the  line  r.  If  the  disc 
incidence  is  i^,  this  is  given  by 


9  ’ ,an''  fc M‘!d  CfM 


V  yjx1  +  y2/  y  v 

7)  Using  the  values  of  Mg  ancj  0,  look  up  each  chart  to  obtain  values  of  the 
harmonic  sound  pressure  level  I  for  n  =  2,  3,  4,  6,  8,  10,  12,  16, 

20,  30,, 40,  and  60-. 

8)  Correct  the  values  obtained  for  thrust,  disc,  toading,and  distance  according 
to 


SPLn=  ln  +  11  +  10logl0 


*(*)] 


dB  re:  0.0002  dyne/cm3 


9)  Plot  SPLn  against  n  and  fit  smooth  cu«ve. 

10)  The  sound  nressure  levels  from  this  curve  for  n  =  B,  2B,  3B,  . . .  give  the 
required  harmonic  level  at  the  point  x,y,z. 

1 1)  The  fundamental  frequency  is  QB/^2ir  (1  -  Mp  cos  0))  Hz. 

Example  —  Calculate  the  rotational  noise  spectrum  1000  ft  from  a  three-blade  rotor  at 
an  angle  of  20  deg  below  the  flight  path  for  the  following  parameters:  T  =’  10,000  lb, 
T/A  -  7  lb/ft2,  V  =  200  ft/sec,  ij  =  5 deg,  QR  =  600  ft/sec,  and  aQ  =  1117  ft/sec 

1)  r  -  1000  ft 

2)  M  =  0.8  x  600/1117  =  0.429 

3)  Mp  =  200/1 1 1 7  =  0 . 1 79 


4)  31  -  20  deg 


5)  Mp 


0.429 

1  -  . 1 79  x  .938 


0.516 


6)  9  =  20°  -  5°  -  15° 
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C&mporison  with  experiment 


To  demonstrate  the  usefulness  ond  accuracy  of  this 'simpl  ified  method  spme  ccse.bav*v- 
been  calculated  to  compare  with  the  measured  data  from  Reference  12,  figures  54-— 
end  55  show  the  results.  For  comparison,  the  Theoretical  results  obtained  by  Schlepe! 
ef  cl using  a  computer  program,  are  also  included .  Lt  can  be  seep  that  the  simple  ~ 
technique  produces  results  which  on  the  avercge  are  more  accurate  than  the  computed 
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,1c  Levels,  Results  From  Stuckey  and  Goddatd 


Microphone 


Figure  8.  UH-1A  External  Noise  Spectrum.  Derived  From  Reference 
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Figure  11  .  Contribution  of  Helicopter  Noise  Components  Via  Various 
Methods  of  Data  Presentation. 
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SPL  -  dB  re:  Measured  Overall  Level  SPL  -  dB  re:  Overall  Level 


Oi 


Octave  Band  Center  Frequency  -  Hz 


c)  Octave  Band  Levels. 


Octave  Band  Center  Frequency  -  Hz. 

d)  Subjectively  Corrected  Octave  Band  Levels. 

Figure  11.  (Continued). 
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Figure  22.  Rotcr  Hormoric  Loadings  as  a  Function  of  Spanwise  Station. 
Data  From  Scheimari.2l 
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Figure  23.  Rotor  Harmonic  Loadings  at  0.85  Span  in  Several  Rough  Running 
Cases  for  H-34  Helicopter .  Data  From  Scheiman. 21 


Direction  of  Flight 


arded  Position"  of  Four-Blade  Kigid  Rotor  Showing  Blade 
ent  Location  at  Retarded  Times  for  Sound  ot  Point  50  Feet 
ont  of  Rotor . 
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BcrmOn?  c  Number 


Figure  28.  Fourier  Coefficients  of  Various  Pulse  Shape* 


Sound  Pressure  Level  -  d3  (Arbitrary  Datum) 
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Loading  rtam>or.ic  N<xnt>*«  a 

Figure  30.  Acoustic  Contribution  of  Looking  Horrnonics  10°  From  Sotor  Axis. 
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Figure  34.  Contribution  of  Force  Component!  to  Fourth  Harmonic 
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Sound  Pressure  Level  -  dB  (Arbitrary  Datum) 
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Sound  Harmonic  Number  m 

t- ; gyre  36.  Effect  of  Rotational  Mach  Number  M  on  Sound  Radiation 
10°  From  Rotor  Axis  in  Hover. 
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Sound  Pressor*  t*v*»  *  dB  (Arbitrary  Datum)  So^nd  Pi 


o)  Z*rt»th  Power  Law  f  o*  Loading  Harmonics. 


c)  In Square  Power  Low  for  Loading  Harmonics, 


d)  Invsns  Cut'*  Pwrsr  Low  for  Loodlrtj  Harmonics. 


Figure  37. 


Effect  of  Forward  Velocity  (Mp  =  0.125)  on  Sound  Radiation 
10°  Below  Rotor  Disc.  Rotational  Mach  Number  M  =  0.5. 
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-  2  4  6  8  10  20 

Sound  Hormonfc  Order  m 

e)  invert#  $<yjor»  Power  Low  for  Loodlny  Hormonict. 


1  2  4  6  8  10  20 

Sound  Hormonic  Order  m 

d)  Invert#  Cub#  Powor  Low  for  Looding  Mormon  let. 


gure  38.  Effect  of  Forward  Velocity  (Mp  =  0.25)  on  Sound  Radiation 
10°  Below  Rotor  Disc.  Rotational  Mach  Number  M  =  0.5. 


Sound  Pre»*jre  Level  -  dB  (Arbitrary  Dotum)  $<X>nd  Pre**u'e  Level  -  dB  (Arbitrary  Datum) 


b)  trrvene  FI  rtf  Power  Law  lor  Looijiftfl  Horwoolct. 


Figure  39.  Effect  of  Forward  Veiocity  (Mp  -  0.125)  on  Sound  Radiation 
10°  Below  Rotor  Disc.  Rotational  Mach  Number  M  =  0.75. 
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Loading  Harmonic  X 

Figure  41 .  Accuracy  of  Effective  Velocity  Approximation 


Distance  -  dia 

Figure  46.  Near-Held  Effects  10  Degrees  Below  Rotor  Disc 
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Figure  47.  Effect  of  loading  Distribution  on  Sound  Field 
for  Randomized  Load; . 
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^  |  Measured  Data  (Reference  12) 

- —-  Reference  12  Theory 

- Present  Theory 

-  Presen*  Theo-y.  230  ft  to  Pert 


250  ft  to  Starboard 


Helicopter  Position  With  Respect  to  Observer  -  ft 
a)  First  Harmonic 


b)  Second  Harmonic 

I  jure  48.  Comparison  of  Theoretical  and  Measured  Harmonic  Sound 
Pressure  Levels  for  H-34  Main  Rotor.  Aircraft  Altitude 
200  ft.  Data  From  Reference  12. 
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Figure  50.  Comparison  of  Measured  Internal  and  External  Harmonic  Sound 
Levels  for  CH-47B  Chinook.  Data  From  Reference  7. 


Recorded  Over  Grasi.  Octa  From  Reference  7 
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tt.  SUMF.EMINTARV  NOTE* 

11  ABSTRACT 

1*.  BRONBORINO  MIUTARr  ACTIVITY 

U.S.  Army  Aviation  Materiel  Laboratories 

Fort  EustJs,  Virginia 

A  co.  prehensive  study  of  the  problem  of  helicopter  noise  radiation  is  presented.  After  a  review  of 
the  basic  features  of  the  noise,  the. limited  experimental  data  are  reviewed  in  some  detail,  and 
empirical  laws  are  proposed.  An  exact  theoretical  expression  for  the  noise  is  derived.  This 
expression  has  been  used  as  the  basis  for  the  development  of  a  comprehensive  computer  program  to 
calculate  helicopter  noise  at  any  field  point,  including  all  effects  of  fluctuating  airloads  and  all 
possible  rigid  and  flexible  blade  motions.  Under  very  reasonable  approximations,  an  analytic 
expression  has  been  found  for  the  sound  field  far  from  the  helicopter,  and  computations  based  on 
this  expression  hcs-e  been  made.  The  results  show  all  the  higher  harmonics  of  the  loading,  at  least 
up  to  the  sixtieth,  to  be  significant  for  noise  generation.  A  study  of  the  harmonic  airloads  is 
presented.  Comprehensive  acoustic  results  from  the  theory  include  the  effect  of  various  loading 
inputs,  thrust,  tip  velocity,  number  of  blades,  blade  motion,  and  forward  speed.  Fair  agreement 
with  experiment  is  found  for  overall  levels,  and  good  agreement  Is  found  with  experimental  trends. 
Design  charts  are  presented  which  enable  routine  calculation  to  be  made  of  noise  radiated  from  any 
helicopter  in  hover  or  forward  flight,  i 
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